Figure 1. Cortisol levels at baseline and after low-dose dexamethasone suppression. A, Baseline cortisol levels in individuals with and without posttraumatic stress disorder (PTSD) stratified by the rs9296158 A carrier status. In individuals carrying rs9296158 GG, PTSD symptom severity was negatively correlated with baseline cortisol level (Pearson correlation = −0.328), and in individuals carrying AA/AG, PTSD symptom severity was positively correlated with baseline cortisol level (Pearson correlation = 0.157). B, Levels of serum cortisol suppression after low-dose (0.5-mg) dexamethasone suppression in individuals with and without PTSD stratified by rs9296158 A carrier status. Posttraumatic stress disorder was associated with increased dexamethasone suppression only in risk allele A carriers. The triple interaction of rs9296158 A carrier status and PTSD status on cortisol suppression was significant (P = .03). To convert cortisol to nanomoles per liter, multiply by 27.588.
Figure 2. Validation and replication of FKBP5 expression differences in posttraumatic stress disorder (PTSD) using quantitative polymerase chain reaction (qPCR). The PTSD Symptom Scale (PSS) total scores are plotted on the x-axis against the expression levels as measured by log2 expression for the microarray data and 2-ΔCT (change in crossing thresholds) for the qPCR validation and replication on the y-axis. The Pearson correlation coefficients for the rs9296158 A carrier status groups are shown. Post hoc genotype-dependent differences of the correlation coefficients were established by converting them into Fisher z scores.10 The z score differences between the 2 genotype groups were significant for the microarray expression (z = 2.17; significant at P = .01), qPCR validation (z = 2.12; significant at P = .05), and qPCR replication (z = 2.87; significant at P = .01). RT-PCR indicates reverse transcriptase–polymerase chain reaction.
Mehta D, Gonik M, Klengel T, Rex-Haffner M, Menke A, Rubel J, Mercer KB, Pütz B, Bradley B, Holsboer F, Ressler KJ, Müller-Myhsok B, Binder EB. Using Polymorphisms in FKBP5 to Define Biologically Distinct Subtypes of Posttraumatic Stress DisorderEvidence From Endocrine and Gene Expression Studies. Arch Gen Psychiatry. 2011;68(9):901-910. doi:10.1001/archgenpsychiatry.2011.50
Author Affiliations: Max-Planck Institute of Psychiatry, Munich, Germany (Drs Mehta, Gonik, Klengel, Menke, Pütz, Holsboer, Müller-Myhsok, and Binder and Mss Rex-Haffner and Rubel); Department of Psychiatry and Behavioral Sciences, Emory University School of Medicine (Drs Bradley, Ressler, and Binder), VA Medical Center (Dr Bradley), and Yerkes National Primate Research Center (Dr Ressler), Atlanta, Georgia; and Howard Hughes Medical Institute, Chevy Chase, Maryland (Ms Mercer and Dr Ressler).
Context Polymorphisms in the gene encoding the glucocorticoid receptor (GR) regulating co-chaperone FKBP5 have been shown to alter GR sensitivity and are associated with an increased risk to develop posttraumatic stress disorder (PTSD).
Objective To investigate interactions of the FKBP5 single-nucleotide polymorphism rs9296158 and PTSD symptoms on baseline cortisol level, low-dose dexamethasone suppression, and whole-blood gene expression.
Design Association of FKBP5 genotypes and PTSD symptoms with endocrine measures and genome-wide expression profiles.
Setting Waiting rooms of general medical and gynecological clinics of an urban hospital at Emory University.
Participants The 211 participants were primarily African American (90.05%) and of low socioeconomic status and had high rates of trauma and PTSD.
Main Outcome Measures Baseline and post–dexamethasone suppression cortisol measures and gene expression levels.
Results In our endocrine study, we found that only risk allele A carriers of rs9296158 showed GR supersensitivity with PTSD; in contrast, baseline cortisol levels were decreased in PTSD only in patients with the GG genotype. Expression of 183 transcripts was significantly correlated with PTSD symptoms after multiple testing corrections. When adding FKBP5 genotype and its interaction with PTSD symptoms, expression levels of an additional 32 genes were significantly regulated by the interaction term. Within these 32 genes, previously reported PTSD candidates were identified, including FKBP5 and the IL18 and STAT pathways. Significant overrepresentation of steroid hormone transcription factor binding sites within these 32 transcripts was observed, highlighting the fact that the earlier-described genotype and PTSD-dependent differences in GR sensitivity could drive the observed gene expression pattern. Results were validated by reverse transcriptase–polymerase chain reaction and replicated in an independent sample (N = 98).
Conclusions These data suggest that the inheritance of GR sensitivity–moderating FKBP5 polymorphisms can determine specific types of hypothalamic-pituitary-adrenal axis dysfunction within PTSD, which are also reflected in gene-expression changes of a subset of GR-responsive genes. Thus, these findings indicate that functional variants in FKBP5 are associated with biologically distinct subtypes of PTSD.
Posttraumatic stress disorder (PTSD) is a severely debilitating psychiatric condition characterized by persistent symptoms of intrusive reexperiencing, avoidance, and hyperarousal following exposure to a traumatic event. Although a lifetime trauma incidence of 40% to 90% in the general population has been reported, the overall lifetime prevalence for PTSD ranges between 7% and 12%.1,2 This suggests that individuals with a specific genetic susceptibility might have a higher risk of developing PTSD after experiencing traumatic events. Twin and family studies have provided evidence for a genetic component in the development of PTSD.3- 5
While studies investigating main genetic effects associated with PTSD have yielded inconsistent findings,6- 8 there is growing evidence supporting the role of gene × environment interactions in PTSD.9 For example, polymorphisms in the serotonin transporter gene as well as in the glucocorticoid receptor–regulating co-chaperone FK506 binding protein 5 (FKBP5) have been found in more than 1 cohort to have significant interactions with environmental exposure on risk for PTSD.9- 11 Additionally, the syndrome of PTSD, as with many other psychiatric disorders, is clinically variable. The identification of biologically distinct subsets of patients with PTSD would provide an important tool for improving our mechanistic, diagnostic, and therapeutic approaches to PTSD. This study is predicated on the possibility that there are subtypes of PTSD that can be differentiated based on differential physiological stress reactivity that may, in part, be genetically determined.
The stress hormone system, or hypothalamus-pituitary-adrenal (HPA) axis, is likely an important mediator of such gene × environment interactions because it is activated by traumatic events and dysregulated in a number of psychiatric disorders, including PTSD.12 In fact, HPA axis dysfunction with higher glucocorticoid receptor (GR) sensitivity is described to be one of the endocrine hallmarks of PTSD.13 A number of studies have shown altered stress hormone/cortisol levels in PTSD and enhanced GR sensitivity as measured by a low-dose dexamethasone suppression test as well as ex vivo binding assays.14- 17
Overall, these studies reflect a sensitization of the HPA axis in response to exposure to stressors, with changes in GR function playing a critical role in the pathophysiology of PTSD. The function of the GR is dependent on a large molecular complex, including chaperones and co-chaperones as well as coactivators and corepressors.18,19 One of the co-chaperones strongly regulating GR function is FKBP5, which is part of the mature GR heterocomplex and is known to regulate GR sensitivity.20 When bound to the steroid receptor, FKBP5 decreases its affinity for the ligand and prevents translocation to the nucleus.20,21 By being a modulator of GR sensitivity, FKBP5 is an interesting candidate gene for PTSD. In fact, 2 studies so far have reported interactions of functional polymorphisms within FKBP5 with early trauma to predict adult PTSD.10,11 In addition to these interaction effects on PTSD symptoms, we have also reported FKBP5 genotype–dependent difference in GR sensitivity in PTSD.10 Only patients with PTSD carrying the risk alleles displayed increased GR sensitivity, suggesting that endocrine changes observed with PTSD might be restricted to a genetically defined subset of patients.
In our previous report, rs9296158 showed the strongest interaction with early trauma on PTSD symptoms10 and was chosen as a tagging single-nucleotide polymorphism (SNP) for a functional haplotype that includes rs1360780 and rs3800373. Within this haplotype, risk alleles for PTSD are associated with an increased induction of FKBP5 messenger RNA by glucocorticoids, increased FKBP5 protein levels, and, as a consequence, a decrease in GR sensitivity and negative feedback regulation of the HPA axis in healthy controls.10,22,23 Genotype-dependent endocrine and messenger RNA effects in PTSD are less well studied.
We hypothesized that functional polymorphisms within FKBP5 are associated with differential GR sensitivity, which together may define biologically distinct subsets of patients with PTSD. To test this, we investigated interactions of the FKBP5 SNP rs9296158 and PTSD symptoms on baseline cortisol level, low-dose dexamethasone suppression, and whole-blood gene expression in a highly traumatized, low-income, primarily African American cohort with high rates of PTSD.24,25
The participants were a subset of a larger study investigating the contribution of genetic and environmental factors in PTSD10,25,26 and they had each experienced at least 1 traumatic event. Individuals were recruited between 2007 and 2010 at Grady Memorial Hospital, Atlanta, Georgia.24 All study procedures were approved by the institutional review boards of Emory University School of Medicine and Grady Memorial Hospital and all subjects gave written informed consent to the study.
Of the 219 participating individuals, 211 had available serum samples for baseline cortisol level only, while 115 also had cortisol measures after low-dose dexamethasone suppression. All individuals had DNA samples from whole blood or saliva and whole-blood RNA. Of 219 individuals, 10 samples were excluded from RNA analysis because of a low RNA integrity number (<5.0).27 Of 209 individuals, 111 individuals served as discovery samples and 98 served as replication samples for gene expression (eFigure 1 and eTable 1).
The PTSD Symptom Scale (PSS) was used as a measure of PTSD symptoms.28 The PSS frequency items were summed to obtain a continuous measure ranging from 0 to 48. The mean (SD) PSS score across all 219 individuals was 13.12 (12.33) points, with values more than 20.0 usually considered to show clinically significant PTSD symptoms. Application of the DSM-IV criteria for PTSD to the PSS frequency items was used to create a proxy variable for PTSD diagnosis status.
The Clinician-Administered PTSD Scale (CAPS) score was available for a subset of 156 of the 219 individuals.29 The mean (SD) current CAPS score in the 156 individuals was 39.95 (40.38). There was a significant correlation between PSS score and current CAPS score (Pearson r = 0.657; P < .01).
The Beck Depression Inventory was administered to measure depression symptoms.30 In this sample, the Beck Depression Inventory score ranged from 0 to 55 (mean [SD], 14.26 [11.63]). Values more than 18.0 are considered to show clinically relevant depressive symptoms.
The Trauma Events Inventory was the primary measure of non–child abuse trauma in this study.31 Total numbers of different types of non–child abuse trauma among the participants were summed. The mean (SD) score for non–child abuse trauma reported within this study sample was 3.90 (2.61).
The Childhood Trauma Questionnaire was used as a continuous measure of child abuse.32 The Childhood Trauma Questionnaire score ranged from 25 to 111, with a mean (SD) of 42.20 (18.59).
Open-ended questions as well as the Kreek-McHugh-Schluger-Kellogg Scale were used to quantify self-exposure to opiates, cocaine, alcohol, and tobacco33 and to assess overall current and past substance use.
The subjects were characterized using a low-dose dexamethasone suppression test. After patients fasted, blood was collected between 8 and 9 AM for baseline serum cortisol measurements. The subjects received 0.5 mg of dexamethasone orally at 11 PM. Blood was collected between 8 and 9 AM the next day for serum cortisol measurement after dexamethasone stimulation. Serum cortisol concentration was measured with a commercial radioimmunoassay kit (Diagnostic Systems Laboratories, Webster, Texas).
DNA was extracted from whole blood or saliva collected in Oragene saliva kits (DNA Genotek, Ottawa, Ontario, Canada) using the Qiagen M48 automated extraction system (Qiagen, Valencia, California). We genotyped rs9296158 using a TaqMan-based assay, as described before.10
After patients fasted, whole blood was collected between 8 and 9 AM in Tempus RNA tubes (Applied Biosystems, Darmstadt, Germany) for baseline cortisol measurements. RNA was isolated using the Versagene kit (Gentra Systems, Minneapolis, Minnesota) and quantified using the NanoPhotometer (Implen, Westlake Village, California) and quality checks were performed on the Agilent Bioanalyzer (Agilent Technologies, Santa Clara, California). RNA was not globin reduced since our preliminary data showed that although globin reduction increased the sensitivity, it altered the expression of nonglobin transcripts in a nonreproducible manner (data not shown). Two hundred fifty nanograms of total RNA were reverse transcribed to complementary DNA, converted to complementary RNA, and biotin labeled using the Ambion kit (AMIL1791; Applied Biosystems). Seven hundred fifty nanograms of complementary RNA were hybridized to Illumina HT-12 version 3.0 arrays (Illumina, San Diego, California) and incubated overnight (16 hours) at 55°C. Arrays were washed, stained with Cy3-labeled streptavidin, dried, and scanned on the Illumina BeadScan confocal laser scanner.
For validation and replication of the results, complementary DNA was synthesized from 200 ng of total RNA using Superscript II Reverse Transcriptase (Invitrogen, Darmstadt, Germany). Quantitative polymerase chain reaction (qPCR) was performed using the Universal Probe Library for FKBP5, IL18R1, and TBP on the Roche LightCycler 480 (Roche Applied Science, Penzberg, Germany). Assays were designed using Probe Finder Software (Roche Applied Science) and run in duplicate according to the manufacturer's protocol, except for a total reaction volume of 10μL (eTable 2).
Raw microarray scan files were exported using the Illumina BeadStudio program and loaded into R for downstream analysis (http://www.R-project.org). The data were transformed and normalized using the variance-stabilizing normalization.34 Of 48 804 probes present on the Human HT-12 version 3.0 assays, probes were selected that fulfilled the criteria of the Illumina probe detection P value of<.01 in 5% of the individuals. A total of 15 583 probes passing the filter criteria were used for subsequent analysis. Association analysis was performed in R. General linear models were constructed using the quantitative PSS scores, rs9296158 A carrier genotypic model, and interactions between the 2 as predictors, controlling for age and sex. The significance of association was estimated by a 2-tailed P value corresponding to the t ratio based on a t reference distribution. To assess the contributions of model covariates, pairwise comparison between the full and reduced models was performed using an analysis of variance F test. The results were corrected for multiple testing using the permutation of regressor residuals test implemented in the R package glmperm (http://cran.r-project.org/web/packages/glmperm/index.html). For each transcript, 10 000 permutations were performed. The permutation of regressor residuals test permutes the residuals from the linear regression models using the shuffle-Z method for permutation.35 Repeated-measures analysis was performed in SPSS (SPSS Inc, Chicago, Illinois) using the analysis of variance design. The Huynh-Feldt method was applied to correct for potential violations of the sphericity assumption.36
The BiblioSphere data mining software from Genomatix was used for pathway analysis of transcripts.37 BiblioSphere is a data-mining solution for extracting and studying gene relationships from literature databases and genome-wide promoter analysis.38 The signal transduction association filter with the highest specificity was used where 2 genes were co-cited in a sentence containing a pathway-associated term, with at least 1 of the co-cited genes bearing a Genomatix signal transduction pathway annotation used for the pathway (http://www.genomatix.de). The JASPAR CORE database via the cREMaG interface (http://184.108.40.206/cremag/) was used to assess the distribution of steroid receptor binding sites.39
Raw data from the Roche LightCycler 480 System were extracted using LinRegPCR software.40 The crossing thresholds (CT) from the technical replicates were averaged across all samples. The transcripts were normalized to housekeeping gene TBP. NormalizedΔCT values were converted into the linear form by 2-ΔCT transformation.41 The transformedΔCT values were regressed against the PSS score, rs9296158, and their interaction, while adjusting for age and sex, in R. For calculation of correlations between expression and PSS scores, negative correlations of theΔCT were used to make the direction of effect comparable with the microarray data.
Phenotypic characteristics of the 211 individuals used for endocrine analysis and 111 individuals used for microarray analysis are shown in Table 1. For comparisons, individuals were separated into groups of PTSD and non-PTSD based on DSM-IV criteria for PTSD using the responses to the PSS interview.
No differences in age, sex, ethnicity, and substance abuse were observed between the groups (P > .05) in both samples. There were no significant differences in adult trauma exposure (P = .11) in the microarray data set, but there was a significant difference in the endocrine data set (P = .001). Differences in child abuse were significant (P < .001 for the endocrine data set; P = .004 for the microarray data set). The severity of the PTSD symptoms measured by the PSS and current CAPS was significantly higher among individuals in the PTSD group (P < .001 for PSS score and P < .001 for CAPS score in the endocrine data set; P < .001 for PSS score and P < .001 for CAPS score in the microarray data set). Individuals with PTSD had higher Beck Depression Inventory scores than those without PTSD (P < .001 for the endocrine data set and P < .001 for the microarray data set). To determine whether the comorbidity between PTSD and depression would lead to multicollinearity, diagnostic tests for variance inflation factor and tolerance were estimated using SPSS. The variance inflation factor was 1.1 and the tolerance was 0.903, indicating no evidence for multicollinearity.
Previously, we reported a possible genotype-dependent difference in the effects of PTSD on GR sensitivity (n = 80), with higher GR sensitivity only in patients with PTSD carrying the risk allele.10 We now reinvestigated the effects of both PTSD symptoms and rs9296158 genotype in a second larger independent sample (n = 211 for baseline cortisol level and n = 115 for the dexamethasone suppression test).
We first regressed baseline cortisol levels against the PSS score, FKBP5 SNP rs9296158 A carrier status, and their interaction term, adjusting for age and sex. Significant main effects of PSS scores (P = .003) and rs9296158 (P = .007) and significant interaction between the 2 (P = .002) on baseline cortisol levels were observed. Posttraumatic stress disorder symptom severity was negatively correlated with baseline cortisol level only in individuals carrying rs9296158 GG (Pearson correlation = −0.328 for GG and 0.157 for AA/AG, both significant at P = .01) (Figure 1A). The association remained significant after correcting for possible confounding effects of child abuse and adult trauma severity (PSS scores: P = .03; rs9296158: P = .04; and PSS score × rs9296158: P = .02). Overall, FKBP5 genotype was not associated with the level of child abuse (P = .55).
Using repeated-measures analysis of variance, we next investigated the interaction of rs9296158 A carrier status and PTSD status on cortisol suppression. We confirmed a significant triple interaction of the genotype and PTSD on the change in cortisol levels (P = .03; F114,1 = 4.963).
We observed a significant PTSD status × rs9296158 genotype carrier status interaction (P = .04), but no main genotype or disease effect, on the change in cortisol levels. As observed before, PTSD was associated with increased dexamethasone suppression only in risk allele A carriers. The triple interaction of the genotype and PTSD status on cortisol suppression remained significant after adjusting for child abuse and adult trauma (P = .04) (Figure 1B).
To analyze the effects of PTSD symptom severity on gene expression, the continuous PSS scores were regressed against the filtered 15 583 expression profiles. After age and sex adjustment and multiple testing correction using the permutation of regressor residuals, 183 transcripts were significant with P values<.01, none of which had been previously reported to be associated with PTSD. Of these, 85 transcripts remained significant after adjusting for possible confounding effects of child abuse and adult trauma (eTable 3 and eFigure 2).
To investigate whether the observed FKBP5 genotype-dependent difference in GR sensitivity in PTSD would influence gene expression patterns, we used linear regression analyses to test the contribution of the continuous PTSD symptom severity, rs9296158 A carrier status, and their interaction to expression levels. Forty-one transcripts were significant for the genotype × PTSD symptom interaction after correcting for multiple testing using the permutation of regressor residuals test and an interaction P value of<.05 (Table 2 and eFigure 3). Nine of these 41 transcripts had already been identified using only PTSD symptoms as a predictor as described earlier, so that 32 novel transcripts were identified using this model. These 32 genes included FKBP5 and IL18R1, where either the molecule itself or the ligand had previously been reported to be regulated by PTSD.42- 44
Analysis of gene expression data in the form of gene networks can provide a more comprehensive view of gene interactions underlying complex traits of interest.45 To investigate possible common components of regulatory pathways, the 41 transcripts were characterized using the BiblioSphere software.37 Pathway analysis of these 41 transcripts revealed a network that directly included 15 of the tested genes as well as associated transcription factors. These included CREB1, HDAC1, FOS, and STAT5B, which had previously been reported to be associated with PTSD symptoms42,43 (eFigure 4). In fact, STAT5B expression showed a significant correlation with FKBP5 expression in this data set, with a Pearson correlation of 0.25 and 2-tailed significance of P = .007.
We next investigated the distribution of steroid hormone biding sites, including GR response elements, androgen receptor response elements, and progesterone receptor response elements, within the 32 transcripts regulated by the SNP × PTSD symptom interaction vs the 183 transcripts associated with PTSD symptom severity alone. For this, the JASPAR CORE database was probed to test for a possible overrepresentation of these steroid receptor response elements. Comparisons between the 183 transcripts identified using the PTSD symptom scores and 32 transcripts significant for the rs9296158 × PSS score interaction revealed a highly significant overrepresentation of steroid response elements (GR response elements and androgen receptor response elements) among the 32 transcripts but not the 183 transcripts associated with PTSD symptom severity alone (eTable 4). No progesterone receptor response elements were predicted in any of the tested transcripts.
Because FKBP5 and IL18 had been reported in previous PTSD gene expression studies, we proceeded to validate FKBP5 and IL18R1 expression results in the 111 individuals from the microarray analysis using qPCR. For both FKBP5 and IL18R1, significant main effects for PSS score as well as significant interaction between PSS and rs9296158 were observed (FKBP5: PSS total score, P = .002552; PSS score × rs9296158, P = .006932; IL18R1: PSS total score, P = .012253; PSS score × rs9296158, P = .018549) in the qPCR data set.
Next, we proceeded to replicate our findings in an additional set of 98 individuals not used in the microarray analysis. We replicated the significant main effects for PSS score as well as the significant interaction between the 2 for FKBP5 and IL18R1 after adjusting for age and sex in the independent replication sample (FKBP5: PSS total score, P = .01967; PSS score × rs9296158, P = .00319; IL18R1: PSS total score, P = .03826; PSS score × rs9296158, P = .00408).
The genotype-dependent differences in correlation of FKBP5 messenger RNA levels and PSS score severity in the microarray data set, the qPCR validation, and the replication data set are depicted in Figure 2.
These data from endocrine measures as well as gene expression profiling suggest that the functional FKBP5 polymorphism rs9296158 is not only associated with an increased risk for PTSD following childhood trauma10,11 but might also define biologically distinct subtypes of this disorder. Endocrine measures and gene expression data were available in more than 200 traumatized individuals. In fact, we could replicate that only risk allele carriers of this SNP show the previously reported association of GR supersensitivity with PTSD.10,13 On the other hand, baseline cortisol levels were decreased in PTSD only in patients with the nonrisk GG genotype. This suggests that the inheritance of GR sensitivity–moderating FKBP5 polymorphisms10,22 can determine specific types of HPA-axis dysfunction within PTSD that were also reflected in gene-expression changes of a subset of GR-responsive genes.
A number of studies have reported both decreased baseline cortisol levels as well as increased GR sensitivity46 in patients with PTSD, although the findings are not uniformly consistent across studies.14,17,47 While we also observed these 2 endocrine features with PTSD, these findings are restricted to specific FKBP5 SNP genotype groups, which might be an explanation for some of the inconsistent findings in the literature. In this article, we replicated a relative GR resistance in nonaffected individuals carrying the risk allele, which is in line with our previous report.10 These data also support the previously reported decrease in negative feedback capacity in these individuals.23 Also, as previously reported, the genotype-dependent association with GR function seems to be state dependent, with increased GR sensitivity in risk allele carriers with PTSD. We are currently investigating the molecular basis of this state dependence and speculate that in addition to possible expression changes of other GR chaperones and other systems influencing the HPA axis, this might be mediated by epigenetic changes. Interestingly, the protective genotype, not the risk, was associated with decreasing baseline cortisol levels with increasing PTSD symptoms. In our sample, decreased baseline cortisol level and increasing GR sensitivity therefore do not seem to be strongly correlated, suggesting that other mechanisms besides efficiency of negative feedback of the HPA axis, possibly including stress reactivity, regulate baseline cortisol levels in PTSD. The possibility that these 2 measures are somewhat disjunct is also supported by the fact that while the finding of increased GR sensitivity is relatively consistent across studies, it is less so for baseline cortisol changes, with decreases, no changes, and increases reported in PTSD.16,48 A previous study had also reported that PTSD-related decreases in baseline cortisol level might be dependent on a functional GR polymorphism, a possible additional confounder in this association.49
The genotype-dependent endocrine differences with PTSD symptoms were also paralleled by genotype-dependent differences in whole-blood gene expression pattern, where expression levels of 183 transcripts were significantly associated with PTSD symptoms alone. Within these, none had been previously identified to be associated with PTSD. Adding the FKBP5 polymorphism to the model identified 32 transcripts not detected in the first model that showed a significant interaction between FKBP5 genotype and PTSD symptom severity. In fact, these transcripts showed strong genotype-dependent differences for the correlation of gene expression levels with PTSD symptom severity (eFigure 3), which likely has precluded their identification when analyzing the association with PTSD symptom severity alone. Among the top transcripts exhibiting the most significant gene expression interaction effects was FKBP5 itself. This is not surprising given the previously discussed influence of PTSD and FKBP5 on GR function and the fact that FKBP5 appears to be one of the most sensitive readouts of GR stimulation. FKBP5 has been shown to be one of the genes with the strongest induction by GR activation in a number of studies, including our own.50,51 Most notably, Yehuda and colleagues43 have also described a downregulation of FKBP5 gene expression in whole blood with PTSD, as we observed for rs9296158 GG carriers. Within these 32 genes, the IL18R1 transcript is also an interesting candidate. Zieker et al44 reported interleukin 18 to be differentially expressed in the peripheral blood of patients with PTSD, concordant with the differential regulation of its receptor IL18R1 in our cohort.
We also found that 15 of the genes regulated by the interaction between FKBP5 genotype and PTSD can be integrated into the same pathway (eFigure 4) using the BiblioSphere software.37 When expanding the annotation to include related transcription factors, we identified another 4 previously reported PTSD candidates: CREB1, FOS, HDAC1, and STAT5B42,43 to be part of this common regulatory pathway. We could confirm the significant correlation between expression levels of FKBP5 and STAT5B in this study, both of which had been previously shown to be regulated with PTSD by Yehuda et al.43 Our data suggest that a data mining approach could be informative in revealing coexpressed gene networks and might help identify transcripts that otherwise would be missed when not passing statistical thresholds of significance because of the subtle gene expression differences.
While we did identify a number of genes regulated by PTSD symptom severity itself, our data suggest that a network of genes is differentially regulated depending on FKBP5 genotype and related differences in GR sensitivity. In fact, there was a significant overrepresentation of GR response elements and androgen receptor response elements among these 32 transcripts significant for the interaction term compared with the 183 genes regulated by PTSD alone. Altogether, these findings suggest that the combined model of FKBP5 SNP genotype and PTSD symptoms seems most relevant for GR-dependent transcripts, while other predictive models might be necessary to capture differences in the expression of other gene classes.
While we were able to replicate some of the genes and pathways identified in previous studies, others could not be replicated. This is likely because of a number of reasons, including false-positive and false-negative associations in either study, differences in clinical study design, sample characteristics, type of microarrays, and statistical analysis.
A strength of our microarray gene expression study (n = 111 in total) is the large sample size, with 41 individuals with a probable PTSD diagnosis, which, to our knowledge, is larger than any of the previously reported expression array studies for this disorder. In the overall sample, we investigated 75 individuals with a probable PTSD diagnosis. We acknowledge that our main PTSD definition is symptom driven; using the PSS might have included individuals with other anxiety disorders also showing high PSS scores.52 However, about 70% of our cohort had CAPS-based current PTSD diagnoses and these concorded well with the PSS scores (Pearson r = 0.657, significant at a P = .01). The use of the PSS allowed us to include all individuals and to increase the power of the study by using a quantitative predictor vs using diagnostic cutoffs, a strategy now advocated by a number of investigators for biomarker discovery.42,53 A limitation of this highly traumatized sample is the heterogeneity regarding the timing of the trauma as well as the types of traumatic events. Despite this, our gene expression results appear to be comparable with other studies with much more defined index trauma, as described earlier.16,42
In this sample, all subjects had experienced at least 1 significant trauma; nonetheless, we noted significant differences in overall trauma exposure as well as the level of child abuse between individuals with high vs low PTSD symptoms. When correcting the analysis for adult trauma and child abuse exposure, about 50% of the associations with PTSD severity alone remained significant at a level of P = .05, indicating the importance of increased trauma exposure on gene expression, even when using traumatized individuals as controls (eTable 3). Differences in severity and timing of trauma exposure might also account for the slight differences in the PTSD × SNP interaction on FKBP5 messenger RNA expression in the discovery vs the replication sample (Figure 2 and eTable 1).
Independent replication is extremely important to assess the validity of gene expression findings, because more often than not specific changes reported for psychiatric disorders could not be replicated.54,55 In this study, we chose FKBP5 and IL18R1 for validation and replication, because they had previously been implicated in this disorder. We could, in fact, validate and replicate the association with genotype-dependent gene expression differences of both 2 transcripts using qPCR.
In summary, using a large sample of 219 traumatized individuals, we examined the effects of FKBP5 rs9296158 genotype, previously reported to increase the risk for PTSD with early trauma, on putative biomarkers for this disorder. In this study, we found that rs9296158 may be used to identify biologically different subtypes of PTSD in that the genotype groups differed with respect to PTSD-related changes in GR sensitivity. This was reflected in genotype- and PTSD-dependent differences in the expression of GR-dependent transcripts in whole blood. Our findings support the importance of using genetic markers and biological measures to increase the predictive validity of biomarkers in psychiatric disorders.56
Correspondence: Elisabeth B. Binder, MD, PhD, Max Planck Institute of Psychiatry, Kraepelinstrasse 2, Munich 80804, Germany (firstname.lastname@example.org).
Submitted for Publication: October 6, 2010; final revision received January 14, 2011; accepted March 4, 2011.
Published Online: May 2, 2011. doi:10.1001/archgenpsychiatry.2011.50
Financial Disclosure: Part of this work was funded by NeuroNova. Drs Binder, Menke, and Holsboer have filed a patent application as inventors of“Means and Methods for Diagnosing Predisposition for Treatment Emergent Suicidal Ideation (TESI)” (European application 08016477.5; international application PCT/EP2009/061575), and Drs Binder, Müller-Myhsok, and Holsboer have filed a patent application as inventors of“FKBP5: A Novel Target for Antidepressant Therapy” (international publication WO 2005/054500). Dr Holsboer is a founder and shareholder of Affectis Pharmaceuticals. Within the last 3 years, Dr Ressler has received research funding support from Lundbeck, and he has an unrelated role as cofounder of Extinction Pharmaceuticals for development of N -methyl-D-aspartate–based therapeutics. Dr Binder has current grant support from PharmaNeuroBoost.
Funding/Support: This work was supported by the Max Planck Society and National Institute of Mental Health grant MH071537. Dr Bradley receives grant support or has received awards from the American Foundation for Suicide Prevention and the American Psychoanalytic Association Psychoanalytic Research Fund. Within the last 3 years, Dr Ressler has received research funding support from the National Institute of Mental Health, the National Institute on Drug Abuse, the Burroughs Wellcome Foundation, and the National Alliance for Research on Schizophrenia and Depression. Dr Binder has current grant support from the National Institute of Mental Health and the Doris Duke Charitable Foundation.