Association of FKBP5 Polymorphisms and Childhood Abuse With Risk of Posttraumatic Stress Disorder Symptoms in Adults FREE

Posttraumatic stress disorder (PTSD) is a debilitating stress-related psychiatric disorder, with prevalence rates of at least 7% to 8% in the US population, and with much higher rates among combat veterans and those living in high-violence areas.^1- 3 Initially viewed as a potentially normative response to traumatic exposure,⁴ it became clear that not everyone experiencing trauma develops PTSD. Thus, a central question in research on PTSD is why some individuals are more likely than others to develop the disorder in the face of similar levels of trauma exposure.^5- 8 Although PTSD is the single disorder within the Diagnostic and Statistical Manual of Mental Disorders (Fourth Edition) (DSM-IV)⁹ that requires a specific environmental insult within its diagnostic criteria, it is becoming increasingly clear that there are critical roles for predisposing genetic and environmental influences in differentially mediating psychological risk to the traumatized individual.^10- 13

Child abuse occurs at disturbingly high rates and is a major public health problem.¹⁴ Despite the resilience of many abused children, child abuse significantly increases risk for impaired physical and psychological health and decreases adaptive functioning in adulthood.^15- 17 A matter of central importance to public health is the identification of factors related to risk and resilience in the wake of child abuse. Of particular relevance to this study is the well-established relationship of child abuse with adult PTSD.^18- 22 A number of factors may account for this relationship. PTSD in adults may represent a prolonged symptomatic reaction to traumatic child abuse.

The experience of child abuse is associated with an increased number of traumatic experiences across the lifespan. Child abuse may also increase vulnerability for the later development of PTSD by altering psychological (eg, attachment) and biological (eg, hypothalamic-pituitary-adrenal [HPA] axis function) developmental processes, including interaction with genetic factors. Although other non–child abuse types of traumatic experiences in childhood²³ (eg, a house fire or being in a vehicle crash) might be expected to negatively affect development, the most robust research to date points to child abuse and related family/interpersonal stressful life events in predicting a wide range of later psychological and physical health problems.

The reasons for this are not fully understood but some possible explanations are (1) compared with other types of traumatic events, child abuse is more likely to occur in the family context²⁴; (2) any 1 type of child abuse is associated with an increased likelihood of exposure to other types of abuse and to increased levels of family-related stressful events/parental dysfunction (eg, parental substance abuse)²⁵; and (3) compared with some other types of trauma exposure, child abuse may be more likely a repeated experience rather than a single event (eg, multiple incidents of sexual abuse by the same perpetrator over a number of years).²⁵ Our study focuses on the interplay between child abuse and polymorphisms in the FKBP5 gene, which is involved in the glucocorticoid-mediated stress response, in the prediction of adult PTSD.

The literature to date on the genetics of PTSD has recently been reviewed,^12- 13 with the resulting suggestion that gene × environment studies are needed to focus more on distinct endophenotypes and influences from environmental factors. A number of studies suggest that genetic factors contribute to the development of PTSD, with heritability estimates ranging from 30% to 40%.^11,26- 29 Candidate gene studies, however, have been inconclusive so far, usually limited by extremely low power to detect any but the strongest possible genetic effects (current published studies include sample sizes of <100).

A recent review¹² covering candidate genes in the serotonin, dopamine, glucocorticoid receptor (GR), γ-aminobutyric acid, apolipoprotein E, brain-derived neurotrophic factor, and the neuropeptide Y system finds that (1) the support for a relationship between the serotonin transporter gene and PTSD exists only in research on the interaction of this system with stressful life events in predicting depressive symptoms; (2) results on the dopamine system are inconsistent; (3) there is a lack of evidence for relationships with brain-derived neurotrophic factor, neuropeptide Y, or GR polymorphisms and PTSD with the exception of a finding between GR genotype and basal cortisol levels in a subgroup of patients with PTSD³⁰; and finally, (4) limited evidence for the γ-aminobutyric acid system and the apolipoprotein E system.

One of the first studies finding an interaction between a genetic polymorphism and child abuse in predicting psychopathology was the study by Caspi et al,³¹ which found that maltreated children with a monoamine oxidase A (MAOA) genotype conferring low levels of MAOA expression were more likely to develop conduct disorder and antisocial-personality disorder and to commit violent crimes as adults compared with those children with the high-activity MAOA genotype. A recent study³² has replicated this result, and a second study³³ replicated the effect in white participants but not in other participants (blacks, Hispanics, American Indians, Pacific Islanders, and others) in the sample. The largest group of genotype × environment studies has examined the interaction between variation at the 5HTTLPR (a complex-repeat polymorphism in the 5′ upstream region of SLC6A4, which encodes the serotonin transporter the 5HTTLPR) and stressful life events, including child abuse, in predicting depression.^34- 40 Very recently, Kilpatrick et al⁴¹ reported an interaction of this polymorphism with severity of trauma and level of social support with the development of PTSD following hurricane exposure as outcome, supporting the relevance of gene × environment interactions for this disease.

From a developmental perspective, HPA axis genes are strong candidates with respect to altering susceptibility to PTSD. Exposure to trauma and stress increase HPA axis activity, and PTSD has been associated with long-lasting alterations in HPA axis reactivity^42- 43 and specifically higher GR sensitivity.^44- 45 Polymorphisms in genes regulating GR activity may impact the acute effects of trauma on the HPA axis and thereby possibly impact long-term HPA axis regulation affecting the development of PTSD. A number of studies suggest that child abuse and neglect affect HPA axis functioning (reviewed by Watts-English et al⁴⁶). Several studies suggest that the depression-related HPA axis hyperactivity may be related to early life stress. For example, plasma corticotropin and cortisol, as well as cerebrospinal fluid corticotropin-releasing hormone (CRH) concentrations, correlate with perceived early life stress more than with current depression.^47- 48 Preclinical studies indicate that the persistent hyperactivity of the HPA axis associated with early life stress is mediated by a hyperactive CRH receptor 1 (CRHR1) system, with chronic overactivity of CRHR1 in limbic brain regions.^49- 50

FKBP5 is a co-chaperone of hsp90. It directly interacts with hsp90, which binds to the GR. FKBP5 also has been shown to regulate GR sensitivity. FKBP5 is part of the mature GR heterocomplex.⁵¹ On hormone binding, FKBP5 is replaced by FKBP4, which then recruits dynein into the complex, allowing translocation into the nucleus where the complex regulates expression of glucocorticoid-responsive genes by functioning as a transcription factor.⁵² FKBP5 expression is induced by glucocorticoids as part of an intracellular ultrashort negative feedback loop for GR activity.⁵³

Overexpression of human FKBP5 in vitro reduces hormone binding affinity⁵⁴ and nuclear translocation of GR.⁵⁵ Naturally occurring overexpression of FKBP5 causes GR resistance in New World monkeys,^54,56 which is accompanied by increased plasma cortisol levels. Furthermore, potentially functional single-nucleotide polymorphisms (SNPs), or SNPs in very strong linkage disequilibrium with a functional variant, appear to alter FKBP5 function. The rare homozygous genotypes of FKBP5 SNPs (rs4713916, rs1360780, and rs3800373) were associated with higher FKBP5 expression in human blood monocytes as well as with a stronger induction of FKBP5 messenger RNA (mRNA) by cortisol. This was accompanied by less corticotropin release measured in patients who were depressed with the combined dexamethasone-CRH test.⁵⁷

The same alleles of rs3800373 and rs1360780 were associated with increased peritraumatic dissociation in children after medical trauma.⁵⁸ Higher levels of peritraumatic dissociation have been shown to be predictors of PTSD in adults.¹⁰ In addition, the extent of up-regulation of FKBP5 mRNA in peripheral blood mononuclear cells only hours after a trauma has been shown to correlate with the development of PTSD at 4 months.⁵⁹ These data suggest that FKBP5 could be an important candidate gene in trauma-related HPA axis disturbances. We therefore hypothesized that the putative functional SNPs in FKBP5 moderate the development of PTSD. Because early trauma, PTSD, and FKBP5 SNPs have all been shown to influence GR resistance, we also hypothesize that variants in this gene may alter the impact of early trauma or PTSD on GR sensitivity and address this by investigating the dexamethasone suppression test (DST) in a subsample of individuals.^44,57,60

Our study addresses the role of polymorphisms in FKBP5 in predicting PTSD, as well as the PTSD symptom–associated changes in GR sensitivity, in a highly traumatized, inner city sample. Specifically, we address whether these polymorphisms interact with increasing levels of child abuse and increasing levels of non–child abuse trauma exposure to predict PTSD symptomatology during adulthood.

Data were collected to investigate the roles of genetic and environmental factors in predicting the development of PTSD in a population of urban, low-income, predominantly black men and women. To determine the race/ethnicity composition of the sample as part of the screening procedures, we asked participants to self-identify their race/ethnicity. Their response was coded into 5 common categories (black, white, Hispanic or Latino, Asian, mixed) or other (participants checked “other” when they thought their race/ethnicity was not included in the other 5 categories) (Table 1). The mean (SD) age in the sample was 40.8 (13.8) years, ranging from 18 to 81 years.

Screen interviews were completed with 900 participants approached while in the waiting rooms of primary care or obstetrical-gynecological clinics of Grady Memorial Hospital in Atlanta, Georgia, between 2005 and 2007. Approximately 58% of those approached to participate in the study agreed to do so. Participants completed a battery of self-report measures that took 45 to 75 minutes to complete (dependent in large part on the extent of the participant's trauma history and symptoms). All measures were obtained by verbal interview. Each person was paid US $15.00 for participation in this phase of the study. Each participant also provided a saliva sample for DNA extraction (participants were paid the same amount for completion of screening data even if they declined to provide a DNA sample). Written and verbal informed consent was obtained for all participants, and all procedures in this study were approved by the institutional review boards of Emory University School of Medicine and Grady Memorial Hospital, Atlanta, Georgia.

Modified PTSD Symptom Scale . The modified PTSD Symptom Scale (PSS) is a psychometrically valid 17-item self-report scale assessing PTSD symptomatology^61- 65 over the prior 2 weeks. Consistent with prior literature, we summed the PSS frequency items (0 indicates not at all to 3 indicates ≥5 times a week) to obtain a continuous measure of PTSD symptom severity ranging from 0 to 51. For this sample, the PSS frequency items had standardized α = .90 (mean [SD], 13.81 [11.96]). No clearly established PSS cutoff score for PTSD diagnosis has been established; however, DSM-IV criteria for PTSD can be applied to PSS frequency items to create a proxy variable for PTSD diagnostic status.

Clinician Administered PTSD Scale . The Clinician Administered PTSD Scale (CAPS) was also administered to a subset of 240 participants within 2 to 6 weeks after completing the screening assessment. We found a significant difference (F_1,239 = 56.55, P < .001) between average PSS score (mean [SD], 18.20 [12.82]) for participants positive for current PTSD based on the CAPS (applying DSM-IV decision rules with a symptom considered as present with a CAPS frequency score of ≥1 and intensity score of ≥2) compared with those participants not meeting current CAPS PTSD criteria (mean [SD], 7.51 [9.23]). In addition, 70% of those participants identified as PTSD positive by using the PSS-based proxy variable were also positive for current PTSD at the time of CAPS administration (the criterion A traumatic experiences used for CAPS diagnosis was not necessarily the same one used for obtaining PSS data). Study participants were asked to respond to the PSS items based on the trauma exposure (inclusive of child physical and sexual abuse and other life trauma exposure) that they believed had impacted them the most. Because we recorded data on the type of trauma identified for PSS score but not the age at which it occurred and because many study participants reported several incidents of the same type of trauma across the lifespan, we were unable to determine if the PSS data was based on a traumatic event occurring in childhood vs adulthood.

Beck Depression Inventory. The 21-item Beck Depression Inventory (BDI)⁶⁶ is a psychometrically validated, commonly used measure of depressive symptoms.⁶⁷ In our sample, the BDI had a standardized α coefficient of .99 and a mean (SD) of 14.4 (13.2).

Traumatic Events Inventory. The traumatic events inventory (TEI)^64- 65 assesses lifetime history of trauma exposure and is our primary measure of both child abuse and non–child abuse trauma. The TEI assesses past experience and frequency of 13 separate types of traumatic events as well as feelings of terror, horror, and helplessness with such events.

For the measure of child abuse, 2 of the TEI questions assessed physical abuse and sexual abuse occurring before age 14 years. Based on these questions, 17.6% of the sample reported a history of childhood physical abuse and 18.8% reported a history of childhood sexual abuse. With these data, we created a 3-level categorical variable reflecting number of types of child abuse: no child abuse (70.5% of sample), 1 type of either physical or sexual abuse (22.7%), or 2 types of both physical and sexual abuse (6.8%).

For the measure of non–child abuse trauma, the remaining TEI questions addressed other types of trauma exposure (Table 2). To summarize level of exposure to trauma other than child abuse, we summed total number of different types of non–child abuse trauma exposure reported by each participant. The mean number of types of non–child abuse trauma reported in our sample was 2.46 types (SD, 1.92). Using the same data, we created a 4-level categorical variable reflecting number of types of non–child abuse trauma experienced reported with 19.9% of participants in the no non–child abuse trauma exposure category, 22.3% of participants reporting 1 type, 31.7% of participants reporting 2 to 3 types, and 26.7% of participants reporting 4 or more types of non–child abuse trauma. The total number or types of trauma exposure variable was chosen because in our prior research^64- 65 and in other research on impact of trauma exposure⁶⁸ it relates in a predictable and consistent manner with a number of measures of adaptive functioning and trauma-related functioning.

Although the non–child abuse variable assesses traumatic experience other than child abuse across the lifespan, data from our sample suggests that these non–child abuse traumatic experiences primarily occur in adulthood. Specifically, we have data from a subset of study participants (n = 322) asking them the earliest age at which all of the traumatic events assessed by the TEI occurred. Averaged across all of the non–child abuse trauma categories, the mean (SD) for youngest age of trauma exposure in the sample was 25.1 (10.6) years and the median youngest of trauma exposure was 23 years.

Childhood Trauma Questionnaire. The Childhood Trauma Questionnaire (CTQ)⁶⁹ is a 28-item psychometrically validated, self-report inventory assessing self-reported level of child abuse and neglect.⁷⁰ We used the total score from this scale in our secondary analyses as a separate continuous measure of child abuse exposure.

Dexamethasone Suppression Test . A total of 80 participants were characterized using a low-dose DST.⁶⁰ Blood samples in the fasting state were obtained for baseline serum cortisol concentration and were drawn between 8 AM and 9AM on day 1. Participants then received 0.5-mg dexamethasone by mouth at 11 PM before the second day blood draw, and blood samples for serum cortisol concentration was collected again in the fasting state between 8 AM and 9AM for the day 2 collection. Using the cortisol concentration values from day 1 and day 2, the percentage suppression was calculated by 100 × [(cortisol day 1−cortisol day 2)/cortisol day 1]. Serum cortisol concentration was measured by using commercially available radioimmunoassay kits (Diagnostic Systems Laboratories, Webster, Texas) from replicate samples and with interassay quality control measures. Significance of the interaction effect of FKBP5 SNPs and PTSD on percentage suppression was determined by using permutation-based methods with age and sex as covariates (10 000 permutations) to avoid inflated P values due to outlier or small cell sizes. In addition, we used repeated measures analysis of variance to examine the interaction between probable PTSD diagnosis and FKBP5 genotypes on the change in serum cortisol suppression from day 1 to day 2. For cortisol analyses, a probable PTSD diagnosis was determined from the DSM-IV criteria applied to the modified PSS interview, as described previously.^64- 65

The individuals with available DST present only a subset of the whole sample. However, these 80 individuals are not different compared with the rest of the sample in PSS total score for samples with and without DST (mean [SD], 11.30 [11.5] vs 10.15 [11.4]; P = .40), BDI total score (13.55 [12.1] vs 14.57 [12.1], P = .45), race (93.7% vs 96.0% black, P = .92), degree of childhood trauma (28.8% and 10.0% in the 1 and ≥2 types of severe child abuse vs 23.7% and 6.7%, respectively; P = .28), and also non–child abuse trauma and sex. We also found no differences in the monthly household income and relationship status. However, the endocrinologically characterized individuals are older (mean age [SD], 43.3 [11.1] vs 39.3 [11.1]; P = .01). Due to the high comparability of all other variables and because all our analyses are corrected for age and sex, this subsample is representative of the sample as a whole.

DNA Extraction. DNA was extracted from saliva collected into Scope mouthwash (n = 46) or into Oragene saliva kits (DNAGenotek, Ottawa, Ontario, Canada) by using the Qiagen M48 automated extraction system. DNA was available for 762 individuals. One hundred thirty-eight individuals with no genotype information consisted of 11 in which DNA was not collected, 88 in which we attempted to collect DNA using Oragene-spit samples but DNA extraction failed completely in 2 separate extraction trials (these failures were either due to noncompliance of the study participants or failure to correctly break the seal that releases the stabilizing solution in the Oragene saliva kits), and 39 for which DNA had been collected but not extracted at the time of analysis. When comparing the 762 individuals with genotypes to those 138 participants without, we found no significant difference in age, sex, self-described race, income, employment status, relationship status, child abuse, and non–child abuse trauma exposure as well as PTSD and depression severity.

SNP Genotyping . Eight SNPs within FKBP5 were selected from dbSNP (http://www.ncbi.nlm.nih.gov/projects/SNP/) to include the 3 potentially functional SNPs (rs4713916, rs1360780, and rs3800373)⁵⁷; the other 5 SNPs were selected to span 120 kb covering the remainder of the FKBP5locus (NM_004117). We also genotyped 3 SNPs within the CRHR1 gene (rs110402, rs242924, and rs7209436), which have been recently reported to show strong interaction with childhood abuse on adult depression symptoms.⁷¹ All SNPs were genotyped using a TaqMan allelic discrimination assay⁷² on an ABI 7900HT instrument (Applied Biosystems, Foster City, California), using predesigned and validated TaqMan assay reagent kits containing 1 pair of polymerase chain reaction primers and 1 pair of fluorescently labeled probes (Applied Biosystems). Polymerase chain reactions were performed in 5-μL reaction volumes in a 384-well plate and contained 5 ng of DNA. The standard protocol that was provided with the kit was followed. Thermal cycler conditions were 95°C for 10 minutes, 40 cycles of 92°C for 15 seconds, and 60°C for 1 minute. The SDS version 2.2 software (Applied Biosystems) was used for allelic discrimination. For quality control, 9% of the samples were genotyped as duplicates across and within a 384-well plate. Only 0.03% of discordances were recorded and excluded from the analyses. Call rates for SNPs ranged from 93% to 97%, which included some samples that failed all assays and for which DNA quality was inferior. We used Haploview⁷³ to determine the linkage-disequilibrium structure of the SNPs within the FKBP5 gene and to test for Hardy-Weinberg Equilibrium. The SNP identifications, their location, the Hardy-Weinberg Equilibrium test P values, minor allele frequencies, and call rate percentages are shown in Table 3.

Ancestry Informative Markers. Population stratification is a potential confounder in association studies of complex traits. If such allele-frequency heterogeneity within a population is coupled to similar population heterogeneity in the outcome of interest, naive tests of interaction could yield spurious results if the resulting confounding is neglected in analysis. Because 95.2% of our sample participants were self-identified as being of black race (which is known to contain varying degrees of genetic admixture of sub-Saharan African, European, and Native American ancestry), we took steps to ensure our results were robust to potential confounding from population stratification by genotyping 134 ancestry informative markers^74- 75 in a subset of 280 individuals from our sample using single-base primer extension and a Beckman SNPstream instrument (Beckman-Coulter, Fullerton, California). Genotyping of ancestry informative markers and data analysis for estimates of ancestry were performed by DNA Print Genomics, Sarasota, Florida (http://www.dnaprint.com). Maximum likelihood estimates of individual biogeographical ancestry admixture were determined by using methods described previously⁷⁶ and by using a 4-continental population model. The choice of a 4-continential population model was based on previously published hypothesis-free cluster analyses of worldwide populations.^77- 78 This model lent itself to use the terms European (genetic ancestry shared among Europeans, Middle Eastern, and to a lesser extent South Asians), sub-Saharan African, East Asian, and Indigenous American (genetic ancestry shared among American Indians, Latin and South American Indians, and certain Central Asian populations). The names of these parental populations were chosen to describe extant elements of genetic structure and are arbitrary in that they are reflective of modern population distributions—not necessarily the distributions of the original parental populations 15 000 to 50 000 years ago. Individuals with more than 35 failed ancestry informative marker genotypes were excluded from the analysis (7%), because their individual biogeographical ancestry admixture could not be estimated with sufficient certainty.

Primary, Descriptive, and Secondary Analyses . Our primary analyses were the main FKBP5 SNP effects on PTSD symptom severity and their interaction effects with child abuse and non–child abuse trauma (690 participants for interaction with non–child abuse trauma and 678 participants for interaction with child abuse). Secondary analyses included BDI scores as outcome (number same as for PSS outcome), FKBP5 SNP genotype dependent correlation of child abuse severity and PTSD symptom severity (n = 728), and effects on the DST (n = 80), as well as analysis using CRHR1 genotypes (number same as for FKBP5 SNP interactions). Descriptive analyses were run on a sample of 900 individuals with valid PTSD measures.

Regression Analyses . We used a variation of linear regression to examine the effects of trauma exposure, FKBP5 genotypes, and their potential interaction on PSS scores. For FKBP5 genotype, we modeled a participant's genotype at each SNP under an additive model that is equivalent to coding of the number of copies of a reference allele that the participant possesses. For trauma exposure, we included both child abuse (the 3-level child abuse TEI score) and non–child abuse trauma exposure (using the 4-level non–child abuse trauma TEI score) as predictors. We chose to separate these 2 types of trauma exposure based on the above described research, suggesting that exposure to child abuse increases risk for development of PTSD both as a response to the child abuse itself and in response to other non–child abuse stressors. Because child abuse is related to increased likelihood of subsequent exposure to additional traumatic stressors, we wanted to control for the possibility that effects associated with our child abuse variables were not artifacts of increased overall trauma exposure alone. Within the regression model, we further modeled the potential confounding effects of age, sex, level of depressive symptoms (BDI total score), total number of experienced non–child abuse trauma types, and genetic ancestry (by incorporating the ancestry estimates based on analysis of the ancestry informative markers). Due to missing phenotype and genotype information on some participants, a maximum number of 676 individuals and a minimum of 633 individuals were informative in the regression models for the interaction effects. The individuals in the group with 676 and 633 that actually entered the analyses were not different in age, sex, self-described race, income, employment status, relationship status, child abuse trauma exposure, as well as PTSD and depression severity from the remainder of the total 900 individuals. They did, however, have higher child abuse rates. In the individuals entering the analyses, the no abuse group size ranged from 67.5% to 68.6% (and were not different) vs 77.1% to 79.0% in the excluded individuals (P = .01 to P = .14).

Permutation Analyses . As the non-normality of PSS and possible sparse-cell counts in the FKBP5-trauma interaction strata could each invalidate the asymptotic P values produced by regression models, we established the significance of main environment, genotype, or genotype-trauma interaction effects using permutation procedures that randomly assigned the sample PSS scores to participants (sampled without replacement), while holding each participant's genotype and environmental variables fixed.^79- 80 For each analysis, we based inference on 10 000 to 100 000 permutations. We conducted these analyses by using appropriate components of SAS version 9.1 (SAS Institute, Cary, North Carolina). Due to smaller cell sizes, we also used these permutation-based analyses to establish the main genotype or genotype and PTSD interaction effects on the DST.

Correlation Analyses . For descriptive and quantitative examination of genotype × continuous environment (CTQ total score) interaction of PTSD symptom severity, we used partial correlation controlling for age and sex, correlating the semiquantitative child abuse measures of the CTQ total score with the current continuous PSS score in the whole sample as well as the sample stratified by FKBP5 SNP genotype. The significance of the interaction of FKBP5 SNP genotype and CTQ total scores on PSS was established using permutation-based methods, and the post hoc genotype-dependent differences of the correlation coefficients were established by converting the correlation coefficients of each genotype group to Fisher Z scores by Z = ln[(r+1)/(r-1)]/2. We then estimated the standard error of difference between the 2 correlations by SE = {[1/(n₁-3)] + [1/(n₂-3)]}^1/2 and divided the difference between the 2 Z scores by the standard error. If the Z value for the difference was 1.96 or higher, the difference in the correlation was significant at P = .05 level. If the difference was 2.58 or higher, the difference in the correlation was significant at P = .01 level.⁸¹

Correction for Multiple Testing . To ensure an overall significance level of P ≤ .05 for the primary analyses (main FKBP5 SNP effects on PTSD symptoms, interaction of FKBP5 SNPs with child abuse or non–child abuse trauma exposure on adult PTSD symptoms), we corrected for the examination of 24 tests (8 tests of the main effects of the FKBP5 SNPs, 8 tests of interaction between each SNP and non–child abuse trauma exposure, and 8 tests of interaction between each SNP and categorical measure of child abuse). We used the conservative Bonferroni method for such a multiple-testing correction, which yielded a significance threshold of α = .002. In an exploratory approach, we also genotyped 3 CRHR1 SNPs based on previous results from our group in examining depression.⁷¹ When we corrected for this additional set of 9 tests (3 CRHR1 SNPs, main effect and interaction of SNP × child abuse or non–child abuse trauma exposure), the α level after correction was .05/33 = .0015.

Non–Child Abuse Trauma Exposure and PTSD Symptoms . To analyze the differential roles of gene × environment interaction with PTSD in adult participants, we first had to demonstrate that within our participant population, the level of total non–child abuse traumatic exposure contributes significantly to current PTSD symptoms. To perform these analyses, we first used a general linear model (controlling for age and sex), using the PSS total score as the dependent variable and the TEI categorical variable representing number of types of non–child abuse traumatic experiences as the independent variable. We found a significant sex effect (F_1,818 = 29.5, P < .001) and a very robust non–child abuse trauma exposure effect (F_3,818 = 61.9, P < .001) (Table 4). In the zero types of non–child abuse trauma group, the mean (SD) PSS score was 3.58 (6.27) and in the 4 or more types of trauma group, we found continuous increases by more than 5-fold to 16.74 (12.90).

Child Abuse Trauma Exposure and PTSD Symptoms . We next examined whether child abuse exposure also predicted the level of current PTSD symptomatology (Table 5). We performed the same analyses as above with the categorical child abuse variable (none, 1 type [physical or sexual], and 2 types [physical and sexual]). Similar to the effect of non–child abuse trauma, we also found a significant effect in the presence of child abuse (F_2,806 = 50.9, P < .001). The experience of child abuse increased the mean (SD) PSS scores from 8.03 (10.48) in the no child abuse group to 14.65 (11.92) in the 1 type and 20.93 (14.32) in the 2 types of child abuse group.

Non–Child Abuse Trauma Exposure, Child Abuse, and PTSD Symptoms . Using a general linear model that included both non–child abuse trauma and child abuse trauma (Table 6) and their interaction term as predictors of adult PTSD symptoms, we observed significant main effects of the 2 terms (non–child abuse trauma: F_3,806 = 16.2, P < .001; child abuse trauma: F_2,806 = 14.4, P < .001) and sex (F_1,806 = 20.2, P < .001) as described above, but no significant interaction (F_6,806 = 0.51, P = .79) between the 2 main predictors. The presence of 1 type of child abuse increases the PSS score at each non–child abuse trauma level by 4.58 points on average (range, 2.03-6.51) and the presence of 2 types of child abuse increases the PSS score at each non–child abuse trauma level by 9.04 points on average (range, 7.79-10.20), suggesting that exposure to child abuse increases risk for higher levels of PTSD symptoms in response to non–child abuse trauma exposure.

Main Effect of FKBP5Polymorphisms on PTSD Symptoms . Based on previous data supporting a clear role of FKBP5 in modulating the glucocorticoid response to stress, as well as evidence supporting association of FKBP5 variants with risks for and rate of recovery from affective disorders,⁵² we hypothesized that genetic variation at FKBP5 may influence liability to PTSD. We examined 8 SNPs spanning 120 kb of the FKBP5 locus (Table 3) but found no significant main effect of FKBP5 genotypes on PSS total score (Figure 1).

Interaction of FKBP5Polymorphisms With Non–Child Abuse Trauma to Predict PTSD Symptoms . We next examined whether FKBP5 polymorphisms interact with increasing levels of non–child abuse (primarily adult) trauma to predict adult PTSD symptoms. We regressed PSS scores on FKBP5 genotype, non–child abuse trauma (using the 4-level variable), and the interaction between genotype and non–child abuse trauma, adjusting for age and sex. We did not find any significant interaction effects using permutation-based methods (Figure 1). These data demonstrate that FKBP5 polymorphisms do not appear to have a simple role in moderating the effects of non–child abuse trauma exposure on PTSD outcomes.

Interaction of FKBP5Polymorphisms With Child Abuse to Predict PTSD Symptoms . We then regressed PSS scores on FKBP5 genotype, child abuse (3-level variable), and the interaction between genotype and child abuse, adjusting for age and sex. In this analysis, we identified 4 SNPs that showed significant interactions with child abuse after correcting for multiple testing (interaction P < .002) (Figure 1). All 4 associated SNPs are in fairly high linkage disequilibrium, with pair-wise r² values ranging from 0.53 to 0.87 (Figure 1). The most significant SNP, rs9296158, is located in intron 5 with an interaction P < .0004. Two of the associated SNPs, rs3800373 and rs1360780, have previously been reported to associate with differential glucocorticoid-mediated responses by Binder et al.⁵⁷ For all 4 significant SNPs, we observed a similar additive mode of interaction, and there appears to be a gene dose-dependent protection from severe child abuse–associated increases in adult PTSD scores (Table 7). In the group of individuals with 2 types of child abuse, individuals homozygous for the protective G allele of rs9296158 had a mean (SEM) PSS score of 13.54 (3.76); in the heterozygous group, the mean PSS score was 21.25 (2.03); and in the group homozygous for the risk allele A, the mean PSS score was 31.11 (5.37). The PSS scores of more than 20 indicate clinical significant PTSD and higher scores indicate more severe PTSD, with a maximum of 51 points that can be reached in that scale.^61- 63 Interestingly, the genotypes previously associated with a higher number of previous depressive episodes and faster response to antidepressant treatment⁵⁷ were the genotypes with the highest level of adult PTSD symptoms in the presence of child abuse.

FKBP5and Child Abuse Interaction . Given the prior data on FKBP5 and major depression, and the observation that in our sample the 3-level child abuse variable is significantly and positively correlated with our 4-level non–child abuse trauma variable (n = 851, Pearson correlation = 0.271; controlling for age and sex, P < .001), we reexamined the above interaction controlling for BDI score as a continuous measure of depression symptoms or the 4-level non–child abuse trauma variable in addition to age and sex. After adjusting for these variables, the interaction terms of child abuse and FKBP5 SNPs remained significant for 2 of the 4 SNPs (rs9296158 and rs3800373; interaction P = .07 for rs1360780 and interaction P = .06 for rs9470080). In addition, we also conducted a regression analysis with BDI as the dependent variable. None of the FKBP5 SNPs showed a significant interaction of genotype and child abuse to predict depression outcome (controlling for age and sex) after controlling for multiple testing. Thus, these data suggest that the FKBP5 genotype may interact with childhood trauma to predict adult PTSD severity, even when controlling for age, sex, non–child abuse trauma level, and severity of depressive symptoms.

Genotype-Dependent Correlations of Child Abuse and PSTD Symptom Severity . We explored the genotype-dependent correlations of child abuse and PSTD symptom severity interaction effect using an ordinal variable of child abuse, which contains 154 points ranging from a CTQ score of 25 to 120. This allowed us to examine correlation analysis without having to group individuals into only 3 child abuse levels, which can lead to small group sizes. There is also importance in validating these effects with a separate, broader, and more continuous measure of child abuse. Thus, we investigated the interaction of the genotypes of the 4-associated FKBP5 SNPs with the CTQ total score to predict PTSD symptom severity. Permutation-based analyses showed an interaction effect with rs9296158 genotype and the continuous environment measure with P = .01 and P = .03 for rs3800373 and P = .07 for rs1360780. To better describe and quantify these differences, we reinvestigated positive correlations of child abuse severity and PSS and BDI scores stratifying by the genotypes of the 4-associated FKBP5 SNPs. Partial correlation (controlling for age and sex) of CTQ total score and PSS scores yielded highly significant correlations of child abuse and PTSD symptom severity in the whole sample (r = 0.347, P = 3.7 × 10⁻²⁷). However, when stratifying this analysis by the rs9296158 genotypes (Table 7), we observed significantly higher correlation coefficients in the AA vs GG group (z value for difference of correlation coefficient = 3.00, P < .01) and the AG vs GG group (z value for difference of correlation coefficient = 2.35, P < .05) for correlations between CTQ total score and PSS score. Similar effects were observed with 2 other SNPs, rs1360780 and rs3800373 (Table 8). Interestingly, although BDI scores and the CTQ total scores also showed highly significant positive correlations (R = 0.39, P = 6.9 × 10⁻³³), there were no significant differences of the correlation coefficients among the 3 rs9296158 genotype groups.

Genetic Admixture and Population Stratification. Based on statistical analysis of the ancestry informative marker genotype data from the subsample, we found that sub-Saharan African ancestry estimates for participants ranged between 41% and 100%, with a mean of 83.6%; European ancestry estimates ranged between 0% to 54%, with a mean of 10.9%; Indigenous American ancestry estimates ranged between 0% to 19%, with a mean of 3.3%; and East Asian ancestry estimates ranged from 0% to 33%, with a mean of 2.3%. We found that the degrees of ancestry from these subpopulations were not significantly correlated with total PSS score, level of non–child abuse trauma exposure or presence of child abuse (all Pearson correlation coefficients between 0.01 and −0.04). In addition, using the percentage ancestry estimates for these 4 populations as covariate in the interaction analysis did not alter the observed SNP × environment interactions. This evidence suggests that population stratification is not likely confounding in our analyses.

Effects of CRHR1SNPs on PTSD Symptom Severity . Because we have recently reported interactions of CRHR1 SNPs with child abuse on adult depression,⁷¹ we tested the 3 SNPs with the strongest interaction effects (rs110402, rs242924, and rs7209436) for the effect on PTSD symptom severity. None of the 3 SNPs showed a significant main effect or interaction effect with child abuse or non–child abuse trauma exposure on PTSD symptom severity.

DST, PTSD, and FKBP5SNP . To assess whether FKBP5 SNPs had a functional consequence on GR sensitivity changes observed in PTSD, we examined the interaction of probable PTSD diagnosis and the 4 FKBP5 SNPs with significant gene × environment interactions on the change in cortisol concentration before and after 0.5-mg dexamethasone. These data were available for 80 individuals and, for any given genotype, we had full genotype, phenotype, and DST data for 75 to 78 individuals. In this sample, the mean (SD) serum concentrations of cortisol on day 1 were 12.08 (5.06) μg/dL (to convert to nmol/L, multiply by 27.588) and following dexamethasone, mean (SD) concentrations were 3.65 (4.38) μg/dL. When we examined the change from baseline cortisol to the dexamethasone-suppressed cortisol concentration using repeated measures analysis of variance, we found significant interactions between risk allele carrier status and PTSD categorical diagnosis (probable diagnosis based on PSS inventory), with repeated measure of cortisol (rs9296158: F_1,78 = 8.76, P < .004; rs3800373: F_1,76 = 6.03, P = .02; rs1360780: F_1,79 = 3.95, P = .05; and rs9470080: F_1,77 = 5.32, P = .02) (Figure 2), but no significant main effects. These data suggest that the majority of the patients with the risk alleles with PTSD show enhanced suppression to dexamethasone or enhanced GR sensitivity, which has been reported to be a possible endocrine signature of PTSD. In contrast, those individuals with the putative resilience genotype with PTSD appear to show the opposite effect. These results are supported by permutation-based analyses using percentage suppression as outcome and FKBP5 risk allele carrier status and probable PTSD as predictors.

Our results indicate that levels of child abuse and non–child abuse trauma each independently predicted the level of adult PTSD symptomatology. Although polymorphisms in FKBP5 did not directly predict the level of PTSD symptoms or interact with the level of non–child abuse trauma to predict PTSD symptoms, SNPs within the FKBP5 locus robustly interacted with the level of child abuse to predict the level of adult PTSD symptoms.

The most novel and important finding of our study was the interaction between FKBP5 polymorphisms and child abuse history to predict the levels of adult PTSD symptoms. The polymorphisms seem to belong to the same bin of SNPs, all in high linkage disequilibrium, which is associated with functional differences in FKBP5 expression and consequent alterations in GR function. Notably, all 4 SNPs showing a significant interaction effect had either been reported to be associated with higher FKBP5 protein levels or a stronger induction of FKBP5 mRNA by cortisol in healthy probands (rs1360780 and rs3800373),⁵⁷ or were in strong linkage disequilibrium with these SNPs (rs9296158 and rs9470080). The SNP genotypes that were associated with the highest FKBP5 mRNA induction in peripheral blood mononuclear cells by cortisol⁵⁷ were also the ones that were associated with the highest vulnerability to PTSD symptoms following child abuse. Individuals carrying the other allele seemed to be protected from the development of PTSD symptoms in a gene-dose dependent manner. This is in agreement with the finding of Segman et al⁵⁹ who showed that trauma-induced increased levels of FKBP5 mRNA expression in peripheral blood mononuclear cells immediately following medical trauma were predictors of the presence of PTSD at 4 months after the trauma.

FKBP5 expression is induced by glucocorticoids as part of an intracellular ultrashort negative feedback loop for GR activity,⁵³ with increased expression of FKBP5 reducing glucocorticoid binding affinity⁵⁴ and nuclear translocation of the GR,⁵⁵ resulting in resistance to glucocorticoid activation. Thus, the alleles previously associated with high FKBP5 protein/mRNA expression⁵⁷ should be associated with GR resistance. This is precisely what we observed in individuals without PTSD symptoms. The healthy carriers of these alleles showed less dexamethasone suppression and thus more GR resistance.

This functional association appears to be switched in patients with PTSD symptoms. These same alleles were associated with a higher dexamethasone suppression ratio and thus increased GR sensitivity, which is associated with PTSD,^43,45 while the protective genotype was associated with relative GR resistance in patients with PTSD symptoms. This environment-dependent reversal of the functional association may be similar to the effects previously reported in patients with depression, where less cortisol response in the dexamethasone-CRH test (an indication of increased GR sensitivity) was observed with the same alleles that had been associated with more FKBP5 protein and thus presumably GR resistance in healthy controls.⁵⁷

Our study is the first to our knowledge to provide evidence directly supporting a developmental or symptom-dependent difference in the functional consequences of these FKBP5 SNPs on GR sensitivity. We hypothesized that specific FKBP5 alleles may enhance the effects of acutely released cortisol following child abuse on FKBP5 mRNA expression and that abnormal FKBP5 expression leads to maladaptive changes in GR sensitivity. These changes resulted in long-term alterations of HPA axis sensitivity, such as GR hypersensitivity that effect adult response to trauma. Consistent with this notion, alterations of HPA axis responsiveness^43,45 have been previously identified as risk factors for PTSD.

PTSD has been suggested to result in part from initial overconsolidation of traumatic memories^82- 83 or, conversely, abnormal extinction of such memories.^8,84 Thus, alterations in FKBP5 function could conceivably be involved in abnormal GR-mediated signaling in neurons involved with stress response and memory formation. Polymorphisms within the FKBP5 gene that lead to altered GR responsiveness could promote sensitization of the neural systems involved in stress response and emotional memory processing, thereby placing children who have been abused with specific genetic variants at higher risk for PTSD or PTSD-spectrum symptoms when exposed to other types of traumas. This hypothesis may be supported by the finding that the alleles of rs3800373 and rs1360780, which are associated with higher risk for PTSD symptoms following severe child abuse in our study, were also associated with higher levels of peritraumatic dissociation in children after medical trauma.⁵⁸ Notably, peritraumatic dissociation has previously been identified as another risk factor for the development of PTSD.¹⁰

It is presently not clear if the associated polymorphisms represent the actual functional variants that lead to the differential FKBP5 expression patterns or are in linkage disequilibrium with a so far untyped, potential functional variant. Denser fine-mapping of this region combined with resequencing and in vitro and in vivo functional studies may allow definitive identification of the genetic variants mechanistically responsible for the interactions and associations observed herein. The fact that in our study rs4713916, located in a potential regulatory region 5′ upstream of FKBP5, did not show a significant interaction effect, although in Binder et al⁵² it had similar association patterns to rs1360780 and rs3800373, is likely due to the different extent of linkage disequilibrium observed among these SNPs in the present sample of black individuals vs the previous sample of ethnic German individuals (linkage disequilibrium was substantially stronger in the German sample).

Depression and PTSD show a high comorbidity,⁸⁵ so that concurrent depressive symptoms could confound our genetic interaction analysis for PTSD symptom outcome. However, controlling for current depressive symptoms, which were measured with the BDI, did not alter the observed interaction effect on PTSD symptoms and there was no gene × child abuse effect observed when BDI score was used as the outcome. Although the PSS and BDI scores strongly and significantly correlate in our sample (Pearson correlation = 0.68, P < .001), our observation of a PTSD symptom-specific effect suggests that interaction of FKBP5 genotypes with child abuse history might influence PTSD symptoms not captured by the BDI, and thereby not shared with depression. Additional evidence suggesting separate gene × environment interactions for depression and PTSD symptomatology also comes from the lack of an interaction effect between SNPs in the CRHR1 gene and child abuse on PTSD symptom severity, despite the fact that these SNPs show strong interaction effects on the symptom severity of adult depression.⁷¹

One major limitation of our study is that we cannot present an independent replication, so that at this level of validation, it can only be considered hypothesis generating. The supporting evidence from functional effects of the polymorphisms in the DST serve to strengthen our finding that awaits independent replication. Another limitation is the lack of ancestry informative marker data on the complete sample. We do, however, feel the existing ancestry results among our sample of 280 participants genotyped for the 134 ancestry informative markers help to alleviate the concern of confounding by population stratification. The ancestry estimates in our subsample genotyped for the ancestry informative markers are not correlated with PSS or BDI scores.⁷¹ Because confounding due to population stratification only occurs when there is both correlation between outcome and ancestry as well as correlation between SNP genotypes and ancestry, our results are not likely susceptible to bias arising from confounding due to population structure.

Limitations in the phenotype side are that we used retrospective reports of child abuse. It is well-documented that reports of trauma are correlated with PTSD symptoms,^86- 87 and we have to acknowledge that the retrospective reports of child abuse and also non–child abuse trauma can therefore be biased in favor of an association between these reports and PTSD symptoms. A similar limitation is the lack of examination of the time since the index trauma. Furthermore, because we have used PTSD symptom severity and not diagnosis as an outcome measure and our sample has very high levels of multiple trauma exposure, our results may not be comparable with those obtained in studies with PTSD diagnosis as outcome and studies with less variation in trauma type or fewer overall types of trauma exposure. However, several studies have now indicated that multiple trauma exposure across the lifespan is the rule rather than the exception in samples similar to our study^3,88- 89 (eg, low income, urban, high percentage of black or Latino participants) and that multiple trauma exposure is related to increased mental and physical health risk.^65,90- 91 Understanding risk and resilience in response to multiple trauma exposures is of high public health importance.

Our sample was also derived from primarily impoverished, inner-city outpatient primary care clinics and participants were not presenting for treatment for PTSD so that we cannot directly generalize our findings to a clinical or epidemiological setting. Comparability with other studies and populations may be further limited as methods of assessment and definitions of child abuse vary greatly across studies. Finally, the mechanisms accounting for the relationship of child abuse to increased risk for PTSD in adulthood are not well understood. It may be that child abuse directly impacts psychological and biological developmental processes or that child abuse is associated with other variables (child/parental temperament, broader family environment, attachment) that impact psychological and biological development or an interaction of the 2.

Data reported herein suggest that genetic variation in the FKBP5 gene, which is involved in glucocorticoid signal transduction, may alter sensitization of the stress-response pathway during development, placing those individuals who have had significant child abuse at significant risk for PTSD in the face of other traumatic experiences. These genotypes potentially serve as predictors of both risk and resilience for adult PTSD among survivors of child physical and sexual abuse.

Our genetic results support the hypothesis that the glucocorticoid response system moderates the effects of early life stress on adult PTSD symptoms and that GR hypersensitivity may be important in the pathophysiology of this disorder. These results suggest the possibility that heritable differences in glucocorticoid-mediated neural functioning exacerbate or dampen the effects of child abuse on the stress hormone system, thus altering HPA axis sensitivity and risk for PTSD in adulthood.