Deployed military personnel are at risk for development of posttraumatic stress disorder (PTSD). Among U.S. marines and soldiers, PTSD prevalence after deployment to Iraq and Afghanistan is estimated at 6%-13% (1). Development of PTSD is hypothesized to be influenced by preexisting vulnerability factors. Pole et al. (2) observed in police academy cadets that increased skin conductance and decreased skin conductance habituation in response to startling sounds under experimental threat were prospectively related to PTSD symptoms 1 year later. No blood-derived biomarkers representing vulnerability factors for development of PTSD have been identified as of yet.
PTSD appears to be associated with altered activity of the hypothalamic-pituitary-adrenal (HPA) axis (3). However, a meta-analysis revealed that hypocortisolism is not a consistent finding in PTSD (4). In addition, there is evidence for higher glucocorticoid sensitivity of the HPA axis in PTSD (5, 6), but this may be induced by military deployment independently of the disorder (7, 8). Regulation of immune responses by glucocorticoids also seems to be altered in PTSD patients. Our group previously reported lower glucocorticoid sensitivity of T-cell proliferation and higher glucocorticoid sensitivity of lipopolysaccharide-stimulated cytokine production in PTSD patients (9). Accordingly, higher in vitro inhibition of lipopolysaccharide-stimulated production of cytokines (10) and lysozymes (11) by glucocorticoids has been observed in PTSD patients as well. However, at least part of the relation between PTSD and glucocorticoid regulation of the immune response could be attributed to trauma exposure, since trauma-exposed individuals without the disorder also displayed alterations in glucocorticoid sensitivity of the immune system compared with healthy subjects (9, 11).
Glucocorticoid receptors (GRs) and mineralocorticoid receptors mediate glucocorticoid actions. Mineralocorticoid receptors are saturated at low basal cortisol levels, while GRs play a key regulatory role during periods of increased cortisol (e.g., in stressful conditions) (12). GRs reside in the cytosol and, upon ligand binding, translocate to the nucleus to regulate gene transcription (13). GR sensitivity is influenced by heat-shock proteins (HSP70/HSP90) chaperoning GRs, which influence proper folding and maturation of the GR and subsequent translocation and gene transcription (14). FKBP5, both target gene and cochaperone of the GR-HSP70/90 heterocomplex, can lower GR affinity and thereby alter glucocorticoid binding (15).
Cross-sectional studies on the GR number in peripheral mononuclear cells (PBMCs) in individuals with PTSD have yielded mixed results. Higher (16, 17), lower (18), and similar (5, 6, 8, 9, 11) GR numbers in PTSD patients relative to comparison subjects have been reported. Differences in the duration of PTSD, type of trauma, reference groups (e.g., trauma-exposed comparison subjects or healthy comparison subjects), and techniques used to determine the GR number may all contribute to the reported variability in changes in the GR number. In addition, it remains unknown whether preexisting differences in the GR number are present in individuals prone to developing PTSD after a traumatic experience.
We previously reported a higher GR number in PBMCs before deployment in military personnel who had developed depressive symptoms and severe fatigue after deployment (19). This higher GR number was not caused by preexisting differences in mental health status. We hypothesized, therefore, that a high GR number in PBMCs before deployment may represent a vulnerability factor for development of mental health problems after deployment.
In the present study, we investigated, before deployment, whether military personnel with and without deployment-related onset of PTSD symptoms differed in dexamethasone binding capacity of PBMCs. Since the gluocorticoid:mineralocorticoid receptor ratio in human PBMCs is approximately 10:1 (20), we refer to the results of the binding assay as the GR number or binding. Moreover, the predictive value of the preexisting GR number for development of PTSD symptoms after deployment was investigated. Subsequently, we assessed whether group differences existed in mRNA expression for subtypes of the GR (GR-α, GR-β, and GR-P), GR target genes (glucocorticoid-induced leucine zipper [GILZ] and serum and glucocorticoid-inducible kinase-1 [SGK-1]), and the GR target gene and cochaperone FKBP5.
This study is part of a prospective cohort study on deployment-related disorders in the Dutch Armed Forces. Individuals volunteered to participate prior to a 4-month deployment to Afghanistan. Duties during deployment included combat patrols, clearing or searching buildings, participation in demining operations, and transportation across enemy territory. The combat group was exposed to typical war zone stressors, including expo sure to enemy fire, armed combat, and seriously injured and dead fellow soldiers and civilians. The study was approved by the Institutional Review Board of the University Medical Center Utrecht (Utrecht, the Netherlands). Written consent was also obtained. One to 2 months prior to deployment and approximately 1 and 6 months after deployment, participants completed questionnaires and a heparinized blood sample was drawn between 8:00 and 11:30 a.m.
Participants were selected from 455 men who completed all assessments. Analysis of GR binding and associated measures was performed by individuals blind to the participants' PTSD status. Participants were assigned to the PTSD group when their score on the Self-Rating Inventory for PTSD (21) 6 months after deployment was ≥38. This cutoff score corresponds to the mean plus two standard deviations, which coincides with the 95th percentile of scores before deployment within a population of 704 soldiers from the Dutch Armed Forces (mean: 26.91 [SD=5.34]). The validity of this cutoff score is supported by van Zelst et al. (22), who compared Self-Rating Inventory for PTSD scores and ratings from a diagnostic interview within a community population of older adults.
Since we investigated PTSD onset, eight participants with scores ≥38 before deployment were excluded. All other participants fulfilling the selection criteria for the PTSD group were included (N=34). Subsequently, a comparison group (N=34) was identified by matching, to participants in the PTSD group, for age, period of deployment, function and rank during deployment, education, and number of previous deployments. Scores had to be below the cutoff score at all time points. Since we previously found that depressive symptoms are related to high GR binding in PBMCs, participants with high levels of depressive symptoms before or after deployment (Symptom Checklist-90 depression subscale score ≥24 ) were excluded from the comparison group.
Matched comparison subjects were compared with non-selected participants without high levels of PTSD or depressive symptoms before deployment or 6 months after deployment. Matched comparison subjects were somewhat younger (matched comparison subjects: mean age=26.68 years [SD=9.28]; nonselected comparison subjects: mean age=29.42 years [SD=9.06]; t=−1.679, df=380, p=0.09) and had experienced fewer previous deployments (matched comparison subjects: mean=0.65 [SD=0.85]; nonselected comparison subjects: mean=1.06 [SD=1.34]; t=−2.547, df=50.71, p=0.01) than nonselected comparison subjects. There were no group differences in year of deployment, rank and function during deployment, and education.
PTSD symptom level over the past 4 weeks was assessed with the 22-item Self-Report Inventory for PTSD. A higher score indicates more PTSD symptoms (range: 22-88). This inventory has good concurrent validity with other PTSD measures, such as the Clinician Administered PTSD Scale and the Mississippi Scale for PTSD (21). Level of depressive symptoms was assessed using the Dutch version of the 16-item Symptom Checklist-90 depression subscale (23). A higher score on this measure indicates more depressive symptoms (range: 16-80). Exposure to potentially traumatic experiences before the age of 18 was assessed using the Dutch version of the short-form Early Trauma Inventory-Self-Report (24). The questionnaire consists of 27 dichotomous items. The total score represents the number of experienced events. Exposure to potentially traumatic deployment stressors was assessed with a 13-item checklist specifically developed for the present study (see the data supplement accompanying the online version of this article).
For determination of the capacity of PBMCs to bind glucocorticoids, a validated whole cell single-point binding assay was used, as described previously (19). This method provides a reliable estimate of Bmax, as determined using a classical binding assay with 3-200 nM 3H-dexamethasone (r2=0.92) (19). Briefiy, PBMCs were isolated from whole blood using Ficoll-Paque (Pharmacia and Upjohn, Uppsala, Sweden), and 107 cells were frozen in dime-thyl sulfoxide. After thawing and 60 minutes of equilibration in a culture medium, cells were washed twice, resuspended in assay buffer (RPMI-1640, with 5% fluorescence correlation spectroscopy), and incubated in duplicate with 100 nM 3H-dexamethasone (Amersham, Buckinghamshire, United Kingdom) in the presence or absence of excess unlabeled dexamethasone (Sigma-Aldrich, Steinheim, Germany) for 1 hour at 37°C. Cell-bound radioactivity was quantified by liquid scintillation analysis.
Total RNA was isolated from PBMCs with Trizol (Invitrogen, Carlsbad, Calif.). One μg of total RNA was used to synthesize cDNA with SuperScript Reverse Transcriptase (Invitrogen, Carlsbad, Calif.). Real-time polymerase chain reactions were performed with an iQ5 Real-Time PCR Detection System (Bio-Rad Laboratories, Hercules, Calif.) (see Table 1 for primer sequences). Data were normalized for glyceraldehyde-3-phosphate dehydrogenase (GAPDH) and β-actin expression.
Primer Sequences Used in Quantitative Real-Time Polymerase Chain Reaction
| Add to My POL
Plasma cortisol levels were measured by electrochemiluminescence immunoassay on the Modular E170 (Roche Diagnostics, Mannheim, Germany). Lower limit of detection was 3 nmol/l; interassay variation was <3%; and reference values (7:00-10:00 a.m.) were 170-540 nmol/l.
Leukocyte Subset Analysis
Leukocyte subsets were assessed by dual-color fluorescence analysis with a Becton Dickinson Calibur flow cytometer (Becton, Dickinson and Company, Franklin Lakes, N.J.) to quantify cluster of differentiation (CD) 14+ (monocytes), CD3+ (T cells), CD4+ (T helper/inducer), CD8+ (T suppressor/cytotoxic), and CD19+ (B cells), as described previously (19).
Analyses were performed using SPSS, Version 15.0 (SPSS, Inc., Chicago). A p value <0.05 (two-tailed) was considered significant. Variables were tested for normality and log-transformed when necessary. Outliers were removed if z values fell outside the range of ±3.29 (GR-α: comparison group, N=1; GR-β: PTSD group, N=1; GR-P: comparison group, N=1; FKBP5: PTSD group, N=2; GILZ: PTSD group, N=1; SGK-1: PTSD group, N=2). Group differences were tested with t tests and analysis of covariance for continuous parametric variables and with chi-square tests for categorical variables. Mann-Whitney tests were used for nonparametric log-transformed variables. Logistic regression analysis with group as a dependent variable, controlling for potentially traumatic childhood experiences and predeployment questionnaire scores for PTSD and depressive symptoms, was performed to test the predictive value of the preexisting GR number for PTSD status after deployment. This analysis provided the odds ratio for inclusion in the PTSD group with each GR increase of 1,000. Linear regression analysis was performed with the preexisting GR number as a dependent variable to test whether the GR number was associated with predeployment participant characteristics. For regression analyses missing variables for predeployment Symptom Checklist-90 depression sub -scale scores (PTSD group: N=1), the number of previous deployments (PTSD group: N=1), and body mass index (PTSD group: N=3) were imputed in SPSS using the option "linear trend at point."
Participant characteristics are depicted in Table 2. Exposure to potentially traumatic events during deployment was similar for both groups (t=−0.751, df=62, p=0.46). After deployment, the PTSD group had more depressive symptoms than the comparison group (t=−7.572, df=66, p<0.001). There were no significant group differences for the matching variables. Participants in the PTSD group already reported somewhat more PTSD symptoms before deployment than matched comparison subjects (t=−4.284, df=66, p<0.001). Participants in the PTSD group also had more depressive symptoms before deployment (t=−4.463, df=65, p<0.001). However, the level of these symptoms was within the normal range. On average, participants in the PTSD group had experienced more potentially traumatic experiences during childhood (t=−3.161, df=66, p=0.002). The two groups did not differ in the number of smokers, alcohol use, and body mass index before deployment. Medication use was very limited (angiotensin converting enzyme inhibitors and diuretics [N=1]; glucocorticoids, local use [N=1]; and antibiotics [N=1]). Excluding these participants from the analyses did not change the results.
Pre- and Postdeployment Questionnaire Scores and Predeployment Matching Variables and Characteristics of Participants in the PTSD and Matched Comparison Groupsaa
| Add to My POL
|Characteristic||PTSD Group (N=34)||Comparison Group (N=34)||Analysis|
|Self-Rating Inventory for PTSD total score after deployment||44.18||5.59||25.15||4.05||<0.001|
|Symptom Checklist-90 depression subscale score after deployment||24.03||5.83||16.94||1.67||<0.001|
|Self-Rating Inventory for PTSD total score before deployment||28.35||4.18||24.65||2.99||<0.001|
|Sympton Checklist-90 depression subscale score before deployment||19.82||4.49||16.56||1.13||<0.001|
|Exposure to potentially traumatic deployment stressors||5.97||2.48||5.50||2.51||0.46|
|Age (years) during deployment||27.35||10.09||26.68||9.28||0.77|
|Number of previous deployments||0.70||1.16||0.65||0.85||0.87|
|Early Trauma Inventory score (number of experiences)||4.62||3.12||2.62||2.81||0.002|
|Body mass index (kg/m2) before deployment||25.45||3.66||24.59||2.53||0.61|
|Smoking before deployment (yes)||20||60.6||16||47.1||0.33|
|Period of deployment||1.000|
|Rank during deployment||0.72|
| Noncommissioned officer||5||14.7||7||20.6|
|Function during deployment||0.79|
| Both equally||6||17.6||4||11.8|
|Alcohol consumption per week before deployment||0.12|
| No alcohol||4||12.1||2||6.1|
| 1-20 units per week||26||78.8||31||93.9|
| >20 units per week||3||9.1||0||0.0|
Predeployment GR Number in PBMCs
We quantified the predeployment dexamethasone binding capacity of PBMCs, reflecting the GR number in these cells, by determining the specific binding of 100 nM 3H-dexamethasone in a whole cell binding assay. The data demonstrate a higher GR number in PBMCs before deployment in participants who developed a high level of PTSD symptoms after deployment (t=−4.343, df=66, p<0.001) (Figure 1). Moreover, logistic regression analysis revealed that the predeployment GR number significantly predicted inclusion in the PTSD group after deployment, controlling for potentially traumatic childhood experiences and predeployment questionnaire scores for PTSD and depressive symptoms (W=10.077, p=0.002). The logistic regression showed that the risk for inclusion in the PTSD group after deployment increased 7.5-fold with each GR increase of 1,000 (unadjusted odds ratio=7.53; 95% confidence interval=2.16-26.20). We examined whether the GR number before deployment was related to predeployment age, rank, number of previous deployments, smoking, alcohol use, body mass index, potentially traumatic childhood experiences, and questionnaire scores for PTSD and depressive symptoms (Table 3). None of these variables were associated with the GR number.
Predeployment Group Differences in the Number of Glucocorticoid Receptor (GR) Binding Sites in Peripheral Blood Mononuclear Cells (PBMCs) of Participants in the PTSD and Matched Comparison Groupsa
a Analyses were performed in duplicate in one predeployment sample per individual. Data are presented as means and standard deviations. Group differences were tested using t tests. Graph A indicates the number of GR binding sites in PBMCs in participants with PTSD symptoms 6 months after deployment (N=34; range: 519-3,816) and in matched comparison subjects (N=34; range: 471-2,939). Statistical significance was achieved (p<0.001). Graph B indicates the number of GR binding sites in PBMCs in participants in the PTSD group with (N=15; range: 519-3,816) and without (N=19; range: 946-3,255) comorbid depressive symptoms 6 months after deployment.
Predictive Value of Predeployment Characteristics for Predeployment Number of Glucocorticoid Receptor Binding Sites in Peripheral Blood Mononuclear Cells in the PTSD and Matched Comparison Groupsaa
| Add to My POL
|Number of previous deployments||−0.006||−0.031||0.98|
|Body mass index (kg/m2)||0.024||0.160||0.87|
|Early Trauma Inventory score (number of experiences)||0.025||0.166||0.87|
|Self-Rating Inventory for PTSD total score||0.006||0.035||0.97|
|Sympton Checklist-90 depression subscale score||0.179||1.124||0.27|
Comorbid Depressive Symptoms
Previously, we reported a higher preexisting GR number in PBMCs among deployed military personnel who had developed fatigue and depressive symptoms 6 months after deployment (19). Therefore, we investigated whether the observed group difference in the GR number in the present study was related to development of comorbid depressive symptoms within the PTSD group. However, the predeployment GR number did not differ between PTSD participants with (N=15) and without (N=19) comorbid high levels of depressive symptoms after deployment (Figure 1).
Five participants in the PTSD group and five comparison subjects were also included in our previous study (19). Therefore, we investigated whether the observed group difference in the GR number in PBMCs was influenced by the inclusion of these participants. The group difference in this GR number remained highly significant after removal of these participants (t=−3.215, df=56, p=0.002).
GR Subtypes and GILZ, SGK-1, and FKBP5 mRNA Expression
Using quantitative real-time polymerase chain reaction analysis, we assessed whether the observed group difference in the predeployment GR number in PBMCs was associated with group differences in predeployment GR mRNA expression. There were no significant group differences in the level of predeployment mRNA for GR-α, GR-β, and GR-P. Furthermore, there were no significant group differences in the level of predeployment GILZ and SGK-1 mRNA expression. In addition, there were no group differences in predeployment FKBP5 mRNA expression levels in PBMCs.
Cortisol and Composition of PBMC Population
The observed group difference in the GR number in PBMCs was not associated with group differences in plasma cortisol or the subset composition of PBMCs (Table 4). However, in the PTSD group, a higher percentage of CD19+ B cells as well as a lower percentage of CD8+ T cells nearly reached significance. The group difference in dexamethasone binding capacity of PBMCs remained significant after controlling for these variables (F=15.138, df=1, 64, p<0.001).
Predeployment Plasma Cortisol Levels, Absolute Leukocyte Numbers, and Relative Distribution of Peripheral Blood Mononuclear Cell Population of Participants in the PTSD and Matched Comparison Groupsaa
| Add to My POL
|Variable||PTSD Group (N=34)||Comparison Group (N=34)||p|
|CD19+ B-lymphocytes (%)||12.01||3.90||10.16||3.88||0.54|
|CD3+ T-lymphocytes (%)||68.70||8.08||69.60||8.22||0.65|
|CD8+ T-lymphocytes (%)||33.83||6.49||37.48||8.52||0.05|
|CD4+ T-lymphocytes (%)||40.74||8.23||38.57||8.66||0.29|
We investigated whether the observed predeployment group difference in the GR number in PBMCs was still present after deployment. The PTSD group had a higher GR number at 1 (t=−2.117, df=66, p=0.04) and 6 (t=−2.366, df=66, p=0.02) months after deployment relative to the comparison group.
This study reveals that high preexisting dexamethasone binding capacity of PBMCs, reflecting a high GR number in these cells, represents a vulnerability factor for development of PTSD symptoms after military deployment. Before deployment, a higher GR number in PBMCs was present in military personnel who developed a high level of PTSD symptoms 6 months after deployment. The predeployment GR number was a strong predictor of inclusion in the PTSD group after deployment; with each increase of 1,000 GRs, the odds ratio for the presence of PTSD symptoms after deployment increased 7.5-fold. The observed group difference was still present after deployment.
Although the PTSD group reported more predeployment symptoms of PTSD and depression, the predeployment GR number was not influenced by these predeployment symptoms. Furthermore, group differences in the PBMC composition also did not influence the results. In addition, group differences in the GR occupancy by endogenous glucocorticoids most likely did not influence the observed group differences, since plasma cortisol levels did not differ between groups. Moreover, cells were equilibrated in a culture medium and extensively washed before the addition of 3H-dexamethasone to prevent possible effects of binding of endogenous glucocorticoids.
We previously reported a higher preexisting GR number in PBMCs among military personnel with depressive symptoms and severe fatigue after deployment (19). In the present study, almost one-half of the participants in the PTSD group reported high levels of depressive symptoms after deployment, which is consistent with comorbidity estimates for PTSD and major depressive disorder in war veterans (25). However, there was no difference in the pre-deployment GR number between participants with and without depressive symptoms in the PTSD group. Therefore, the present finding is related to development of PTSD symptoms after deployment and not to development of depressive symptoms in a subgroup of participants.
Although changes in GR functioning in PTSD and major depressive disorder are supposed to be different (3, 26), a high preexisting GR number in PBMCs is related to both PTSD and depressive symptoms after deployment (19). Thus, a high preexisting GR number in PBMCs may be a predictor of mental health problems specifically induced by stress and/or trauma.
A limitation of the present study is that mental health status was not assessed with structured clinical interviews. Therefore, it is uncertain whether participants in the PTSD group fulfilled DSM-IV diagnostic criteria for PTSD. Individuals with a subclinical level of PTSD symptoms, therefore, might also have been included in the PTSD group. However, it is reasonable to assume that the questionnaires are a reliable reflection of the presence of PTSD symptoms (21).
Currently, various subtypes of the GR have been identified. GR-α is the most abundant variant and is transcriptionally active (27), while GR-β is less ubiquitously expressed, does not bind glucocorticoids, and has limited transcriptional capability (28). Additionally, GR-P is widely expressed and increases GR-α activity in cell lines (29, 30). To our knowledge, there are no published data on GR mRNA expression levels in PTSD. However, within major depressive disorder, the peripheral GR number (for review see reference 26) and peripheral GR-α mRNA expression (31) both appear to be decreased. Our results show that the observed group difference in the GR number in PBMCs was not accompanied by group differences in GR mRNA expression level. This result implies that the higher GR number in the PTSD group is mediated by posttranscriptional mechanisms. Additionally, we cannot completely exclude that a group difference in mineralocorticoid receptor level influenced the observed higher dexamethasone binding capacity of PBMCs in the PTSD group.
Alterations in FKBP5 expression influence glucocorticoid binding capacity of PBMCs (15). Segman et al. (32) reported that the amount of up-regulation of FKBP5 mRNA expression after trauma predicted the presence of PTSD 4 months later. In addition, Yehuda et al. (33) found decreased FKBP5 mRNA expression levels in PTSD patients. In our study, however, there was no group difference in FKBP5 mRNA expression before deployment. We also did not observe group differences in the GR target genes GILZ and SGK-1. It may well be possible, however, that compensation by other pathways that regulate these genes compensate for a putative difference in GR activity.
The observed higher GR binding in the PTSD group may be associated with the presence of specific GR polymorphisms, which could, for example, alter posttranscriptional mechanisms influencing GR protein expression or directly affect the binding capacity of the GR. Bachmann et al. (34) reported no association between PTSD and the presence of GR polymorphisms BC1I and N363S. Currently, however, more polymorphisms in the GR gene have been identified (for review see reference 35), which were not yet investigated in relation to development of PTSD at the time of the study. Whether genetic factors are involved in the higher GR number in PBMCs observed in the PTSD group will be investigated in a larger sample.
In our study, participants in the PTSD group experienced almost twice as many potentially traumatic events during childhood than comparison subjects. Although the number of traumatic experiences during childhood was not related to the GR number in PBMCs, it is possible that a higher GR number in PBMCs and exposure to childhood trauma together increase the vulnerability for PTSD development after deployment. In this context, Bet et al. (36) reported that the co-occurence of both the GR polymorphisms 22/23EK and 9β with childhood trauma is associated with increased risk for development of major depressive disorder.
Peripheral GRs are often considered to be an accessible model for GR expression in the brain. Nonetheless, whether the mechanisms of regulation of central and peripheral GRs are similar has not been studied extensively. In rats, neuronal and lymphoid cytosolic GRs are similar in glucocorticoid affinity and specificity (37), and cytosolic GRs in the brain and peripheral immune tissues are down-regulated after chronic corticosterone administration following adrenalectomy (38, 39). It remains unknown whether and how the observed higher predeployment GR number in PBMCs is involved in the pathophysiology of PTSD. We speculate that our finding of a higher GR number in PBMCs may be paralleled in the brain. If a suitable radiotracer could be found, positron emission tomography neuroreceptor mapping of GRs could provide valuable information about GR-mediated abnormalities in the brain (40).
The present study is the first, to our knowledge, to show that a high GR number in PBMCs before military deployment was a strong predictor of the presence of PTSD symptoms after deployment and was maintained until at least 6 months after deployment. We propose that a high GR number in PBMCs may be a biomarker of increased risk for development of PTSD symptoms after trauma exposure.
The authors thank Col. C. IJzerman and the commanders and troops for their time and effort. The authors also thank Kim Kroezen, Anne Muilwijk, Lt. Maurits Baatenburg de Jong, Jessie Smulders, and Sgt. Maj. Martijn Derks for organizing the data acquisition as well as Karima Amarouchi, Linda Schild, Zabi Mohklis, Marijke Tersteeg-Kamperman, Esther Rudolph, and Jitske Zijlstra for technical assistance in the framework of the PRISMO project.