Abstract
A stressor-related disorder wherein traumatic experience precipitates protracted disruptions to mood and cognition, post-traumatic stress disorder (PTSD) is associated with wide-ranging abnormalities across the body. While various methods have investigated these deviations, only proton magnetic resonance spectroscopy (1H MRS) enables noninvasive measurement of small-molecule metabolites in the living human. 1H MRS has correspondingly been employed to test hypotheses about the composition and function of multiple brain regions putatively involved in PTSD. Here we systematically review methodological considerations and reported findings, both positive and negative, of the current 1H-MRS literature in PTSD (N = 32 studies) to communicate the brain regional metabolite alterations heretofore observed, providing random-effects model meta-analyses for those most extensively studied. Our review suggests significant PTSD-associated decreases in N-acetyl aspartate in bilateral hippocampus and anterior cingulate cortex with less evident effect in other metabolites and regions. Model heterogeneities diverged widely by analysis (I2 < 0.01% to 90.1%) and suggested regional dependence on quantification reference (creatine or otherwise). While observed variabilities in methods and reported findings suggest that 1H-MRS explorations of PTSD could benefit from methodological standardization, informing this standardization by quantitative assessment of the existing literature is currently hampered by its small size and limited scope.
Keywords
Introduction
In post-traumatic stress disorder (PTSD), one or more exposures to a strong stressor precedes chronic disruption to sensory, emotional, and cognitive function. As currently defined in the DSM-V, PTSD takes on four clusters of symptoms that persist for at least one month following either direct or extensive indirect exposure to a traumatic event, including “actual or threatened death, serious injury, or sexual violence”: intrusive memories, dreams, sensory flashbacks, or triggered associations to psychological distress or physiological reactivity; avoidance of internal or external reminders of the trauma; negatively altered cognition and mood; and significantly altered arousal and reactivity. 1
A 2017 survey over 24 countries reported trauma exposure in 70.4% of respondents, with conditional PTSD risk from 0.1% to 19.0%, depending on trauma type. 2 Because trauma exposure thus does not necessitate PTSD development, a persistent question in especially post-trauma case-control analyses thereof is the degree to which group differences reflect PTSD predisposition as opposed to downstream effects of either the trauma itself or the ensuing illness. Studies involving both trauma-exposed individuals with PTSD and their genetically identical trauma-unmatched twins have been useful in unraveling these distinctions in putative causality. 3 For example, hippocampal volume reductions observed in PTSD patients have also been seen in their trauma-unexposed twins, suggesting that a smaller hippocampus may predispose one to develop PTSD post-trauma.4,5 On the other hand, reduced extinction learning retention relative to trauma-exposed controls has been observed in PTSD patients but not their trauma-unmatched twins, suggesting that this factor does not reflect a genetically determined illness predisposition. 6 In the absence of identical twin cohorts, some studies also attempt to separate the influences of trauma itself from those of either predisposition to or downstream effect of PTSD by distinguishing between trauma-unmatched and trauma-exposed controls. Of course, predispositions to and effects of trauma and/or PTSD are not necessarily mutually exclusive phenomena. 3 For example, though reduced hippocampal volume may be a PTSD risk factor,4,5 even trauma-exposed individuals who did not develop PTSD have exhibited reduced hippocampal volumes relative to trauma-unmatched controls, suggesting an additional role for trauma itself in this phenotype. 7
Whether cause and/or effect of illness, multitudinous differences from the molecular to the behavioral have been observed between PTSD patients and various controls (Figure 1). The following introduction is intended only as a brief integrative summary of PTSD's biological context, and the reader is referred elsewhere for comprehensive reviews of PTSD 8 and its genetic associations, 9 transcriptionally visible sex differences, 10 interaction with brain development, 11 pharmacologically manipulable manifestations, 12 functional correlates in the brain, 13 and biological studies in both human and animal models. 3
Twin studies support PTSD heritability estimates ranging from 24% to 71% depending on cohort, with women suggesting higher heritability than men, and on trauma type, with PTSD following interpersonal trauma demonstrating higher heritability than that after accidents or natural disasters. 9 Multiple disease-associated loci encompass both intergenic and potentially protein-coding sequences, supporting a polygenic origin for predisposition 9 possibly modulated by developmental environment. 14 PTSD-associated abnormalities in gene expression15,16 as well as methylation17,18 have also been observed.
Phenotypically, individuals with PTSD have demonstrated differential reactivity to both trauma-related and unrelated stimuli, 3 including enhanced heart rate, skin conductance,19,20 ocular electromyographic responses, 19 diastolic blood pressure, self-reported subjective distress, 20 and cortical event-related potentials 21 relative to various controls.
PTSD-associated changes in arousal state may also manifest in assorted endocrine abnormalities, including smaller pituitary gland volume; 22 reduced plasma cortisol, prolactin, thyroid-stimulating hormone, 23 and GABA 12 but higher cerebrospinal fluid (CSF) cortisol; 24 higher lymphocytic glucocorticoid receptor density;3,25 higher baseline plasma corticotropin-releasing hormone (CRH); 26 lower baseline plasma adrenocorticotropic hormone (ACTH) but increased ACTH following cognitive stress 27 and reduced ACTH following dexamethasone injection;12,28 heightened ACTH release in a patient subset following exposure to opioid receptor antagonist naloxone, thought to disinhibit hypothalamic CRH from facilitating ACTH secretion; 29 lower baseline plasma 30 and CSF 31 neuropeptide Y (NPY), a putative sympathetic nervous system modulator and endogenous anxiolytic; blunted plasma NPY 30 and orbitofrontal and prefrontal cortex glucose metabolism rate response to α2-adrenergic receptor antagonist yohimbine; 32 decreased plasma and CSF hypothalamic neuropeptide orexin-A;12,33 higher metabolic androgen precursor and GABAA receptor (GABAAR) antagonist dehydroepiandrosterone (DHEA); 3 reduced CSF GABAAR modulators allopregnanolone and pregnanolone in women; 34 increased arginine vasopressin (AVP), a hormone possibly potentiating ACTH release, in men; 35 and higher whole-brain cannabinoid receptor type 1 (CB1R) density. 36 While different patient subsets may manifest abnormality in disparate neurotransmitter systems, 37 many of these measures, including plasma cortisol, 23 baseline plasma NPY, 30 lymphocytic glucocorticoid receptor density,3,25 and CSF orexin-A,12,33 have correlated with PTSD symptom severity.
PTSD-associated aberrations in autonomic homeostasis are further accompanied by observations of anatomical, functional, and metabolic abnormalities in multiple brain regions, including those targeted by the in vivo proton magnetic resonance spectroscopy (1H-MRS) experiments forming the basis of the following review: frontal cortex including anterior cingulate cortex, basal ganglia, occipital cortex, parietal cortex and posterior cingulate gyrus, insula, temporal cortex, hippocampus, and amygdala (Supplementary Information).
Despite the unequivocal involvement of manifold biological systems in PTSD predisposition and/or evolution, still no definitive biomarker, MR-visible or otherwise, exists for its diagnosis, prognosis, or treatment monitoring. One result thereof is continued diagnostic uncertainty: Since PTSD entered the American Psychiatric Association Diagnostic and Statistical Manual (DSM) third edition in 1980, 38 the definitions surrounding qualifying traumas and symptoms have changed, most notably between the DSM-IV and DSM-V, resulting in diagnostic discrepancies in as many as 30% of patients in at least one analysis.8,39 Another consequence is limited pharmacotherapeutic options: Only two medications, sertraline and paroxetine, both selective serotonin reuptake inhibitors (SSRIs), are currently FDA-approved for the condition, 12 though research continues on other candidates, including SSRI fluoxetine; serotonin-norepinephrine reuptake inhibitor (SNRI) venlafaxine; dopamine and serotonin antagonist risperidone; multimodal norepinephrine modulator quetiapine; GABA and glutamate modulator topiramate; 40 glutamate modulators lamotrigine and riluzole; N-methyl-D-aspartate receptor (NMDAR) antagonist ketamine; and NMDAR agonist D-cycloserine. 41
Especially while rodent models of PTSD exhibit limited correspondence to human illness as defined by the DSM-V, particularly in female rodents 42 despite the greater published prevalence of PTSD in women than men,43,44 as well as to specialized human prefrontal brain regions implicated in the disease, 3 currently the only means of directly and noninvasively investigating the metabolic manifestation of PTSD is 1H MRS. The present study therefore serves as an update and extension to two previously published surveys of the 1H-MRS literature in PTSD, including a 2010 meta-analysis of findings in hippocampus and anterior cingulate cortex 45 as well as a 2019 systematic review without meta-analysis of 24 research reports published between 1998 and 2017. 46 The present report includes a systematic review of all metabolite comparisons with control and correlations published in located investigations on PTSD, as well as meta-analyses in all regions and metabolites involving at least four reported case-control comparisons therein.
Special attention is accorded in this review to study design decisions including participant sex, age, and comorbidity, as well as other methodological details like scanner field strength, voxel placement, and metabolite focus. Alongside age, sex appears to be an especially important mediator of PTSD development following trauma exposure and is therefore a salient consideration for investigations into possible mechanisms thereof. Lifetime prevalence of PTSD may be twice as high in women as men,42,43 potentially due to greater lifetime exposure to trauma, though disputed; 47 greater exposure to trauma types, like sexual violence, disproportionately likely to precede PSTD development;2,9,48 and/or differential rates of conversion from trauma exposure to PSTD even once trauma frequency and type are controlled.8,49
Comorbidities, including to mood disorders and substance use, are also important to consider in PTSD research. The prevalence of, for example, major depressive disorder (MDD) in PTSD patients may be as high as 50%. 50 These and other comorbidities may influence putative disease findings: While increased plasma cortisol has generally not been observed in PTSD, for example, this may not apply to individuals with comorbid MDD; 3 inversely, reductions in anterior cingulate cortex 5-hydroxytryptamine (5-HT) binding in PTSD patients relative to trauma-unmatched controls may reflect depression. 51 Similarly, comorbidities like alcohol abuse 52 and borderline personality disorder 53 have, like PTSD, themselves been associated with reduced hippocampal volume, 54 with variable confound by the former on the latter.3,55
With such methodological features under explicit consideration, the following review attempts to systematically investigate and report the current state of knowledge regarding 1H-MRS-visible in vivo brain metabolic abnormalities associated with a diagnosis of post-traumatic stress disorder.
Methods
Literature review was conducted in November 2020 and updated in April 2021 according to search terms described previously, 46 including all combinations of “post-traumatic stress disorder” or “PTSD” with “magnetic resonance spectroscopy” or “MRS” on PubMed®. This search yielded thirty-two unique publications that applied magnetic resonance spectroscopy to at least one individual reported to have post-traumatic stress disorder; this set was manually cross-checked with the previously mentioned 2019 systematic literature review 46 and 2010 meta-analysis 45 to ensure that it encompassed all studies included in these two latter reports. Studies were manually excluded if they did not report experiments (N = 22), did not use MR spectroscopy (N = 21) or examine PTSD patients (N = 16), humans (N = 11), or include a matched non-PTSD control group (N = 4), or if they involved a non-psychiatric comorbidity (N = 1) (Figure 2). One exception to these exclusion criteria was a case study included for examination of methodological variables. 56 Review protocols were not registered in any database for public availability.

Post-traumatic stress disorder has been associated with observed disparities from control over multiple biological levels of detail. Post-traumatic stress disorder (PTSD) has been associated with allelic variations at multiple genetic loci; DNA methylation patterns and gene transcription patterns have also been reported to differ with PTSD status and/or stress exposure. Disparate levels of multiple hormones and other biochemicals have been observed in saliva, blood, or cerebrospinal fluid (CSF), and differences in autonomic reactivity to either trauma-related or -unrelated stimuli have been suggested via measures of electromyography (EMG), skin conductance, and heart rate. Brain regions additionally examined by proton magnetic resonance spectroscopy (1H MRS), such as the prefrontal cortex, hippocampus, and amygdala, have demonstrated abnormalities in volume, neurotransmitter receptor binding, and patterns of blood oxygen-level-dependent (BOLD) contrast in PTSD. In assessing these reported differences, as in 1H-MRS findings, it should be noted that depending on study design, biological differences observed between PTSD patients and control may elucidate predisposition to develop PTSD, effect of trauma exposure or PTSD pathology, or all of the above. TLS: total life stress; CRH: corticotrophin-releasing hormone; ACTH: adrenocorticotropic hormone; TSH: thyroid-stimulating hormone; AVP: arginine vasopressin; AlloP: allopregnanolone; CSF: cerebrospinal fluid; EMG: electromyogram; rCBF: regional cerebral blood flow; ACC: anterior cingulate cortex; 5-HT: 5-hydroxytryptamine (serotonin); CB1R: cannabinoid receptor type 1; rACC: rostral anterior cingulate cortex; dACC: dorsal anterior cingulate cortex. This figure was created with BioRender.com.
Data review was conducted manually by a single investigator. For completeness all reported significant and null findings in metabolite comparisons with control are described; for simplicity only significant correlations are reported. Total N-acetyl aspartate, total creatine, and total choline are named by single metabolites for simplicity.
Random-effects meta-analysis 57 for the influence of PTSD on regional metabolite concentration was conducted for all metabolites and brain regions exhibiting at least four published comparisons against either trauma-exposed or trauma-unmatched groups, with or without comorbidities. When PTSD groups were compared with both trauma-unmatched and trauma-exposed controls in the same study, only comparisons with trauma-exposed controls were included for meta-analysis. When PTSD groups with and without additional comorbidities (eg, alcohol use disorder) were compared against distinct matched controls, both comparisons were included in meta-analysis. All statistics were conducted in R version 4 (R Foundation for Statistical Computing, Vienna, Austria).
Results
Methodological Overview
N = 33 PTSD versus control comparisons from N = 28 studies were assessed for average age of the PTSD group as well as age-matching between PTSD and control. Three studies did not provide patient group age means58–60 and one was a case study with no control; 56 these were excluded from age analysis. PTSD group average ages most frequently represented the 30 to 40 year range and demonstrated a marginal trend to increase relative to control (random-effects model standardized mean difference 0.10 ± 0.062, P = .10) (Figure 3A; Supplementary Tables 1 and 2). While N = 10 studies not including the case study on a single male patient 56 included all-male PTSD cohorts, N = 3 studies included all-female PTSD cohorts, with a median 2% decrease in control versus PTSD percent women among those N = 13 studies that did not exhibit sex-matched groups (Figure 3B; Supplementary Tables 3 and 4). N = 4 of 32 studies acquired data at 4 Tesla field strength, while none had been conducted at 7 Tesla or above (Figure 3C; Supplementary Table 5). N-acetyl aspartate (N = 93), choline (N = 48), and creatine (N = 35) comprised the majority of metabolite concentrations reported, particularly in non-cortical voxels, with fewer comparisons involving glutamate and/or glutamine (N = 34), myoinositol (N = 26), GABA (N = 18), glutathione (N = 2), or other resonances including fucoses and lipids (N = 5) (Figure 3D).

Records processing pipeline for systematic review of published proton magnetic resonance spectroscopy (1H-MRS) investigations of post-traumatic stress disorder (PTSD). Initial records (N = 107) were located via designated search terms on PubMed based on Quadrelli et al 46 following which abstracts were screened for exclusion on the basis of not being a human 1H-MRS case study of PTSD or comparative analysis with at least one healthy non-PTSD control. Remaining studies (N = 32) were cross-checked against an older systematic review 46 and meta-analysis 45 to confirm inclusion of all relevant studies in the present investigation. SARS: Severe Acute Respiratory Syndrome.
Explicitly reported consideration of comorbidities was also evaluated, including exclusion or other controls for comorbid psychiatric conditions (Table 1) and past or current use of alcohol, nicotine, or other substances (Table 2) in participants with PTSD.
1H MRS in PTSD Study Design Consideration of Psychiatric Comorbidities.
1H MRS: proton magnetic resonance spectroscopy; PTSD: post-traumatic stress disorder.
1H MRS in PTSD Study Design Consideration of Comorbid Alcohol, Nicotine, or other Substance Use.
1H MRS: proton magnetic resonance spectroscopy; PTSD: post-traumatic stress disorder.
Frontal Cortex (Non-ACC)
In humans the frontal cortex comprises the six cortical layers of the brain's frontal lobe and is thought to be involved in numerous executive processes integrating and mediating activity in other brain regions. Prefrontal cortex (PFC) regions in particular have been implicated in PTSD for their connections with the hippocampus and amygdala that together regulate fear learning and memory, 13 and it has been hypothesized that glutamatergic excitotoxicity in this area may play a role in PTSD pathogenesis. 61
Complementarily to a number of non-1H-MRS findings in the frontal cortex, in PTSD multiple examinations of non-anterior cingulate prefrontal cortex N-acetyl aspartate : creatine (NAA/Cr),61,68 choline : creatine (Cho/Cr)

Overview of current literature on 1H MRS in PTSD by study design, hardware, and reported significant or null findings by metabolite and brain region. Currently published reports on 1H MRS in PTSD cover a broad range of patient ages, with a smaller relative number of investigations into pediatric or geriatric participants, and generally demonstrate small disparities in age between groups under study for analyses reporting experimental group ages (random-effects standardized mean difference SMD 0.10, 95% confidence interval CI −0.02-0.23, P = .08; N = 33 comparisons from 28 studies) (A). The number of 1H-MRS analyses examining male-only participants is higher than the reverse, with an additional skew toward more male-saturated control groups in the few studies that demonstrated uneven between-group sex-matching (B). The majority of existing 1H-MRS investigations of PTSD have been conducted at 1.5 and 3T, field strengths that do not typically enable the straightforward separation of signals from neurotransmitter glutamate from those of its metabolic byproduct glutamine (C). Especially in hippocampus, for which voxel placement and size may contribute to challenging spectral quantification of lower-SNR metabolites, the 1H-MRS literature on PTSD is currently dominated by N-acetyl aspartate (NAA), free cholines (Cho), and creatine (Cr), with fewer investigations on neurotransmitters glutamate (Glu) and GABA or endogenous antioxidant glutathione (GSH) despite putative roles in excitatory/inhibitory balance and oxidative stress potentially important to regional anatomical and functional alterations previously observed in PTSD (D). D based on data reported in Tables 3, 5, and 6; see Supplementary Information for data tables underlying A-C. Gln: glutamine; Glx: glutamate + glutamine. TEC: trauma-exposed control; HC: healthy control; VC: veteran control; AUD: alcohol use disorder.

Findings of significant difference from control in 1H-MRS-visible brain metabolites for patients with post-traumatic stress disorder. General decreases in N-acetyl aspartate (NAA) in PTSD relative to trauma-exposed or -unmatched control groups have been observed predominantly in anterior cingulate cortex and hippocampus, while NAA increases in the amygdala have also been reported by one study. Differences from control in choline (Cho), when observed, have tended to be positive, also in anterior cingulate cortex and hippocampus, while decreases in choline have been observed in the medial temporal lobe. Regional effects of PTSD on common quantification reference metabolite creatine (Cr) have varied by region, with decreases reported in hippocampus and occipital cortex and increases in amygdala. Myoinositol (mIns) has been reported to increase in both anterior cingulate cortex and amygdala. Only a few studies have examined excitatory neurotransmitter glutamate to show increases in medial temporal lobe and hippocampus but decreases in both glutamate + glutamine (Glx) and glutamine (Gln) in anterior cingulate cortex. Reported alterations in GABA have similarly varied by region, with decreases reported in insula, medial temporal lobe, and parieto-occipital cortex; increases in dorsolateral prefrontal cortex; and mixed findings in anterior cingulate cortex. Only two published findings have been reported for glutathione (GSH), both increases, in dorsolateral prefrontal and anterior cingulate cortices. Studies of additional metabolites, such as sugars and lipids measurable by 2D J-resolved proton spectroscopy reported in Quadrelli et al, 79 have been limited. *PTSD versus control groups that included trauma-exposed and/or military veteran individuals. GABA: γ-aminobutyric acid; Glu: glutamate; IMI-1: imidazole from histamine, histidine, and homocarnosine.
Controlled In Vivo 1H-MRS Findings for PTSD in Frontal Cortex and Basal Ganglia Metabolites.
Reported as marginally significant trend.
PTSD group also involved AUD + participants.
1H MRS: proton magnetic resonance spectroscopy; PTSD: post-traumatic stress disorder; PFC: prefrontal cortex; ACC: anterior cingulate cortex; HC: trauma-unmatched control; TEC: trauma-exposed control; VC: veteran control; AUD: alcohol use disorder: NAA: N-acetyl aspartate; Cr: creatine; Cho: choline; mIns: myoinositol; Glx: glutamate + glutamine; Glu: glutamate; Gln: glutamine; GABA: γ-aminobutyric acid; GSH: glutathione.
Correlations with In Vivo 1H-MRS-Visible Frontal, Occipital, and Parietal Cortex Metabolites Reported in PTSD.
Reported as marginally significant trend.
1H MRS: proton magnetic resonance spectroscopy; PTSD: post-traumatic stress disorder; ACC: anterior cingulate cortex; TEC: trauma-exposed control; NAA: N-acetyl aspartate; Cr: creatine; Cho: choline; mIns: myoinositol; Glu: glutamate; GABA: γ-aminobutyric acid; CAPS: clinician-administered PTSD scale; IMI-1: imidazole from histamine, histidine, and homocarnosine.
Anterior Cingulate Cortex
The anterior cingulate cortex (ACC) overlaps with portions of both ventromedial and dorsomedial prefrontal cortex, both of which have demonstrated structural and functional abnormality in PTSD. 13 It is thought that with roles in fear conditioning and extinction the dorsal ACC (dACC) in particular could contribute to PTSD vulnerability and/or pathogenesis, 69 with local functional abnormalities possibly contributing to disruptions in the salience network leading to hypervigilance. 88
In PTSD comparisons against control and/or correlation analyses of anterior cingulate NAA/Cr,58,60–62,70–73 NAA,63,69,71,74–77 creatine,63,69,74–77 Cho/Cr,58,60,69,70,72–77 myoinositol : creatine (mIns/Cr),63,69,72–77 (glutamate + glutamine) : creatine (Glx/Cr), 73 Glx,63,69,74,76,77 glutamate,76,77 glutamine, 69 Glu/Gln,61,69 GABA/Cr,61,62 GABA,69,76,77 GABA/Gln, 69 and GSH/Cr 61 have been reported (Figures 4 and 5; Tables 3 and 4; Supplementary Information).
Basal Ganglia
The basal ganglia are a cluster of subcortical grey-matter nuclei comprising the striatum (caudate and putamen) and the globus pallidus. 89
One previous proton spectroscopy study has targeted this region due to its proximity to the temporal lobes, also linked with a number of PTSD-associated metabolic abnormalities (below) while maintaining greater separation from the hippocampus enabling larger voxel size without partial-volume contamination from this other grey-matter region, as well as for its putative vulnerability to hypoxic damage and exposure to toxic metals. 68 In this study findings against control from left basal ganglia NAA/Cr 68 and Cho/Cr 68 were reported (Figures 4 and 5; Tables 3 and 4; Supplementary Information).
Occipital Cortex
The occipital cortex is a large posterior cortical region thought to predominantly comprise circuitry responsible for visual processing. While PTSD-specific hypotheses regarding this brain area are sparse and it has therefore served as a control region for hypothesis-driven research,72,78 metabolite concentration abnormalities in parieto-occipital regions have been reported in anxiety and mood disorders, also impelling their targeted exploration in PTSD.76,77 In addition, the occipital cortex tends to lie close to radiofrequency receive coils and does not exhibit the steep transitions in magnetic susceptibility between tissue and air cavities that can influence data acquired from temporal and especially frontal regions, justifying its 1H-MRS examination for methodological convenience.
In PTSD comparative and/or correlation analyses of occipital or parieto-occipital cortex NAA/Cr,
72
NAA,76–78 creatine,76–78 Cho/Cr
Controlled In Vivo 1H-MRS Findings for PTSD in Occipital and Parietal Cortex Metabolites.
1H MRS: proton magnetic resonance spectroscopy; PTSD: post-traumatic stress disorder; HC: trauma-unmatched control; TEC: trauma-exposed control; AUD: alcohol use disorder: NAA: N-acetyl aspartate; Cr: creatine; Cho: choline; mIns: myoinositol; Glu: glutamate; Gln: glutamine; GABA: γ-aminobutyric acid; IMI-1: imidazole from histamine, histidine, and homocarnosine.
Parietal Cortex
The parietal cortex, posterior to the frontal lobe and superior to the occipital lobe, is a large area of associative cortex involved in somatosensory integration, motor planning, and attention, among other functions. 90 While, like the occipital cortex, the parietal cortex has been implicated in a minimum of hypothesis-driven 1H-MRS research into PTSD, it has been involved in one study of regional periventricular white matter. 68 It also envelops the posterior cingulate gyrus, the traditionally preferred target of 2D COrrelated SpectroscopY (COSY) for its proximity to the radiofrequency receive hardware employed. 79
In PTSD comparative and/or correlation analyses of parietal cortex or posterior cingulate NAA/Cr, 68 NAA, 79 Cho/Cr, 68 glutamate, 79 glutamine, 79 myoinositol, 79 GABA, 79 the imidazole moiety visible on histamine, histidine, and homocarnosine (IMI-1), 79 fucosylated glycans (including t-fucose, fucose IV, and fucose VI), 79 and lipid HC = CH–CH2–CH2–CH3 79 have been reported (Figures 4 and 5; Tables 4 and 5; Supplementary Information).
Insula
The insula correspond to Brodmann Areas 13 through 16, each located deep in the bilateral Sylvian fissures beneath the corresponding opercula, and are connected with several brain regions demonstrating reported associations with PTSD, including the prefrontal cortex, anterior cingulate gyrus, hippocampus, and amygdala. 91 They are also thought to be involved in emotional recall and have, among other regions, demonstrated regional cerebral blood flow signals correlating positively with severity of PTSD-associated flashbacks, 92 justifying their exploration by 1H MRS. 54
In PTSD comparative and/or correlation analyses of right insular NAA/Cr,54,62 NAA, 54 creatine, 54 and GABA/Cr, 62 as well as left insular NAA/Cr, 54 NAA, 54 and creatine 54 have been reported (Figures 4 and 5; Tables 6 and 7; Supplementary Information).
Controlled In Vivo Proton 1H-MRS Findings for PTSD in Insular, Hippocampal, and Amygdalar Metabolites.
Reported as marginally significant trend.
PTSD group also involved AUD + participants.
1H MRS: proton magnetic resonance spectroscopy; PTSD: post-traumatic stress disorder; HC: trauma-unmatched control; TEC: trauma-exposed control; VC: veteran control; AUD: alcohol use disorder: NAA: N-acetyl aspartate; Cr: creatine; Cho: choline; mIns: myoinositol; Glx: glutamate + glutamine; Glu: glutamate; Gln: glutamine; GABA: γ-aminobutyric acid; BPD: borderline personality disorder.
Correlations with In Vivo 1H-MRS-Visible Insular, Hippocampal, and Amygdalar Metabolites Reported in PTSD.
Reported as marginally significant trend.
1H MRS: Proton magnetic resonance spectroscopy; PTSD: Post-traumatic stress disorder; HC: trauma-unmatched control; VC: veteran control; TEC: trauma-exposed control; GABA: γ-aminobutyric acid; NAA: N-acetyl aspartate; Cr: creatine; Cho: choline; mIns: myoinositol; Glu: glutamate; CAPS: clinician-administered PTSD scale.
Temporal Lobe (Right)
The temporal lobe is the seat of both the hippocampus and the amygdala, subcortical structures that have both demonstrated a number of structural and functional abnormalities in PTSD (Supplementary Information).
The subcortical hippocampus within the medial temporal lobe has generally been the target of proton spectroscopy voxels located in the temporal lobe, including within one group choosing the nomenclature “temporal cortex” instead.76,77 In PTSD comparative and/or correlation analyses of right temporal lobe NAA/Cr,64,65,80 NAA,76,77 creatine,76,77 Cho/Cr,65,80 choline,76,77 mIns/Cr, 65 myoinositol,76,77 glutamate,76,77 and GABA76,77 have been reported (Figures 4 and 5; Tables 6 and 7; Supplementary Information).
Temporal Lobe (Left)
In PTSD comparative and/or correlation analyses of left temporal lobe NAA/Cr,64,65,80 Cho/Cr,65,80 and mIns/Cr 65 have been reported (Figures 4 and 5; Tables 6 and 7; Supplementary Information).
Hippocampus (Right)
The hippocampus is a subcortical structure implicated in the encoding of spatial information as well as episodic memories. 93
Complementary to a number of other non-1H-MRS investigations in this region, comparative and/or correlation 1H-MRS analyses of right hippocampal NAA/Cr,54,58,60,66,67,71,81–83 NAA,54,71,75,78,82,84,85 NAA/Cho, 83 creatine,54,75,78,84 Cho/Cr,58,60,66,81 choline,75,78,84 myoinositol, 75 Glu/Cr, 82 Glu/NAA, 82 and glutamate 82 have been reported (Figures 4 and 5; Tables 6 and 7; Supplementary Information).
Hippocampus (Left)
In PTSD comparative and/or correlation analyses of left hippocampal NAA/Cr,54,58,60,66,67,71,81–83 NAA,54,63,71,75,78,82,84,85 NAA/Cho, 83 creatine,54,75,78,84 Cho/Cr,58,60,66,81 choline,75,78,84 myoinositol, 75 Glx, 63 Glu/Cr, 82 Glu/NAA, 82 and glutamate63,82 have been reported (Figures 4 and 5; Tables 6 and 7; Supplementary Information).
Amygdala
The amygdala comprises two almond-shaped subcortical structures located in the bilateral medial temporal lobes and is implicated in the integration of sensory information, predicted threats, and ensuing emotional and behavioral responses. As such, the amygdalae exhibit connections with the hippocampus and multiple regions in the prefrontal cortex, with which they are thought to regulate the learning and memory of fear. 13
Complementarily to numerous non-1H-MRS findings in PTSD, comparative and/or correlation analyses of amygdalar NAA,74,86 creatine,59,74,86 choline,74,86 myoinositol,74,86 and Glx 74 have been reported (Figures 4 and 5; Tables 6 and 7; Supplementary Information).
Meta-Analysis of Anterior Cingulate Cortex and Bilateral Hippocampus Metabolite Findings
We conducted meta-analyses for all regions and metabolites for which at least four previously published comparisons with explicitly reported group statistics were present, which included NAA (seven), NAA/Cr (six), Glx or glutamate (four), creatine (four), myoinositol (four), and Cho/Cr (four) in the anterior cingulate cortex as well as NAA (seven, six) or NAA/Cr (seven, six) in the left and right hippocampus, respectively.
Random-effects models demonstrated significantly lower anterior cingulate cortex NAA in PTSD when referenced to creatine (k = 6 studies; pooled standardized mean difference −0.60, 95% confidence interval (CI) −1.16 to −0.05, P = .04) but not otherwise (k = 7 comparisons in 6 studies). No significant effect of PTSD was discerned for anterior cingulate cortex Glx or glutamate (k = 4 studies), creatine (k = 4 studies), myoinositol (k = 4 studies), or Cho/Cr (k = 4 studies). Between-study heterogeneities appeared to be higher for NAA (NAA/Cr Q = 14.76, P = .01, I2 = 66.1%; NAA Q = 39.77, P < .0001, I2 = 84.9%) reports than those for Glx or glutamate findings (Q = 0.88, P = .8, I2 < 0.01%) or creatine (Q = 3.81, P = .3, I2 = 21.2%) regardless of quantification reference. Because study pools from Meyerhoff et al 76 and Pennington et al 77 exhibited significant overlap, only values from Pennington et al 77 were reported here (Figure 6; Supplementary Tables 6-11).

Findings of significant correlation or regression coefficient between 1H-MRS-measured brain metabolites with each other or other measures in post-traumatic stress disorder (PTSD). Because proton magnetic resonance spectroscopy (1H-MRS) experiments measure multi-functional small molecules in large (>1 cm3) heterogeneous voxels containing many different types of tissue, cell, and associated functional circuitry, correlations between measured metabolite concentrations and anatomical or cognitive measures or clinically relevant variables like symptom reports are useful for interpreting other putative metabolite findings in 1H-MRS studies of disease. The majority of such reported correlations have pertained to clinician-administered PTSD scale (CAPS) indices like arousal, intrusion, symptom re-experience, or total score. Other correlates have included self-reported insomnia indices, magnetic resonance imaging (MRI)-measured brain region volumes, memory performance, and other metabolites, though the latter may also sometimes indicate relationships introduced by quantification method rather than findings of biological significance. *In combined PTSD and control groups including trauma-exposed and/or military veteran individuals; **In combined PTSD and trauma-unmatched control groups. NAA: N-acetyl aspartate; Cr: creatine; mIns: myoinositol; GABA: γ-aminobutyric acid; Glu: glutamate; GSH: glutathione; Cho: total choline; IMI-1: imidazole from histamine, histidine, and homocarnosine.
In addition, random-effects models demonstrated significantly lower left hippocampus N-acetyl aspartate (NAA) in PTSD when referenced to creatine (Cr) (k = 7 studies; pooled standardized mean difference −0.99, 95% CI −1.94 to −0.04, P = .04) with a marginally significant decrease on the right (k = 6 studies; pooled standardized mean difference −0.75, 95% CI −1.54 to 0.04, P = .06), while NAA referenced otherwise demonstrated significant effect for both left (k = 7 studies; pooled standardized mean difference −0.78, 95% CI −1.11 to −0.45, P = .001) and right (k = 6 studies; pooled standardized mean difference −0.73, 95% CI −1.25 to −0.22, P = .01) hippocampus. Between-study heterogeneities appeared to be higher for studies of hippocampal NAA when creatine was used as a reference (left hippocampus NAA/Cr Q = 36.21, P < .0001, I2 = 83.4%; right hippocampus NAA/Cr Q = 22.96, P = .0003, I2 = 78.2%) than when it was not (left hippocampus NAA Q = 6.10, P = .4, I2 = 1.7%; right hippocampus NAA Q = 8.34, P = .1, I2 = 40.0%) (Figure 7; Supplementary Tables 12-15).

Random-effects models of group difference between post-traumatic stress disorder (PTSD) and control in proton magnetic resonance spectroscopy (1H-MRS)-measured anterior cingulate cortex metabolite concentrations examined by four or more independent published investigations. Random-effects models across the listed studies demonstrated significantly lower anterior cingulate cortex N-acetyl aspartate (NAA) in PTSD when referenced to total creatine (Cr) (pooled standardized mean difference −0.60, 95% confidence interval CI −1.16 to −0.05, P = .04) but not otherwise. No significant effect of PTSD was discerned for anterior cingulate cortex glutamate (Glx, Glu), creatine (Cr), myoinositol (mIns), or creatine-referenced choline (Cho/Cr). Between-study heterogeneities were higher for NAA reports than those for glutamate or creatine findings regardless of quantification reference. Because study pools from Meyerhoff et al 76 and Pennington et al 77 exhibited overlap, only values Pennington et al 77 were reported here. When PTSD groups were compared with both trauma-unmatched and trauma-exposed controls (TEC) in the same study, only comparisons with TEC were included in this analysis. Supporting group statistics reported in Supplementary Tables 6-11. AUD: alcohol use disorder.
Discussion
Here we have reviewed the proton magnetic resonance spectroscopy literature targeted to comparisons of individuals with PTSD and either trauma-exposed or trauma-unmatched controls to find N = 32 publications, including 31 studies reporting group 1H-MRS data collection and one case report included for methodological considerations, fitting our inclusion criteria.
While the current 1H-MRS literature in PTSD exhibits better age-matching between experimental and control groups than the currently available 1H-MRS literature of at least one other disease, 94 the distribution of ages previously studied is not uniform (Figure 3A) and the sex distribution is heavily skewed toward male participants (Figure 3B) despite a lifetime PTSD prevalence in women potentially double that of men.42,43 In addition, hardware fields employed by reviewed 1H-MRS analyses of PTSD do not exceed 4 Tesla (Figure 3C), which may contribute to the relative skew of reported metabolite comparisons toward traditionally high-SNR singlets like NAA, creatine, and choline (Figure 3D).

Random-effects models of group difference between PTSD and control in in proton magnetic resonance spectroscopy (1H-MRS)-measured hippocampal metabolite concentrations examined by four or more independent published investigations. Random-effects models across the listed studies demonstrated significantly lower left and a trend to lower right hippocampus N-acetyl aspartate (NAA) in PTSD when referenced to creatine (Cr) (pooled standardized mean difference left −0.99, 95% CI −1.94 to −0.04, P=.04; right −0.75, 95% confidence interval CI −1.54 to 0.04, P=.06); NAA referenced otherwise demonstrates significant effect for left (pooled standardized mean difference −0.78, 95% CI −1.11 to −0.45, P=.001) and right (pooled standardized mean difference −0.73, 95% CI −1.25 to −0.22, P=.01) hippocampus. Between-study heterogeneities were higher for studies of hippocampal NAA when creatine was used as a reference than when it was not. When PTSD groups were compared with both trauma-unmatched and trauma-exposed controls (TEC) in the same study, only comparisons with TEC were included in this analysis. Supporting group statistics reported in Supplementary Tables 12-15. *Descriptive statistics estimated from Figure 2. VC: veteran control; HC: healthy control.
Meta-analyses for all regions and metabolites for which at least four previously published comparisons with control were present included NAA and NAA/Cr in anterior cingulate cortex and bilateral hippocampus as well as Glx or glutamate, creatine, myoinositol, and Cho/Cr in the anterior cingulate cortex; no other metabolites in any other region passed this k = 4 threshold of published studies.
As concluded for a recent similar review of the literature on multiple sclerosis, 94 regional decreases in N-acetyl aspartate are potentially too nonspecific and insensitive to serve as a biomarker for PTSD vulnerability, development, or treatment response. Abnormal decreases in anterior cingulate NAA or NAA/Cr have also been reported by meta-analysis of studies on first-episode psychosis 95 and autism spectrum disorders; 96 reductions in hippocampal NAA or NAA/Cr have additionally been observed in meta-analyses on schizophrenia, 95 fibromyalgia, 97 obstructive sleep apnea, 98 mild cognitive impairment, 99 and Alzheimer's disease. 100 The heavy focus of currently published 1H-MRS literature in PTSD on N-acetyl aspartate therefore in principle limits the application of its findings to the development of clinically useful 1H-MRS-based regional single-molecule biomarkers, calling for further targeted research into a more diverse range of metabolites measurable by the method, especially at higher field strengths and using acquisition, processing, and analysis pipelines that minimize both imprecision and inaccuracy in metabolite concentration estimates.
Meta-analysis demonstrated decreased anterior cingulate cortex NAA in PTSD only when referenced to creatine and decreased hippocampal NAA regardless of reference. These findings further support the conclusions of a previous systematic review observing consistent decreases in hippocampal and anterior cingulate cortex NAA across the 1H-MRS literature on PTSD. 46 In particular, our anterior cingulate cortex results also extend previous meta-analytic findings of decreased anterior cingulate NAA/Cr in PTSD relative to healthy controls and update three previously reviewed comparisons from two studies71,75 underlying a significant effect of NAA in anterior cingulate cortex 45 with four more63,69,74,77 that bring this metabolite result to non-significance. Similarly, reported reductions in left hippocampal NAA/Cr and NAA extend previously significant meta-analytic findings in both indices. 45 Our finding of significant right hippocampal NAA decreases in PTSD repeats in a mixed-control cohort a previous meta-analytic result using trauma-exposed controls, 45 while our own marginally significant random-effects model for right hippocampal NAA/Cr did not replicate a previous significant result, potentially because of our choice to exclude three previously incorporated studies that either did not report at least one group standard deviation or reported a voxel in the medial temporal lobe but not specifically the hippocampus.60,65,71
Conflicting meta-analytic results between NAA and NAA/Cr in anterior cingulate cortex and right hippocampus may be suggestive of systematic heterogeneities between studies choosing to employ either referencing technique, ie, groups choosing either creatine or non-creatine references may differ in other methodological details that themselves influence reported metabolite effect size. Indeed, very few groups included in our meta-analyses reported both values, making this a feasible explanation for the differences observed here. NAA results that differ by concentration reference may also suggest the disparate influence of the reference itself either as an unwanted source of error or, on the other hand, as more sensitive to actual biological differences among the participants of each study. It is important to note that these two directions of influence are not necessarily mutually exclusive and may lead to different outcomes in different brain regions, as we observed here for PTSD-associated reductions in anterior cingulate NAA that were definitively significant only when referenced to creatine but right hippocampus NAA significant only when not referenced to creatine. It is of potential note that increased cross-study heterogeneity was observed among reports of left and right hippocampal NAA/Cr relative to NAA, which, when combined with the lower number of included studies targeting right than left hippocampal NAA/Cr, may have contributed to the lack of significance in the former statistic, especially considering the previously mentioned published significant meta-analytic finding in the same region. 45
Meta-analysis of anterior cingulate glutamate or Glx, not previously attempted by the above reference, 45 demonstrated no significant effect of PTSD. So far only one significant effect of PTSD on anterior cingulate cortex glutamate or Glx, referenced to creatine, has been published, 73 though significant negative correlations with clinician-administered PTSD scale (CAPS) arousal76,77 and total scores 77 have been previously observed. In one model of PTSD pathogenesis, regional glutamate excitotoxicity can result from increases in the ratio of extrasynaptic to synaptic glutamate due to glutamate “spillover,” a consequence of inadequate astrocytic clearance of extrasynaptic glutamate, and attendant stimulation of presynaptic mGluR-2/3 inhibiting its further release at the synapse. 41 It is important to note that while 1H MRS may provide clues regarding the biochemical composition over a whole voxel, eg, the possible cross-voxel density of metabolically healthy neurons 75 or glial cells 74 as inferred from N-acetyl aspartate and myoinositol or choline concentrations, respectively, the proposed mechanism of regional excitotoxic damage in PTSD does not necessitate cross-voxel changes in glutamate concentration large enough to be detected by 1H MRS. It remains therefore important to note that the mechanistic reasons behind disease-related metabolic abnormalities visible to 1H MRS are often themselves unknown, complicating the use of this method for conclusive hypothesis testing regarding pathophysiological mechanisms at the cellular or molecular scale.
This currently small literature limits the feasibility of further meta-analysis evaluating the use of 1H MRS as a source of biomarkers for PTSD risk, effect, or both via parsing effect sizes by control group type (trauma-exposed or trauma-unmatched), though such bifurcations have been previously reported. 45 Even the most well studied metabolite (NAA) in the most widely studied regions (anterior cingulate cortex and left hippocampus) represents a maximum of four published analyses of either trauma-exposed or trauma-unmatched controls, each with group sizes well under N = 100.
In addition, psychiatric comorbidity may present a limiting factor for assigning specificity of reported metabolic changes to PTSD alone. For example, shared metabolic effects have been reported between PTSD and anxiety, 101 or abnormalities reported in PTSD groups with comorbid depressive symptoms, 51 the latter of which has a documented association with PTSD. 50 Changes in brain morphology52,55 and neurometabolite concentrations have also been found with both acute 102 and chronic 103 alcohol use, as in tobacco or nicotine use, 104 also possibly confounding putative findings attributed to PTSD status. As such, parsing the specific effects of PTSD on metabolite concentrations apart from those associated with comorbid conditions remains a significant challenge.
A few investigations reviewed here noted associations between certain comorbid symptom axes and the metabolites under study, eg, parieto-occipital GABA and depressive symptoms 76 or right insula GABA/Cr and anxiety symptoms, 62 or showed metabolite differences between PTSD patients with and without, for example, alcohol use disorder. 77 Others statistically controlled their results for comorbidities like nicotine use to find no interactive effect.61,76 Despite such efforts, it is important to note the small number of participants within these studies and attendant limitations in statistical power and generalizability.
Taken together, our findings suggest a small (N = 32 studies) literature using 1H MRS to investigate PTSD that is well balanced in terms of participant ages but heavily skewed toward the male sex, lower field strengths, limited regional targets, and high-SNR metabolites. Regardless, a number of notable findings have been reported, including in brain areas like dorsolateral prefrontal cortex, anterior cingulate cortex, amygdala, and hippocampus previously shown to exhibit structural and/or functional abnormalities in PTSD relative to either trauma-unmatched or trauma-exposed controls (Figure 4, Tables 3, 5, and 6; Supplementary Information), as well as a number of significant correlations between metabolites in these regions and clinically relevant variables like CAPS scores or other measures of symptom experience (Figure 5, Tables 4 and 7; Supplementary Information). As such, the potential for 1H MRS to forge new ground in the understanding of PTSD, as well as play a role in the development of biomarkers for risk, diagnosis, and/or treatment monitoring, is not only as-yet largely untried but has also yielded meaningfully interpretable information where it has been, therefore remaining a useful piece of the puzzle behind the pathogenesis and evolution of this life-changing condition.
Supplemental Material
sj-docx-1-css-10.1177_24705470221128004 - Supplemental material for Proton Magnetic Resonance Spectroscopy in Post-Traumatic Stress Disorder—Updated Systematic Review and Meta-Analysis
Supplemental material, sj-docx-1-css-10.1177_24705470221128004 for Proton Magnetic Resonance Spectroscopy in Post-Traumatic Stress Disorder—Updated Systematic Review and Meta-Analysis by Kelley M. Swanberg, Leonardo Campos, Chadi G. Abdallah and Christoph Juchem in Chronic Stress
Footnotes
Author Contributions
KMS: Data review and analysis, manuscript drafting; LC: Data review, manuscript drafting; CGA: Manuscript editing; CJ: Manuscript editing.
Declaration of Conflicting Interests
CGA has served as a consultant and/or on advisory boards for Aptinyx, Genentech, Janssen, Psilocybin Labs, Lundbeck, Guidepoint, and FSV7, and as editor of Chronic Stress for Sage Publications, Inc. He also filed a patent for using mTORC1 inhibitors to augment the effects of antidepressants (Aug 20, 2018).
Ethical Approval
Not applicable, because this article does not contain any studies with human or animal subjects.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by the Yale Center for Clinical Investigation.
Informed Consent
Not applicable, because this article does not contain any studies with human or animal subjects.
Trial Registration
Not applicable, because this article does not contain any clinical trials.
Supplemental Material
Supplemental material for this article is available online.
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
