|Home | About | Journals | Submit | Contact Us | Français|
Pituitary adenylate cyclase-activating polypeptide (PACAP) is known to broadly regulate the cellular stress response. In contrast, it is unclear if the PACAP/PAC1 receptor pathway has a role in human psychological stress responses, such as posttraumatic stress disorder (PTSD). In heavily traumatized subjects, we find a sex-specific association of PACAP blood levels with fear physiology, PTSD diagnosis and symptoms in females (N=64, replication N=74, p<0.005). Using a tag-SNP genetic approach (44 single nucleotide polymorphisms, SNPs) spanning the PACAP (ADCYAP1) and PAC1 (ADCYAP1R1) genes, we find a sex-specific association with PTSD. rs2267735, a SNP in a putative estrogen response element within ADCYAP1R1, predicts PTSD diagnosis and symptoms in females only (combined initial and replication samples: N=1237; p<2x10−5). This SNP also associates with fear discrimination and with ADCYAP1R1 mRNA expression. Methylation of ADCYAP1R1 is also associated with PTSD (p < 0.001). Complementing these human data, ADCYAP1R1 mRNA is induced with fear conditioning or estrogen replacement in rodent models. These data suggest that perturbations in the PACAP/PAC1 pathway are involved in abnormal stress responses underlying PTSD. These sex-specific effects may occur via estrogen regulation of ADCYAP1R1. PACAP levels and ADCYAP1R1 SNPs may serve as useful biomarkers to further our mechanistic understanding of PTSD.
PACAP was first isolated from ovine hypothalamic extracts based on its ability to stimulate cAMP production in anterior pituitary cells1. It is a highly conserved member of the VIP/secretin/glucagon peptide family with pleiotropic functions in development, cell signaling, metabolism, homeostasis and cell protection2–6. Among these myriad functions are studies demonstrating 1) high expression of PACAP peptide and its selective PAC1 receptor in hypothalamic and limbic structures, 2) PACAP regulation of corticotropin releasing hormone (CRH) and autonomic function, 3) actions of PACAP in stress-related behavior, 4) reduced anxiety-like phenotypes in PACAP and PAC1 receptor null mice, and 5) blunted corticosterone response in knockout animals after emotional stressors. Thus, PACAP/PAC1 receptor signaling is integrally involved in stress mechanisms7,8. We hypothesized that PACAPergic systems may be important mediators of abnormal stress responses following psychological trauma contributing to PTSD, an extreme maladaptive and debilitating psychiatric disorder affecting up to 40% of individuals over lifetime exposure to traumatic events9,10.
Little is known about the biological processes regulating PTSD and other stress-related responses. To examine whether the PACAP/PAC1 pathway is associated with PTSD in a high risk, heavily traumatized population, we analyzed blood levels of PACAP, and genetic variation and methylation of the PACAP (ADCYAP1) and PAC1 receptor (ADCYAP1R1) genes in a cohort of over 1200 highly traumatized subjects with and without PTSD (See Suppl. Tables 1–2 for demographic information).
Using radioimmunoassay, we first examined PACAP peptide levels in peripheral blood samples from a previously described, highly traumatized, at risk population11–13 that had been matched on age, sex, and trauma histories (n = 64, see Suppl. Tables 1–3 for demographics). We found that PTSD symptoms (PTSD Symptom Scale14) were significantly correlated with PACAP38 (PACAP peptide containing 38 residues) blood levels in females (p < 0.005, r = 0.497, Fig 1A), but not in males (p>0.5). Also in females, PTSD diagnosis was associated with PACAP38 levels (p ≤ 0.001), with higher PACAP38 found in the PTSD cohort. Furthermore, PACAP levels (median split, low vs. high) were differentially associated with PTSD symptoms in females (Fig 1B). PACAP38 levels also predicted differential response on all three symptom clusters necessary to fulfill diagnostic criteria for PTSD (intrusive re-experiencing (e.g., trauma flashbacks), avoidance (e.g., avoidance of trauma reminders) and hyperarousal (e.g., increased startle response)) in females but not males (Fig 1C). These analyses were repeated in a second, all female cohort (N=74) with similar findings (Fig 1D; high vs. low PACAP38 levels, controlling for age, substance abuse and total trauma exposure, one-tailed t-tests: total symptoms, p ≤ 0.05, hyperarousal symptoms, p ≤ 0.001; and % with clinically significant symptoms, χ2=4.9, p<0.05). These observations were especially notable since females may be at twice the risk for PTSD as compared to males 9,11, implicating roles for sex hormones, especially estrogen, in the disorder15–17. When we controlled for common stress-related phenomena (depression and history of substance abuse), the effect of PACAP level on PTSD remained (p<0.05). In contrast, there was no effect of PACAP level on depression symptoms or substance abuse when controlling for PTSD.
In addition to the psychological symptoms that define the syndrome, subjects with PTSD have been found to have abnormally high conditioned fear responses. This high level of fear may result from a combination of an inability to habituate to aversive stimuli, a decreased ability to extinguish (learn to inhibit) fear memories, and possibly an over-consolidation of the original fear memory18–22. Hence, we determined the physiological (electromyographic) levels of conditioned fear for 27 participants (16 male, 11 female) with PACAP blood levels. Fear potentiation was determined by measuring the acoustic startle reflex in the presence of startle cues alone, or startle cues combined with stimuli paired (conditioned stimulus, CS+) or unpaired (CS−) with an aversive airblast (Fig 1E). Female (but not male) subjects with high PACAP38 levels demonstrated markedly increased startle reflex responses to both CS+ and CS− cues. This was particularly pronounced during the late acquisition phase when normal subjects had habituated to the fearful stimuli (Fig 1F). In aggregate, these data suggest that PACAP38 peptide is strongly associated with the psychological and physiological symptoms of PTSD in women with a history of trauma.
To assess whether there may be a genetic association of PTSD with polymorphisms in either the PACAP (ADCYAP1) or PAC1 receptor (ADCYAP1R1) locus, we performed a tag-SNP analysis (r2 = 0.8; MAF = 0.1) across both genes with a total of 44 SNPs (14 ADCYAP1 and 30 ADCYAP1R1 SNPs). Using logistic regression, we examined whether each SNP was associated with PTSD diagnosis in this cohort of highly traumatized urban civilian subjects (n = 798)11,12,23, in total, or stratified by gender (females: n = 503; males: n = 295). Only the ADCYAP1R1 receptor SNP rs2267735 (p = 0.0002 in females; NS in males) remained significant after experiment-wide multiple correction for sex and 44 independent tests (Figure 2A–B and Suppl. Fig. 1). No SNPs in the peptide ADCYAP1 gene met experiment-wide criteria for association (Suppl. Fig. 2). Given these striking gender differences and recent data demonstrating that ADCYAP1R1 gene expression may be dynamically modulated by estrogen24, the distribution of estrogen response elements (EREs) within ADCYAP1R1 gene was examined (Suppl. Table 4). We found that rs2267735 was within a predicted ERE (Fig 2C, Genomatix; matrix similarity = 0.877, core similarity = 1.0). Because rs2267735 is positioned within the central variable region of the consensus sequence, in silico analyses do not currently allow us to predict how the ‘C’ vs ‘G’ allele may differentially alter ERE function and further in vitro analyses are warranted.
We next determined if the association between rs2267735 and PTSD diagnosis could be replicated in an additional 439 subjects. These subjects were from the same overall study, but were interviewed and had DNA collected after the original discovery population. Thus they served as a replication source from the same population but distinct in time and with different interviewing staff. The table in Figure 3A shows the logistic regression results for males and females separately in the initial population described in the tag-SNP analysis, the replication sample from the same population, and the combined sample of 1237 individuals. The main effect of the SNP on PTSD diagnosis could be replicated in women (p<0.05) and combining both samples increased the significance of the association (N = 763, p< 0.00002). As in the discovery sample, no effects were observed in males (male combined sample N = 474, p = 0.7).
To further examine ADCYAP1R1 rs2267735 SNP associations with continuous PTSD symptom levels in females, we analyzed both an additive and a dominant model with total PTSD symptoms and symptom subscales using the combined samples (Fig 3B – E). The ‘CC’ allele was most robustly associated with total PTSD symptoms and among subscales, hyperarousal symptoms were the most strongly associated with rs2267735. Notably, even after controlling for childhood trauma history and adult trauma, age and race, (which slightly reduces total N due to missing data) the rs2267735 ‘CC’ allele was associated with higher levels of PTSD hyperarousal symptoms compared to ‘G’ carriers in women (p = 0.0008, Fig 3E), but not men (p=0.51).
We repeated the above analyses with Beck Depression Inventory (BDI) symptoms and history of lifetime substance abuse, and found no associations with these measures and rs2267735 (Suppl. Fig. 3), suggesting that this association may be relatively specific to PTSD. To address whether rs2267735 might be associated with other severe psychiatric illnesses, we performed analyses using bipolar disorder, schizophrenia, and Alzheimer’s disease samples. From the data of the Genetic Association Information Network (GAIN) publicly accessible database (http://www.ncbi.nlm.nih.gov/projects/gap), we analyzed the association of rs2267735 (included on the Affymetrix 6.0 SNP array) with bipolar disorder as well as schizophrenia. We did not observe a significant association of this SNP with these two disorders in subjects with African American (954 cases, 1195 controls) or European (1378 cases, 1351 controls) ancestry. Specifically, we found that all pre-computed p-values for associations of rs2267735 with schizophrenia or bipolar disorder were higher than the multiple-testing correction p-value of 0.01, indicating no major contribution of this variant.
Additionally, we examined the association of rs2267735 and Alzheimer’s in a previously characterized Alzheimer’s disease sample25. In this cohort of 342 subjects, we found no association with rs226735 and Alzheimer’s disease diagnosis using either the additive genetic model (p=0.19) or the dominant/recessive model (p=0.89). These data suggest that we find robust associations with rs2267735 in women, but not men, with PTSD. In contrast, we find no association with depression symptoms, substance abuse, Alzheimer’s disease, bipolar disorder, and schizophrenia across different samples. Note that for all of these negative results, due to the limited sample sizes, we cannot rule out the possibility that rs2267735 may be associated with PTSD in men or with other disorders with a smaller effect size than we see with PTSD in women.
To parallel our results with plasma PACAP38 levels, we next examined whether physiological measures of fear are differentially associated with the ADCYAP1R1 rs2267735 SNP. In PTSD, but not depression18, fear response to an inhibitory CS−, or ‘safety signal’, is exaggerated. The discrimination between CS+ and CS− improves across the training procedure in controls, but not in those with PTSD. We examined whether rs2267735 was associated with impaired fear discrimination late in conditioned acquisition, during the same period noted in Figure 1E. Notably, females with the ‘CC’ genotype were significantly less able to discriminate CS+ from CS− signals (Figure 3F, sex by genotype interaction, p<0.05, and ‘CC’ vs ‘G’ carriers in females, p<0.05).
We next examined whether a difference in dark-enhanced startle, a measure of increased anxiety in humans that is similar to light-enhanced startle in rodents26–28, was differentially associated with rs2267735. Again, we found that females, but not males, with the ‘CC’ genotype showed significantly more startle in the dark compared to the light (Fig 3G, males, N=35, p=0.71; females, N=53, p=0.02). Together, these data suggest that the ADCYAP1R1 rs2267735 SNP may be relatively specific in its association with PTSD psychological and physiological phenotypes. Further, the robust association of rs2267735 with hyperarousal symptoms suggests that the role of PACAP/PAC1 may be specifically involved in the normalization of the stress response, a process which is particularly dysregulated in PTSD.
Environmental, genetic and epigenetic mechanisms likely moderate the long-term effects of trauma exposure. Using the Illumina HumanMethylation27 BeadChip, we interrogated methylation in DNA extracted from peripheral blood at the first site within the ADCYAP1R1 CpG island (Suppl. Fig. 2). Methylation at this site was significantly associated with total PTSD symptoms (Figure 4A, N=94, r = 0.354, p<0.0005) in a sex-independent manner. Further, CpG methylation level (median split) was associated with PTSD diagnosis (Figure 4B, χ2=8.1, p < 0.005), but not depression (p > 0.05, Suppl. Fig. 3E). There was no significant association between methylation of ADCYAP1 and PTSD symptoms. These data suggest that ADCYAP1R1 is regulated, in part, through epigenetic mechanisms that contribute to differential function of the PAC1 receptor in PTSD.
To examine the potential relationship of genotype and brain mRNA expression as described previously29, we utilized a brain mRNA expression data set30 to test whether ADCYAP1R1 rs2267735 is associated with differential gene expression. We first examined whether cortical ADCYAP1R1 and ADCYAP1 mRNA levels were correlated. As shown in Fig 4C, these mRNA levels were significantly inversely correlated (r = −0.219, p < 0.001, including males and females), suggesting that brain levels of PACAP peptide and PAC1 mRNA are tightly regulated.
We next utilized a prior analyzed genome-wide association and brain mRNA expression data30 to examine whether ADCYAP1R1 rs2267735 imputed genotypes were associated with differential ADCYAP1R1 expression in brain. We found a sex x genotype effect (Fig 4D, F(3,99) = 4.3, p < 0.05) with females with the ‘CC’ genotype expressing significantly less ADCYAP1R1 mRNA than males(F(1,33) = 5.5, p < 0.05) or than females who are ‘G’ carriers (one-tailed, F(1, 45)=2.87, p<0.05). Thus, mRNA encoding the PACAP peptide and PAC1 receptor appeared to be tightly regulated within the human cortex, and ADCYAP1R1 mRNA levels were associated with the ADCYAP1R1 rs2267735 SNP.
Despite prior studies examining PACAP/PAC1 receptor function in CNS/PNS development, endocrine homeostasis, metabolism, cellular protection/regeneration, and chronic stress responses2,3,6,31–34, a role for PACAP signaling in fear conditioning has not been evaluated. Given our data implicating PACAP in PTSD, we wondered if ADCYAP1R1 mRNA was differentially regulated in mice using Pavlovian fear conditioning35–40, a means of studying acute fear and trauma responses that has been proposed to model PTSD19,22. We performed classical fear conditioning experiments using male mice, in which a previously neutral tone CS (6kHz) was paired with 10 footshocks (1mA, 0.5sec; Figure 5A). This conditioning paradigm consistently provides robust fear learning in mice leading to changes in gene expression within the amygdala, a region critical for fear learning and expression. Quantitative PCR analyses shows that amygdala ADCYAP1R1 mRNA increased ~1.5-fold during the consolidation of fear (Figure 5B, p<0.05), with a similar trend within the mPFC. When peak freezing was compared with brain mRNA levels, we find a significant correlation between fear learning and ADCYAP1R1 mRNA (Fig 5C, r2=0.49, p<0.05).
To further establish the relationship between PACAP/PAC1 receptors and estrogen in a validated model of sex hormone regulation, we examined estrogen-induced changes in ADCYAP1 and ADCYAP1R1 transcripts in the bed nucleus of stria terminalis (BNST) in female rats. The BNST is a component of the extended amygdala that is subject to significant gonadal hormonal control7,27,28. In rodents, it is critical for emotional behavior, mediating stress responses and the light-enhanced startle response. We examined gene expression in the BNST in ovarectomized (OVX) female rats following 21-day implantation of continuous release estrogen pellets. Compared to control implants, estradiol increased ADCYAP1 transcripts in the dorsal and ventral BNST 2.1- and 3.4-fold, respectively (p ≤ 0.01, Figure 5D). Additionally, estradiol increased ADCYAP1R1 mRNA 1.5-fold in the dorsal BNST samples (p<0.05, Figure 5E), and future studies should also examine estradiol sensitivity of these genes in amygdala. While these rodent studies are complex and have differing experimental designs, these data clearly illustrate dynamic PACAP/PAC1 receptor regulation within central areas mediating fear, stress and estrogen responsiveness.
Since its identification more than 20 years ago, PACAP has been increasingly implicated in diverse cellular stress response pathways and neurotrophic function. However, the organizational role of the PACAP system in orchestrating behavioral stress responses has yet to be clarified. Our data suggest that PACAP/PAC1 receptor expression and signaling may be integrally involved in regulating the psychological and physiological responses to traumatic stress. Further, we report association of an ERE-embedded ADCYAP1R1 SNP with PTSD, and we demonstrate fear- and estrogen-dependent regulation of PACAP systems within stress-responsive regions of the brain. These data may begin to explain sex-specific differences in PTSD diagnosis, symptoms, and fear physiology. Future work targeting the PACAP/PAC1 receptor system may lead to novel and robust biomarkers as well as to further our understanding of the neural mechanisms underlying pathological responses to stress with potential therapeutic targets towards the prevalent and debilitating syndrome of PTSD.
This highly traumatized, civilian, cross-sectional cohort has been previously described in candidate gene-association studies of PTSD and depression11–13. Research interviews, salivary DNA and blood samples were collected from patients receiving services in the primary care clinics at Grady Memorial Hospital (Atlanta, GA). All study procedures have been reviewed and approved by the Emory Institutional Review Board and the Grady Hospital Research Oversight Committee. PTSD measures in this manuscript are based on the PTSD symptom scale 14, which has been validated within this population using the Clinician Administered PTSD Scale. PACAP38 radioimmunoassay (1:30,000, Peninsula Laboratories, Belmont, CA) was performed at University of Vermont, using double antibody immunoprecipitation as previously described 41. For genotyping, pairwise tagging (R2 > 0.8, MAF > 0.1) was used to choose tag-SNPs for both ADCYAP1 and ADCYAP1R1. The coordinates were chr18:885000-906000 and chr7:31048667-31117836 for ADCYAP1 and ADCYAP1R1, respectively (NCBI B36) which includes approximately 10kb upstream and 5kb downstream of the coding regions for both genes. Genotypes for the tag-SNPs were generated using Sequenom iPlex with follow-up analyses using Taqman. For methylation analyses, DNA was whole-genome amplified, fragmented, and hybridized to the HumanMethylation27 BeadChip (Illumina, San Diego, CA). Individual samples were stratified to separate BeadChips according to PTSD status to limit bias. The BeadChips were scanned using a BeadStation 500GX, and the methylation level (beta value) was calculated using the Methylation Module of the BeadStudio software. The eyeblink component of the acoustic startle response was measured by EMG recordings of the right orbicularis oculi muscle with two 5-mm Ag/AgCl electrodes filled with electrolyte gel, as previously described18–19. The mouse fear conditioning and rat estrogen replacement studies are described in detail in Supplementary Methods.
This work was primarily supported by National Institutes Health Grants MH071537, DA019624 (KJR) and HD27468 (VM). Support was also received from Emory and Grady Memorial Hospital General Clinical Research Center, NIH National Centers for Research Resources (M01RR00039 and P20RR16435), the American Foundation for Suicide Prevention (BB) and the Burroughs Wellcome Fund (KJR). RNA samples were run using the Oncogenomics Core Facility at University of Miami. Dr Myers is supported by The National Institute on Aging (AG034504). We thank Drs. Joseph F. Cubells, Yilang Tang, and Karen Conneely for excellent discussions. We thank the Grady Trauma Project, including Drs. Charles F. Gillespie, Ann Schwartz, Aliza Wingo, David A. Gutman, and Tamara Weiss for medical support; Allen Graham, Angelo Brown, Justine Phifer, Daniel Crain, Asante Kamkwalala, James Poole, Dorthie Cross, Negar Fani, and Dr. Ami Smith for clinical research support; and Kristin Schutz, Emily Reiser and Caitlin Fitzgerald for excellent molecular/genetics technical support. Methylation chip assays were performed by the Emory University Biomarkers Service Center.
Supplementary Online Material: Detailed Materials and Methods (13 pages), with 4 Tables and 3 Figures
Reprints and permissions information is available at www.nature.com/reprints.
The authors declare no competing financial interests.
Author Contributions: K.J.R. designed the experiments and wrote the paper. K.J.R., B.B., and V.M. organized collaborations, obtained funding, supervised data collection and analyses, and revised the paper. K.B.M., K.K., and E.B.B. performed the genetics experiments and analyses on the primary, replication, and GAIN cohorts, and revised the paper. T.J. and S.D.N. performed and supervised the human physiology studies and revised the paper. V.K. and A.K.S. performed and supervised the methylation studies and revised the paper. A.M., S.E.H., D.T., K.J.R., and V.M. performed and supervised the rodent fear conditioning and estrogen replacement studies and rtPCR analyses and revised the paper. A.J.M., M.R., and A.E. performed and supervised the human mRNA expression analyses and the Alzheimer’s disease gene association studies. K.M.B. and V.M. performed and supervised the PACAP38 radioimmunoassays and related data analyses.