Skip to main content

Advertisement

Region- and time-dependent gene regulation in the amygdala and anterior cingulate cortex of a PTSD-like mouse model

The Correction to this article has been published in Molecular Brain 2019 12:47

Abstract

Posttraumatic stress disorder is developed by exposure to a threatening and/or a horrifying event and characterized by the presence of anxiety, hyperarousal, avoidance, and sleep abnormality for a prolonged period of time. To elucidate the potential molecular mechanisms, we constructed a mouse model by electric foot shock followed by situational reminders and performed transcriptome analysis in brain tissues. The stressed mice acquired anxiety-like behavior after 2 weeks and exaggerated startle response after 4 weeks. Avoidance latency and freezing behavior were sustained up to 5 weeks post stress and abnormal static behavior was observed during the sleep period. RNA sequencing was performed in two of the emotional regulatory regions, anterior cingulate cortex and amygdala, at 2 and 5 weeks post stress. More than 1000 differentially expressed genes were identified at 2 weeks in both regions. The number of the regulated genes remained constant in amygdala at 5 weeks post stress, whereas those in anterior cingulate cortex were plummeted. Although synaptic remodeling and endocrine system were the most enriched signaling pathways in both anterior cingulate cortex and amygdala, the individual gene expression profile was regulated in a region- and time-dependent manner. In addition, several genes associated with PTSD involved in Hypothalamic-Pituitary-Adrenal axis were differentially regulated. These findings suggested that global gene expression profile was dynamically regulated in accordance with the disease development stage, and therefore targeting the distinct signaling molecules in different region and development stage might be critical for effective treatment to PTSD.

Introduction

Post-traumatic stress disorder (PTSD) is a psychiatric disorder defined by profound disturbances in cognition, emotion, behavior and physiological function in response to a life threatening or a horrifying traumatic event. The symptoms of PTSD are characterized by three long-lasting symptomatic clusters (DSM-5): (1) re-experiencing (intrusive thoughts, recurrent nightmares, flashbacks and reactivity to reminders of the trauma); (2) avoidance and emotional numbness (avoiding stimuli associated with trauma and anhedonia); and (3) hyperarousal (hypervigilance, exaggerated startle response and sleep disturbance). PTSD patients may also exhibit impaired cognition, working memory, and reduced social activities. Another characteristic of PTSD is that those symptoms emerge long after traumatic events in some patients, suggesting that disease can develop after an extended asymptomatic period.

Anxiety disorders including PTSD encompass heterogeneous conditions that are associated with abnormalities in brain regions controlling fear response such as the medial prefrontal cortex, amygdala (AMY), hippocampus, and nucleus accumbens [1, 2]. Changes in activity and functional connectivity of these brain regions lead to the development and maintenance of anxiety disorders [3,4,5]. AMY is a critical brain region involved in the detection of threat, fear learning and expression, and heightened emotional memory. The exaggerated activity in AMY following traumatic events is correlated with the severity of symptoms in PTSD patients. AMY is modulated by the cortical influence from medial prefrontal cortex. The altered activity in the anterior cingulate cortex (ACC), a subregion of medial prefrontal cortex, has been demonstrated in PTSD using functional MRI [6]. With extensive connection to limbic structures including AMY, ACC is thought to be the direct top-down regulator of PTSD susceptibility through modulating AMY activation to threatening stimuli [7, 8]. In particular, dorsal ACC is thought to facilitate the expression of conditioned fear and increased neuronal activity in PTSD during fear conditioning [2, 9]. The volume of ACC was consistently found to decrease in PTSD, however some studies have revealed hyperactivation of the ACC rather than hypoactivation [10]. Therefore, it is important to investigate the functional activity of ACC in relation to AMY hyperactivation in PTSD.

Recent genomic profiling and neuroimaging studies for human PTSD have uncovered several candidate genes for the pathogenesis, however detailed investigation on gene regulation and the underlying mechanisms are still limited [11,12,13]. Molecular and cellular encoding of traumatic events and behavioral responses can be reflected by the changes of synaptic plasticity. Specifically, signaling molecules associated with synaptic transmission and plasticity in connecting to AMY are implicated as primary substrates for fear learning and memory. Alteration in glutamatergic and GABAergic neurotransmission contributes to abnormal behavioral responses [14]. In addition, downstream signaling components, such as MAPK, cAMP and calcium mediated signaling molecules are involved in synaptic plasticity [15]. The alteration in several neuromodulatory systems such as endocannabinoid signaling, Wnt signaling and ErbB signaling can affect fear memory [16] and anxiety-like behavior by synaptic remodeling in medial prefrontal cortex [17] and AMY [18, 19].

Hypothalamic-pituitary-adrenal (HPA) axis plays a pivotal role in stress response upon exposure to emotional or physical stressors, and its feedback loop mediated by cortisol is altered in PTSD [20, 21]. In fact, human genetic studies identified several PTSD-associated genes are involved in HPA axis and cortisol, for instance FKBP5 and corticotropin releasing hormone receptors [22]. Oxytocin and thyroid hormone signaling pathways are modulated by cortisol and the imbalance may cause emotional lability, impatience, anxiety and sleep disturbances [23, 24].

As PTSD diagnostic criteria are not met at least 1 month after traumatic event according to DSM-5, sustained psychological abnormality is one of the prominent features. However, the molecular mechanisms underlying these psychological changes are not well-addressed. In this study we employed genome-wide unbiased gene analysis using RNA-seq to characterize the regulation of differentially expressed genes in a stress induced PTSD-like mouse model and found that several signaling pathways involved in synaptic remodeling and endocrine systems, as well as PTSD-associated genes were regulated in a region- and time-dependent manner.

Materials and methods

Animals

Male, 8 to 10 weeks old C57BL/6J mice were purchased from Jackson Laboratories (Bar Harbor, ME). Animals were maintained under standard conditions of controlled temperature (23 °C) and housed in a 12 h light/dark cycle (light on at 6 AM) with access to water and food ad libitum. They were group housed with 3 to 5 mice per cage and habituated to the vivarium for at least 3 days. Experimental procedures were carried out in accordance with NIH guidelines and approved by the Uniformed Services University Animal Care and Use Committee. One hundred twenty-four mice in total were used for performing various behavioral tests (Fig. 1) and molecular analyses.

Fig. 1
figure1

Timeline of PTSD-like model construction and behavior tests. C57BL/6J mice at 8–10 weeks of age were subjected to acclimation followed by electric foot shock. The control mice underwent the same acclimation procedure without foot shock. Avoidance test and situational reminder (AV/SR) were weekly performed, and the freezing test was done at 5 weeks post stress (PS). Open field test (OF) and elevated plus maze test (EPM) were performed at 2 weeks and 4 weeks PS. Light/dark box test (LDB) was performed around 4 weeks PS. Acoustic startle response test (ASR) was performed around 2 weeks and 4 weeks PS, and locomotive activity test were carried out from 4 weeks to 5 weeks PS

Foot shock induction

Animal model for PTSD-like symptoms was constructed by foot shock induced fear conditioning and re-exposure to the shock chamber as previously described [25, 26]. Briefly, trauma procedure consisted of a single exposure to an inescapable electric foot shock (Day 0) followed by re-exposure to the shock chamber every week until 4 weeks (Fig. 1). Mice received 15 s handling for 3 consecutive days before given foot shock. On the shock day, each animal was placed in the light compartment. After a 2 min adaptation period, the door to the dark shock apparatus made of plexiglass (17 cm × 17 cm × 25 cm) was open to let the animal entering the shock chamber. After another 5 min adaptation in the closed chamber, the animal received inescapable foot shock (1.5 mA × 2 s with 2 min interval for 8 times). The control group received the same procedures without stress. Animals were removed from the apparatus 1 min after the last shock.

Behavioral test

All the behavioral tests were performed during the light on phase from 8 AM to 3 PM. Mice were given one test a day to avoid the potential interference from the other tests (Fig. 1). Animal cages were moved to a testing room at least 1 h prior to each test. After completion of the test session, the behavioral apparatus and chamber were cleansed with 70% ethanol, and then dried by hand-fan completely.

Avoidance latency/situational reminder and freezing tests

For measuring avoidance latency, animals were placed in a compartment adjacent to a foot shock chamber. After 30 s acclimation, the door was open and the time spent before entering foot shock chamber was measured within a period of 2 min. Test was terminated once mouse entered, and the mouse was immediately removed from the shock chamber to a new cage. For measuring freezing time, animals in 5 weeks post stress were placed inside chamber with the same context, and mouse behavior was monitored for 5 min with the digital video system. The data was analyzed by ANY-maze software (Stoelting Co, Wood Dale, IL) to determine freezing time and the number of episodes.

Open field test

Each mouse was placed in the center of a cubic chamber (40 cm × 40 cm × 35 cm) with a black Plexiglas wall under dim light in the testing room. The position and movement of animals were monitored for 10 min by video camera installed above the chamber. Data analyses were carried out by ANY maze software to determine the moving velocity, total travel distance, the time spent in the center zone, the entry number to the center zone. Center zone was assigned in a 20 cm × 20 cm square area in the center of arena.

Elevated plus maze test

The apparatus consisted of four arms (34 cm × 5 cm width), with two arms open and two closed by gray walls (15 cm height) arranged in the opposite side of the same type. The platform was located 50 cm above the floor of the testing room illuminated and four arms were connected in the center platform (5 cm × 5 cm), where animal was placed facing to a closed arm. The position and movement of animal was monitored for 5 min by video camera. An entry was defined if more than half of the animal’s full body entering the open arm. The time spent in the open arms and number of visits to the open arms were analyzed using ANY-maze.

Light/dark box test

The light–dark box was made of white and black opaque Plexiglas (40 × 20 × 30 cm light chamber, 23 × 20 × 30 cm dark chamber). The chambers were connected by a 6 × 6 cm door in the middle of the wall separating the two chambers. Animals were placed in the middle of the dark chamber with the door closed. Video monitoring was started when the door was open and continued for 10 min to measure transition number into the lit area, retention time in lit, and travel distance inside lit using ANY-maze software.

Acoustic startle test

Startle response was measured using ventilated startle chambers (SR-LAB system, San Diego, CA). Sound levels within each chamber were measured using a sound-level meter to ensure consistent presentation. Animals were placed in a Plexiglas cylinder (4.7 cm × 10 cm long) inside a chamber (29 cm × 29 cm × 29 cm). The cylinder resting on a platform was connected to a piezoelectric accelerometer to monitor movement. Each test session was started with a 3 min acclimation period to background white noise of 70 dB, followed by 10 times acoustic stimulus (120 dB × 20 ms) with a 1 min interval. Animals were returned to home cage immediately after test session. Data analyses were carried out by SD software (SR-LAB system) and summarized with Microsoft Excel. The latency to first peak, maximum amplitude, and average startle amplitude from 1 ms to 100 ms were analyzed.

Home-cage immobility test

Video, magnet, photobeam, or infrared sensor system has been used to measure home cage activity [27]. This system has been applied for monitoring sleep-wake behavior instead of polysomnography [28, 29]. Group housed animals were separated individually into a cage (20 cm × 38 cm) and nourished with water gel and food ad libitum for 3 days. One day prior to the test, the cages were placed on a test platform before 6 PM to adapt the environment. The room light was turned off at 6 PM and turned on at 6 AM the following day, in accordance with the previous light-dark cycle. Locomotive activity was monitored immediately after 6 AM with the video camera system and continued for 11.5 h. Data was analyzed numerically by ANY maze software to determine the immobile period. Locomotive activity during the entire test period was dissected with a 2 min interval to assess immobile period which showed no locomotion within each 2 min-section. Total travel distance and total immobile time were measured by ANY maze. Short immobile episode was determined as the number of immobile periods of 1 or 2 sequential immobile sections.

RNA extraction

Mice were sacrificed after deep anesthetization from 10 AM to 2 PM. Brains were immediately fresh frozen on dry ice and stored at − 80 °C. Coronal section of the brain at 300 μm was cut at − 20 °C. Three adjacent brain sections of bilateral amygdala (AMY) (Bregma − 0.94 mm to − 1.84 mm) and those of bilateral anterior cingulate cortex (ACC) (Bregma + 1.34 mm to + 0.44 mm) [30] were collected by punching with 1 mm-diameter brain punch based on the method described by Palkovits [31]. Total RNAs from AMY and ACC were isolated by Trizol (Invitrogen, Carlsbad, CA) and purified with the RNeasy Mini Kit (Qiagen, Germantown, MD). The RNA concentration was measured by Qubit 3.0 Fluorometer (Invitrogen).

RNA-seq

Library preparation for RNA-seq was performed by TrueSeq stranded mRNA sample preparation kit (Illumina, San Diego, CA) starting from 600 ng of total RNA. Sequence data were generated by Illumina 3000 and NextSeq platforms. We converted raw BCL files to FASTQ files and performed multiplexing using Illumina bcl2fastq2 software (version 2.17). Sequencing quality control was done by Fastqc program. We aligned the resulting FASTQ files to mm 10 mouse reference genome by STAR aligner. Quantification to read counts, FPKM and TPM was done by RSEM program. With raw read counts matrix as input, we identified genes that were differentially expressed by DESeq2 R package, using unpaired two-class significance analysis and a false discovery rate threshold of 0.05. Gene expression visualization was plotted by heatmap R package. KEGG Pathway and GO analysis were performed by the Database for Annotation, Visualization and Integrated Discovery (DAVID) v6.8 [32] or “DOSE” R Bioconductor package [33].

Statistics

The GraphPad Prism 7 (GraphPad Software Inc., San Diego, CA) was used for statistical analysis. The weekly avoidance latency test was analyzed by two-way ANOVA followed by Bonferroni’s multiple comparison test. The statistical comparison between the two groups was determined by Welch’s unpaired two tailed t-test. All data are expressed as mean ± SEM. Statistical significance was set at p < 0.05. *, **, and *** depicting p < 0.05, p < 0.01, and p < 0.001, respectively.

Results

Sustained avoidance and fear memory in foot shock stressed mice

For assessment of avoidance behavior, latency till moving into the same foot shock chamber from the light compartment was measured every week. The latency of stressed mice was longer than that of control mice starting from the first week and sustained for at least 4 weeks (F4,104 = 2.782, p = 0.0305) (Fig. 2a). Notably, almost all stressed mice did not move out from light compartment to the stressed chamber except for the first trial (Fig. 2a). When placed inside the same shock chamber at 5 weeks after stress, the stressed mice possessed a longer freezing time compared to the control mice (p = 0.0038) (Fig. 2b), although the frequency of freezing episode was not significantly altered.

Fig. 2
figure2

Latency to avoidance and the freezing behavior test. The latency to enter the foot shock chamber was tested weekly and the stressed mice had a prolonged latency compared to the control mice (a). At 5 weeks PS, the freezing behavior inside the chamber was monitored for 5 min and analyzed by ANY maze software (b). Data are represented as mean ± S.E.M. (n = 14/group), Two-way ANOVA with Bonferroni’s test in (a) and unpaired two-tailed t-test in (b) were performed

Anxiety-like behavior emerged at 2 weeks post stress

As anxiety is a typical behavioral feature in PTSD, we examined if stressed mice acquired and sustained anxiety using several behavioral tests. Open field test was performed in a black walled squire space for 10 min. At 2 weeks PS, the stressed mice showed decreased retention time in the center zone and entry number to the center zone (p = 0.0056) (Fig. 3a, b). Total travel distance was significantly shorter in the stressed mice than control mice (p = 0.0006) (Fig. 3c). At 4 weeks PS, the stressed mice showed significant reduction in retention time in the center zone (Fig. 3d, p = 0.0126), entry number to the center zone (Fig. 3e, p = 0.0304), and total travel distance (Fig. 3f, p = 0.0047). Shorter travel distance was unlikely caused by motor deficit of the stressed mice, since maximum running speed was not impaired (Fig. 3g). In the elevated plus maze test, the stressed mice spent less time in open arms than control mice at 2 weeks PS, despite the significant difference was not reached (Fig. 4a). The entry number to open arms was significantly lower in the stressed mice (p = 0.0095) (Fig. 4b). When examined at 4 week PS, there were no significant differences in both time spent and entry number in open arms between the stressed mice and control mice (Fig. 4c and d). To further examine the anxiogenic behavior at these time points, we also performed a light/dark box test. It was shown that the number of transitions to light area and travel distance in light area during 10 min test period were significantly reduced in the stressed group than the control group (Fig. 5a and c). Although time spent in light area was not significantly different (Fig. 5b), the stressed mice preferred dark protected area and less explored in an unprotected area compare to control mice at 4 weeks PS. Overall, these tests demonstrated that the stressed animals started to show anxiety-like behavior at 2 weeks PS and continued till 4 weeks PS.

Fig. 3
figure3

Anxiety-like behavior assessed by Open field test. Mice were placed in a cubic square under dim light and allowed to move freely for 10 min. Mice movement related to anxiety was assessed based on the time in center zone (a and d), entry number to center zone (b and e), total travel distance (c and f), and maximum speed (g). Panels a, b, and c show the results at 2 weeks PS, panels d-g show those at 4 weeks PS. Unpaired two-tailed t-test was performed. Data are represented as mean ± S.E.M. (n = 24/group at 2 weeks PS, n = 34/group at 4 weeks PS)

Fig. 4
figure4

Anxiety-like behavior assessment by elevated plus maze test, mice were placed in a plus maze apparatus and allowed to move freely for 5 min. Time in open arms (a and c) and entry number to open arms (b and d) were measured at 2 weeks PS (a and b) or 4 weeks PS (c and d). Unpaired two-tailed t-test was performed. Data are represented as mean ± S.E.M. (n = 24/group at 2 weeks PS, n = 34/group at 4 weeks PS)

Fig. 5
figure5

Anxiety-like behavior assessment by light/dark box test. Mice at 4 weeks PS were freely moved in an arena divided by light and dark areas connected with a small door for 10 min. Animal movement in terms of transition frequency (a), retention time in lit area (b), and travel distance in lit are (c) were analyzed. Unpaired two-tailed t-test was performed. Data are represented as mean ± S.E.M. (n = 16/group)

Hyperarousal symptoms in stressed mice

Acoustic startle response against 120 dB stimulus repeated for ten times were monitored by a motion sensor system, and the maximum peak amplitude, latency to the first peak, and average amplitude were measured. As shown in Fig. 6a, both control and stressed mice showed a similar pattern to startle response. The average amplitude of the maximum peak out of 10 trials was not different between control and stressed mice at 2 weeks PS. The average amplitude from 1 ms to 100 ms was also similar despite a slight increase observed in the stressed mice (Fig. 6c). However, at 4 weeks PS the peak amplitude was significantly higher in the stressed mice (p = 0.0048) (Fig. 6d), and the average amplitude was also greater in stressed mice than control mice (p = 0.0440) (Fig. 6e). There was no difference in startle response latency between the stressed and control mice at either test period. These results suggested that startle response was exaggerated in stressed animals at relatively later time point.

Fig. 6
figure6

Acoustic startle response test for control and stressed mice. Mice were acclimated in the chamber under 70 dB background noise for 3 min, followed by start acoustic stimulus (120 dB × 20 msec) 10 times in 1 min interval. Representative startle response profiles were shown in (a) for control as broken line and for stressed mice as solid line. The average of 10 trials of peak response amplitude was measured at 2 weeks PS (b) and 4 weeks PS (d). The average amplitude from 1 ms to 100 ms was measured at 2 weeks PS (c) and 4 weeks PS (e). Unpaired two-tailed t-test was performed. Data are represented as mean ± S.E.M. (n = 24/group at 2 weeks PS, n = 34/group at 4 weeks PS)

Increased frequency of short immobile episodes during the sleep period

Sleep disturbance is a common feature and impacts the pathophysiology of PTSD and other psychiatric diseases [34]. We characterized the sleep-wake behavior of stressed mice with long-term video monitoring in their home cages. There was no significant difference between control and stressed mice in total immobile time during the entire inactive phase (Fig. 7a). However, the number of short immobile episodes was significantly larger in the stressed mice than the control mice (p = 0.0353) (Fig. 7b), suggesting that the sleep/wake cycle was disturbed in stressed mice. The total travel distance showed no significant difference between the control and stressed mice, indicating the locomotive activity was not impaired in stressed animals (Fig. 7c).

Fig. 7
figure7

Characterization of immobile pattern during inactive phase. Locomotive activity in home-cage was monitored by video camera during inactive phase for 11.5 h and analyzed using ANY-maze software. Total immobile time was shown in (a), and number of short immobile episode (immobile duration between 2 min and 4 min) was shown in (b), and total travel distance was shown in (c). Data are represented as mean ± S.E.M. (n = 28/group). Unpaired two-tailed t-test was performed

Differentially expressed genes in ACC and AMY

We performed genome wide RNA expression analysis in the anterior cingulate cortex (ACC), which is involved in the modulation of emotional and fear encoding and consolidation. Of the 3233 differentially expressed genes (DEGs) in ACC at 2 weeks post stress (PS), 1335 were upregulated, and 1898 were downregulated (Additional file 1: Figure S1A). However, at 5 weeks PS the numbers plunged to 43 (Additional file 1: Figure S1B) and only 16 of them were overlapped between these time points (Additional file 1: Figure S1C). In AMY 2 weeks PS, there were 1063 DEGs, in which 873 genes were predominantly down-regulated (Additional file 2: Figure S2A). Different from ACC, there were a substantial number of DEGs in 5 weeks PS (Additional file 2: Figure S2B). Of more than 1000 DEGs, approximately 9% (94 DEGs) were overlapped between 2 weeks and 5 weeks PS as shown in the Venn diagram (Additional file 2: Figure S2C).

Regulation of PTSD associated genes in the mouse model

Several human genomic studies have been conducted recently to identify candidate genes for PTSD patients [35,36,37,38]. We found that ten genes associated with PTSD were regulated in either AMY or ACC (Table 1). Several of these genes were involved in HPA axis and glucocorticoid regulation, which include Crhr1 and 2, Fkbp5, Rora, Sgk1, and Stat5b. Other genes including Adcy8, Cacna1c, Klhl1, and Shank1 were involved in synapse plasticity and regulation. Gene expression regulation in our animal model was correlated, at least to some extent, with that identified in different brain areas or blood samples from patients with PTSD.

Table 1 Regulation of PTSD-associated genes

Differentially expressed genes in AMY

Those DEGs in AMY were applied for Gene Ontology (GO) analysis using DAVID to determine enriched functional annotation in terms of molecular functions (Additional file 3: Table S1) or biological processes (Additional file 4: Table S2). First, to characterize the regulation of global mRNA expression, we searched for the GO terms related to transcriptional regulation and found comparable numbers of DEGs in the Transcription and Negative regulation of transcription between 2 weeks and 5 weeks PS, whereas the Positive regulation of transcription was enriched only in 5 weeks PS (Additional file 5: Figure S3). Of note, the Covalent chromatin modification was enriched only in 2 weeks PS. Consistently, there were four more enrichments related to the chromatin modifications: unmethylated CpG binding, histone deacetylase binding, histone H3-K4 methylation and histone H4-K5 acetylation.

There were several gene regulations in terms of cellular alteration in the central nervous system (CNS) (Additional file 6: Figure S4). Nervous system development was enriched both in 2 weeks and 5 weeks PS, while neuron projection development was only seen in 2 weeks PS. Neurogenesis and Neural precursor cell proliferation were enriched in 5 weeks PS. GOs related to neuron differentiation were enriched both in 2 weeks PS and in 5 weeks PS. For non-neuronal cell differentiation, Oligodendrocyte differentiation and Epithelial cell differentiation were enriched in 2 and 5 weeks PS, respectively.

Gene regulation in AMY

A number of KEGG pathways were involved in synaptic remodeling (Fig. 8) from the enriched pathway analyses (Additional file 7: Table S3). Pathways for neurotransmitters including acetylcholine, dopamine, GABA, glutamate, and serotonin were enriched in 5 weeks PS. We found gene regulation of glutamatergic kainite receptor subunits, Grik3, 4 and 5, metabotropic receptor subunits, Grm4 and 7. In GABAergic synapse, type A GABAergic receptor subunits (Gabrg1, Gabrg3, and Gabrq) and GABA transporter (Slc38a1) were regulated after 5 weeks. Moreover, acetylcholine esterase, Ache, and muscarinic cholinergic receptor, Chrm3, dopamine receptor1, Drd1, and serotonergic receptors, Htr1b, 1d, 2c, and 7 were regulated.

Fig. 8
figure8

Enriched KEGG pathways related to synapse remodeling in AMY. Enriched pathways were categorized into 3 functions, involved in neurotransmitter (upper panel), synapse mediated signaling (middle panel), and synapse formation and activity (lower panel). Each pathway in the two time points was represented as heatmap with color index based on p-value

In terms of synapse mediated signaling, LTP and LTD pathways and several downstream signaling pathways crucial for synapse plasticity such as MAPK, Ras, and PI3K/Akt were enriched. Ras signaling pathways including Mras, Kras, as well as Ras GTPase regulatory factors, and Rap1 signaling pathways including Rap1 GTPase regulatory factors were overlapped in both 2 weeks and 5 weeks PS. And the key molecules in calcium mediated signaling, calmodulin (Calm2), CAM kinase (Camk2a), and CREB binding protein (Crebbp) which are involved in LTP and LTD were also altered. The key components of Wnt signaling, four receptor genes and three ligands including β-catenin, were regulated. Notably, the retrograde endocannabinoid signaling was highly regulated at 5 weeks PS. Moreover, pathways related to synaptic formation and activity, including Adherens junction, Axon guidance, Focal adhesion, Endocytosis, and Regulation of actin cytoskeleton pathways were enriched.

To further characterize the gene regulation of synapse remodeling, the pathways overlapped in 2 weeks and 5 weeks PSs were characterized by individual DEGs. Most of DEGs were regulated either after 2 weeks or 5 weeks PS, but there were a few DEGs overlapped in both time periods, particularly involved in the Endocytosis pathway (Table 2). These results suggest that individual DEG was differentially regulated at 2 weeks and 5 weeks PS. Consistent with the differential regulation of synapse remodeling genes in AMY, DEGs related to transcription showed only a few overlapped genes between 2 and 5 weeks PS, suggesting a diverse transcriptional regulation.

Table 2 DEGs in four representative overlapped pathways and transcription in AMY

GO analysis and pathway analysis revealed that several GOs and pathways were involved in hormone and neuroendocrine, particularly in 5 weeks PS but not at 2 weeks PS (Fig. 9), among them are Oxytocin, Thyroid hormone, Gonadotropin releasing hormone, and Estrogen signaling pathways. Oxytocin (Oxt) was regulated at both 2 and 5 weeks PS, while oxytocin receptor (Oxtr) was only in 2 weeks PS. Thyroid hormone receptors and transporter (Thra, Rxrg, and Slc16a2) were regulated in 5 weeks PS.

Fig. 9
figure9

Enriched KEGG pathways related to hormone and neuroendocrine system in AMY. Nine pathways that were enriched based on DEGs in AMY were represented as heatmap with color index based on p-value

Gene regulation in ACC

There were substantial number of GOs and pathways enriched in 2 weeks PS. As seen in AMY, we found several GOs and KEGG pathways related to hormone and neuroendocrine system and synapse remodeling. Figure 10 shows enriched KEGG pathways including Estrogen signaling, GnRH signaling, Insulin signaling, Oxytocin signaling, Prolactin signaling, and Thyroid hormone signaling. Substantial number of KEGG pathways were also involved in synapse remodeling. MAPK signaling pathway was the only pathway seen in 5 weeks PS among them. Thus, synapse remodeling and hormone and neuroendocrine system were the mainly regulated pathways in 2 weeks PS, and they were mostly terminated in 5 weeks PS.

Fig. 10
figure10

KEGG pathways enriched in ACC 2 weeks and 5 weeks PS. Substantial number of DEGs in ACC showed enrichment in hormone and neuroendocrine system or synapse remodeling. The enriched pathways from the DEGs were represented in heatmaps, and these pathways included hormone & neuroendocrine (upper panel), neurotransmitter (upper middle panel), synapse mediated signaling (lower middle panel) and synapse formation & activity (lower panel). Color index represents the level of significance (p-values)

Gene regulation related to PTSD-like behavior in AMY and ACC

GOs associated with PTSD-like behavior such as fear response, startle response, social behavior were predominantly seen after 5 weeks in AMY and 2 weeks in ACC (Fig. 11). In AMY at 5 weeks PS, Behavioral fear response did not manifest a strong enrichment, though Startle response, Response to drug, Sensory perception of pain, Memory, Learning, and Social behavior were significantly enriched. Relating to the Response to drugs, several GOs were selected including Response to cocaine, amphetamine, morphine, nicotine and ethanol (data not shown). On the other hand, Social behavior and Sleep were listed in 2 weeks PS, despite with a weak enrichment (p > 0.05). In ACC at 2 weeks PS, Sensory perception of pain, Social behavior, Memory, and Learning were enriched. Overall, GOs associated with PTSD-like behavior were regulated selectively in AMY at 5 weeks PS and in ACC at 2 weeks PS.

Fig. 11
figure11

GOs related to PTSD-like behavioral traits. Several GOs in AMY (left columns) and ACC (right columns) were associated with behavioral and cognitive abnormality in PTSD, including fear memory, startle response, drug comorbidity, pain sensitivity, memory and learning, social interaction, and sleep. Color index represents level of significance (p-values)

Discussion

Behavioral characterization of the foot shock induced stress model

Electric foot shock has been widely used as fear conditioning, because it can provide consistent stress to animals and is easy to assess fear memory. Severe foot shock has also been used as a stressor for PTSD-like animal model [39]. In previous reports, animals receiving electric foot shock was subsequently re-exposed to the same contextual environment [40]. This procedure which is called situational reminder seems critical to acquire long-term PTSD-like abnormal behaviors including long-lasting hyperarousal [26], sustained freezing behavior [25] as well as higher evoked field potential in basolateral AMY [41]. In this study, animals were given the weekly avoidance test and we observed that, except for the first week, all the stressed animals avoid to enter the previously foot-shock chamber, suggesting these stressed animals were threatened under this situation. Together with freezing test at 5 weeks PS, we found that the stressed animals acquired sustained fear memory for more than 1 month. Since anxiety-like behavior is a prominent characteristic of PTSD, we carried out three different behavioral tests (Figs. 3, 4, and 5) and found that the stressed mice acquired anxiety-like behavior, which emerged after 2 weeks and sustained until 4 weeks PS. In contrast, exaggerated acoustic startling response of the stressed mice was only acquired after 4 weeks PS (Fig. 6). Thus, two abnormal behavioral traits in our model were differentially evolved. Moreover, short immobile episodes were more frequent in stressed mice than control animals, suggesting a sleep disturbance in the stressed animals and a reminiscence of the abnormal sleep-wake cycle in PTSD patients.

Differential susceptibility of animals to stress has been recognized as an important trait of the animal model because of the existence of human resilience or vulnerability to trauma [42, 43]. In fact, we observed a wider sample variance in stressed group than in control group in the acoustic startle test at 4 weeks PS (Fig. 6c and e), but not at 2 weeks PS (Fig. 6b and d), suggesting that differential susceptibility was created when stress induced hyperarousal was prominent. Cohen’s group has evaluated individual rodent susceptibility based on the acoustic startle response [44]. Currently it is unknown whether the wide variation was derived from individual difference occurred during the development of PTSD-like abnormality, since there is lack of consistency for the stressed animals to behave as vulnerable or resilient in different behavioral assays.

Gene regulation associated with PTSD

Genome-wide association studies of human PTSD have identified the disease specific polymorphism of several associated genes in their gene loci. Some of the associated genes are involved in HPA axis and glucocorticoids to regulate stress response and emotional control [13]. In our model, several genes were regulated in either AMY or ACC, including Crhr1 and 2 [45,46,47], Rora [48], Sgk1 [49], and Stat5b [50] (Table 1). Among them, FKBP5 has been intensively investigated since it is a co-chaperone for glucocorticoid (GC) receptor to suppress its transcriptional activity [51]. In addition, we found genes involved in synaptic plasticity: Adcy8, whose deficiency led to stress-induced anxiety [52] implying a role in long-term potentiation and synaptic plasticity for fear related memory [38], Cacna1c (Cav1.2 subunit of L-type calcium channel), which mediates long-term synaptic plasticity through calcium signaling, Klhl1, the actin binding protein modulating calcium current [53] and neurite outgrowth [54], and Shank1, the anchoring protein playing a role in glutamatergic synaptic plasticity [55]. These results suggest that those PTSD-associated genes can play an important role on the disease development.

Time-dependent gene regulation in AMY

Numerical gene regulation was comparable between 2 weeks and 5 weeks PS since there were more than 1000 DEGs in each condition (Additional file 2: Figure S2). However, DEGs in the overlapped pathways were predominantly in either 2 weeks or 5 weeks PS, but only a few DEGs overlapped, consistent with the transcription related DEGs (Table 2). Alteration of drastic transcriptional regulatory genes resulted in dynamic changes of global gene expression including synapse remodeling or hormone/neuroendocrine system. In addition, histone mediated epigenetic regulation may contribute to time-dependent regulation (Additional file 5: Figure S3), as recently postulated that epigenetic alteration is deeply involved in the pathogenesis as well as in the disease susceptibility [43, 56, 57]. Of note, the glutamatergic, GABAergic and endocannabinoid signaling were all highly enriched in 5 weeks PS (Fig. 8), consistent with the notion that the excitatory and inhibitory neurotransmission is tightly regulated by the endocannabinoid signaling under stressful conditions [58]. In addition, cholinergic and serotonergic input can modulate neurotransmission at 5 weeks PS. Ras, Rap1, and MAPK pathways which play an important role in synaptic transmission, including fear memory encoding and expression [59, 60] and calcium mediated signaling, including LTP were also predominantly regulated (Fig. 8). Moreover, hormone and neuroendocrine system were also predominantly enriched in 5 weeks PS (Fig. 9). However, we did not find significantly regulated expression of the core HPA axis genes including corticotropin releasing hormone, its receptor, or GC receptor, which was different from a previous study [61]. It is possible that the differences might be derived from the different animal species or experimental procedures. PTSD-like model constructed by foot shock followed by situational reminder may not clearly provoke HPA axis activation in male animal as reported previously [62]. Instead, we found that oxytocin signaling, which is known to exert anxiolytic and antidepressant effect [63, 64], and also as an important negative modulator of HPA axis and GC signaling [65], was highly enriched in 5 weeks PS (Fig. 9). In fact, oxytocin was upregulated in 2 weeks PS, then plummeted in 5 weeks PS. Oxytocin receptor gene was also upregulated in 2 weeks PS, while returned to basal level in 5 weeks PS. Thus, there is a possibility that extra-hypothalamic GC signaling is modulated by regulatory neuroendocrine system such as oxytocin, although further investigation is necessary to clarify the regulation of HPA axis. Moreover, several GOs that are related to the PTSD-like behavioral traits such as startle response, response to drug, and sensory perception to pain were enriched at 5 weeks PS (Fig. 11), suggesting that those DEGs play a key role in the molecular mechanism to disease development. Taken together, gene regulation involved in synapse remodeling and hormone/neuroendocrine system in AMY is likely culminated in 5 weeks PS, which might be correlated to the development of PTSD-like abnormal behaviors.

Termination of gene regulation after 5 weeks in ACC

ACC is reported to alter the neuronal activity in PTSD patients during fear conditioning test as well as in the resting condition [66]. ACC has extensive afferent connections to emotion regulatory limbic regions including AMY and is involved in emotional learning and certain social behavior [67]. Given its proximity and connection to the limbic structures, ACC is thought to be the direct top-down regulator of PTSD susceptibility through modulating AMY activation to threatening stimuli [7, 8]. Despite that the volume of ACC is reduced in PTSD, some studies have shown hyperactivation of the ACC rather than hypoactivation [10]. Thus, regulation of AMY activation as well as PTSD development by ACC remains to be enigmatic. Our gene analyses in the region found that more than 3000 genes were regulated at 2 weeks PS, then plunged to 43 DEGs at 5 weeks PS, which suggests that gene regulation was mostly terminated at 5 weeks PS (Fig. 6). In addition to synapse remodeling and neuroendocrine and hormone pathways, 18 DEGs related to oxidative phosphorylation were all upregulated specifically in 2 weeks PS (data not shown), suggesting that neuronal activity as well as gene regulation was accelerated together with enhanced mitochondrial activity. Furthermore, a possible anxiogenic gene Crhr1 [47] was up-regulated whereas the anxiolytic gene Crhr2 [45, 46] was down-regulated (Table 1), implying the anxiogenic state at 2 weeks PS.

To date, genome-wide transcriptomic studies for PTSD-like model animals have been performed by several groups [68,69,70]. Muhie et al. reported using aggressor exposed animal models that genes involved in synaptic plasticity, neurogenesis, inflammation, obesity and cardiac infarction pathways were differentially regulated [71]. Even their model construction and sampling procedures were different, enriched KEGG pathways were remarkably similar to our results, especially for synaptic plasticity and neuroendocrine systems. On the other hand, human genomic analyses have identified genes with single nucleotide mutations associated with PTSD. Thus, the results generated from our PTSD-like mouse model constructed by a combination of inescapable foot shock (physical stressor) and repeated situational reminders (psychological stressor) is consistent with the previous reports from other laboratories. Nonetheless, we are aware that similar to the other animal models of human diseases, our current model system may only represent certain aspects of PTSD symptoms. We have not completely characterized our model animals to fulfill face validity, for example the Sucrose preference test to assess anhedonia, Social interaction test for social numbness, and Morris water maze test for memory impairment. Despite the limitation of the current model system, the RNA-seq analysis can partially reflect the PTSD-specific gene regulation, which is likely to provide a clue for molecular basis of temporal alteration in the brain regions crucial for emotional control and fear memory. This study may also help to identify target genes for intervention at particular brain regions and specific stages during the development of PTSD.

Change history

  • 09 May 2019

    Following publication of the original article [1], the authors reported that one of the authors’ names was spelled incorrectly.

Abbreviations

ACC:

Anterior cingulate cortex

AMY:

Amygdala

CNS:

Central nervous system

DAVID:

Database for Annotation, Visualization and Integrated Discovery

DEG:

Differentially expressed genes

GO:

Gene ontology

GR:

Glucocorticoid receptor

HPA:

Hypothalamic-pituitary-adrenal

PFC:

Prefrontal cortex

PS:

Post stress

PTSD:

Post-traumatic stress disorder

References

  1. 1.

    Zelikowsky M, Hersman S, Chawla MK, Barnes CA, Fanselow MS. Neuronal ensembles in amygdala, hippocampus, and prefrontal cortex track differential components of contextual fear. J Neurosci. 2014;34:8462–6.

  2. 2.

    Pitman RK, Rasmusson AM, Koenen KC, Shin LM, Orr SP, Gilbertson MW, Milad MR, Liberzon I. Biological studies of post-traumatic stress disorder. Nat Rev Neurosci. 2012;13:769–87.

  3. 3.

    Belzung C, Turiault M, Griebel G. Optogenetics to study the circuits of fear- and depression-like behaviors: a critical analysis. Pharmacol Biochem Behav. 2014;122:144–57.

  4. 4.

    Shin LM, Rauch SL, Pitman RK. Amygdala, medial prefrontal cortex, and hippocampal function in PTSD. Ann N Y Acad Sci. 2006;1071:67–79.

  5. 5.

    Etkin A, Wager TD. Functional neuroimaging of anxiety: a meta-analysis of emotional processing in PTSD, social anxiety disorder, and specific phobia. Am J Psychiatry. 2007;164:1476–88.

  6. 6.

    Rougemont-Bucking A, Linnman C, Zeffiro TA, Zeidan MA, Lebron-Milad K, Rodriguez-Romaguera J, Rauch SL, Pitman RK, Milad MR. Altered processing of contextual information during fear extinction in PTSD: an fMRI study. CNS Neurosci Ther. 2011;17:227–36.

  7. 7.

    Etkin A, Egner T, Peraza DM, Kandel ER, Hirsch J. Resolving emotional conflict: a role for the rostral anterior cingulate cortex in modulating activity in the amygdala. Neuron. 2006;51:871–82.

  8. 8.

    Shin LM, Whalen PJ, Pitman RK, Bush G, Macklin ML, Lasko NB, Orr SP, McInerney SC, Rauch SL. An fMRI study of anterior cingulate function in posttraumatic stress disorder. Biol Psychiatry. 2001;50:932–42.

  9. 9.

    Milad MR, Pitman RK, Ellis CB, Gold AL, Shin LM, Lasko NB, Zeidan MA, Handwerger K, Orr SP, Rauch SL. Neurobiological basis of failure to recall extinction memory in posttraumatic stress disorder. Biol Psychiatry. 2009;66:1075–82.

  10. 10.

    O'Doherty DC, Chitty KM, Saddiqui S, Bennett MR, Lagopoulos J. A systematic review and meta-analysis of magnetic resonance imaging measurement of structural volumes in posttraumatic stress disorder. Psychiatry Res. 2015;232:1–33.

  11. 11.

    Daskalakis NP, Rijal CM, King C, Huckins LM, Ressler KJ. Recent genetics and epigenetics approaches to PTSD. Curr Psychiatry Rep. 2018;20:30.

  12. 12.

    Girgenti MJ, Duman RS. Transcriptome alterations in posttraumatic stress disorder. Biol Psychiatry. 2017;83:840–8.

  13. 13.

    Nees F, Witt SH, Flor H. Neurogenetic approaches to stress and fear in humans as pathophysiological mechanisms for posttraumatic stress disorder. Biol Psychiatry. 2018;83:810–20.

  14. 14.

    Mahan AL, Ressler KJ. Fear conditioning, synaptic plasticity and the amygdala: implications for posttraumatic stress disorder. Trends Neurosci. 2012;35:24–35.

  15. 15.

    Schafe GE, Nader K, Blair HT, LeDoux JE. Memory consolidation of Pavlovian fear conditioning: a cellular and molecular perspective. Trends Neurosci. 2001;24:540–6.

  16. 16.

    Maguschak KA, Ressler KJ. A role for WNT/beta-catenin signaling in the neural mechanisms of behavior. J Neuroimmune Pharmacol. 2012;7:763–73.

  17. 17.

    Chen YH, Lan YJ, Zhang SR, Li WP, Luo ZY, Lin S, Zhuang JP, Li XW, Li SJ, Yang JM, et al. ErbB4 signaling in the prelimbic cortex regulates fear expression. Transl Psychiatry. 2017;7:e1168.

  18. 18.

    Bi LL, Sun XD, Zhang J, Lu YS, Chen YH, Wang J, Geng F, Liu F, Zhang M, Liu JH, et al. Amygdala NRG1-ErbB4 is critical for the modulation of anxiety-like behaviors. Neuropsychopharmacology. 2015;40:974–86.

  19. 19.

    Patel S, Roelke CT, Rademacher DJ, Hillard CJ. Inhibition of restraint stress-induced neural and behavioural activation by endogenous cannabinoid signalling. Eur J Neurosci. 2005;21:1057–69.

  20. 20.

    Yehuda R, Cai G, Golier JA, Sarapas C, Galea S, Ising M, Rein T, Schmeidler J, Muller-Myhsok B, Holsboer F, et al. Gene expression patterns associated with posttraumatic stress disorder following exposure to the world trade center attacks. Biol Psychiatry. 2009;66:708–11.

  21. 21.

    Yehuda R, Southwick SM, Krystal JH, Bremner D, Charney DS, Mason JW. Enhanced suppression of cortisol following dexamethasone administration in posttraumatic stress disorder. Am J Psychiatry. 1993;150:83–6.

  22. 22.

    Gillespie CF, Phifer J, Bradley B, Ressler KJ. Risk and resilience: genetic and environmental influences on development of the stress response. Depress Anxiety. 2009;26:984–92.

  23. 23.

    De Leo S, Lee SY, Braverman LE. Hyperthyroidism. Lancet. 2016;388:906–18.

  24. 24.

    Neumann ID, Slattery DA. Oxytocin in general anxiety and social fear: a translational approach. Biol Psychiatry. 2016;79:213–21.

  25. 25.

    Korem N, Akirav I. Cannabinoids prevent the effects of a footshock followed by situational reminders on emotional processing. Neuropsychopharmacology. 2014;39:2709–22.

  26. 26.

    Pynoos RS, Ritzmann RF, Steinberg AM, Goenjian A, Prisecaru I. A behavioral animal model of posttraumatic stress disorder featuring repeated exposure to situational reminders. Biol Psychiatry. 1996;39:129–34.

  27. 27.

    Tang X, Orchard SM, Sanford LD. Home cage activity and behavioral performance in inbred and hybrid mice. Behav Brain Res. 2002;136:555–69.

  28. 28.

    Fisher SP, Godinho SI, Pothecary CA, Hankins MW, Foster RG, Peirson SN. Rapid assessment of sleep-wake behavior in mice. J Biol Rhythm. 2012;27:48–58.

  29. 29.

    Pack AI, Galante RJ, Maislin G, Cater J, Metaxas D, Lu S, Zhang L, Von Smith R, Kay T, Lian J, et al. Novel method for high-throughput phenotyping of sleep in mice. Physiol Genomics. 2007;28:232–8.

  30. 30.

    Paxinos G, Franklin KBJ. The mouse brain in steriotaxic coordinates. San Diego: Academic Press; 2001.

  31. 31.

    Palkovits M. Isolated removal of hypothalamic or other brain nuclei of the rat. Brain Res. 1973;59:449–50.

  32. 32.

    Huang d W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44–57.

  33. 33.

    Yu G, Wang LG, Yan GR, He QY. DOSE: an R/Bioconductor package for disease ontology semantic and enrichment analysis. Bioinformatics. 2015;31:608–9.

  34. 34.

    Germain A. Sleep disturbances as the hallmark of PTSD: where are we now? Am J Psychiatry. 2013;170:372–82.

  35. 35.

    Ashley-Koch AE, Garrett ME, Gibson J, Liu Y, Dennis MF, Kimbrel NA, Veterans Affairs Mid-Atlantic Mental Illness Research E, Clinical Center W, Beckham JC, Hauser MA. Genome-wide association study of posttraumatic stress disorder in a cohort of Iraq-Afghanistan era veterans. J Affect Disord. 2015;184:225–34.

  36. 36.

    Duncan LE, Ratanatharathorn A, Aiello AE, Almli LM, Amstadter AB, Ashley-Koch AE, Baker DG, Beckham JC, Bierut LJ, Bisson J, et al. Largest GWAS of PTSD (N=20 070) yields genetic overlap with schizophrenia and sex differences in heritability. Mol Psychiatry. 2018;23:666–73.

  37. 37.

    Flaquer A, Baumbach C, Ladwig KH, Kriebel J, Waldenberger M, Grallert H, Baumert J, Meitinger T, Kruse J, Peters A, et al. Mitochondrial genetic variants identified to be associated with posttraumatic stress disorder. Transl Psychiatry. 2015;5:e524.

  38. 38.

    Wolf EJ, Rasmusson AM, Mitchell KS, Logue MW, Baldwin CT, Miller MW. A genome-wide association study of clinical symptoms of dissociation in a trauma-exposed sample. Depress Anxiety. 2014;31:352–60.

  39. 39.

    Borghans B, Homberg JR. Animal models for posttraumatic stress disorder: an overview of what is used in research. World J Psychiatry. 2015;5:387–96.

  40. 40.

    Schoner J, Heinz A, Endres M, Gertz K, Kronenberg G. Post-traumatic stress disorder and beyond: an overview of rodent stress models. J Cell Mol Med. 2017;21:2248–56.

  41. 41.

    Segev A, Korem N, Mizrachi Zer-Aviv T, Abush H, Lange R, Sauber G, Hillard CJ, Akirav I. Role of endocannabinoids in the hippocampus and amygdala in emotional memory and plasticity. Neuropsychopharmacology. 2018;43:2017–27.

  42. 42.

    Goswami S, Rodriguez-Sierra O, Cascardi M, Pare D. Animal models of post-traumatic stress disorder: face validity. Front Neurosci. 2013;7:89.

  43. 43.

    Zovkic IB, Meadows JP, Kaas GA, Sweatt JD. Interindividual variability in stress susceptibility: a role for epigenetic mechanisms in PTSD. Front Psychiatry. 2013;4:60.

  44. 44.

    Matar MA, Zohar J, Cohen H. Translationally relevant modeling of PTSD in rodents. Cell Tissue Res. 2013;354:127–39.

  45. 45.

    Bale TL, Contarino A, Smith GW, Chan R, Gold LH, Sawchenko PE, Koob GF, Vale WW, Lee KF. Mice deficient for corticotropin-releasing hormone receptor-2 display anxiety-like behaviour and are hypersensitive to stress. Nat Genet. 2000;24:410–4.

  46. 46.

    Kishimoto T, Radulovic J, Radulovic M, Lin CR, Schrick C, Hooshmand F, Hermanson O, Rosenfeld MG, Spiess J. Deletion of crhr2 reveals an anxiolytic role for corticotropin-releasing hormone receptor-2. Nat Genet. 2000;24:415–9.

  47. 47.

    Timpl P, Spanagel R, Sillaber I, Kresse A, Reul JM, Stalla GK, Blanquet V, Steckler T, Holsboer F, Wurst W. Impaired stress response and reduced anxiety in mice lacking a functional corticotropin-releasing hormone receptor 1. Nat Genet. 1998;19:162–6.

  48. 48.

    Frederic F, Chianale C, Oliver C, Mariani J. Enhanced endocrine response to novel environment stress and lack of corticosterone circadian rhythm in staggerer (Rora sg/sg) mutant mice. J Neurosci Res. 2006;83:1525–32.

  49. 49.

    Luca F, Kashyap S, Southard C, Zou M, Witonsky D, Di Rienzo A, Conzen SD. Adaptive variation regulates the expression of the human SGK1 gene in response to stress. PLoS Genet. 2009;5:e1000489.

  50. 50.

    Sarapas C, Cai G, Bierer LM, Golier JA, Galea S, Ising M, Rein T, Schmeidler J, Muller-Myhsok B, Uhr M, et al. Genetic markers for PTSD risk and resilience among survivors of the world trade center attacks. Dis Markers. 2011;30:101–10.

  51. 51.

    Mehta D, Gonik M, Klengel T, Rex-Haffner M, Menke A, Rubel J, Mercer KB, Putz B, Bradley B, Holsboer F, et al. Using polymorphisms in FKBP5 to define biologically distinct subtypes of posttraumatic stress disorder: evidence from endocrine and gene expression studies. Arch Gen Psychiatry. 2011;68:901–10.

  52. 52.

    Schaefer ML, Wong ST, Wozniak DF, Muglia LM, Liauw JA, Zhuo M, Nardi A, Hartman RE, Vogt SK, Luedke CE, et al. Altered stress-induced anxiety in adenylyl cyclase type VIII-deficient mice. J Neurosci. 2000;20:4809–20.

  53. 53.

    Aromolaran KA, Benzow KA, Cribbs LL, Koob MD, Piedras-Renteria ES. T-type current modulation by the actin-binding protein Kelch-like 1. Am J Physiol Cell Physiol. 2010;298:C1353–62.

  54. 54.

    Seng S, Avraham HK, Jiang S, Venkatesh S, Avraham S. KLHL1/MRP2 mediates neurite outgrowth in a glycogen synthase kinase 3beta-dependent manner. Mol Cell Biol. 2006;26:8371–84.

  55. 55.

    Mao W, Watanabe T, Cho S, Frost JL, Truong T, Zhao X, Futai K. Shank1 regulates excitatory synaptic transmission in mouse hippocampal parvalbumin-expressing inhibitory interneurons. Eur J Neurosci. 2015;41:1025–35.

  56. 56.

    Kuan PF, Waszczuk MA, Kotov R, Marsit CJ, Guffanti G, Gonzalez A, Yang X, Koenen K, Bromet E, Luft BJ. An epigenome-wide DNA methylation study of PTSD and depression in world trade center responders. Transl Psychiatry. 2017;7:e1158.

  57. 57.

    Rutten BPF, Vermetten E, Vinkers CH, Ursini G, Daskalakis NP, Pishva E, de Nijs L, Houtepen LC, Eijssen L, Jaffe AE, et al. Longitudinal analyses of the DNA methylome in deployed military servicemen identify susceptibility loci for post-traumatic stress disorder. Mol Psychiatry. 2018;23:1145–56.

  58. 58.

    Hill MN, Campolongo P, Yehuda R, Patel S. Integrating endocannabinoid signaling and cannabinoids into the biology and treatment of posttraumatic stress disorder. Neuropsychopharmacology. 2018;43:80–102.

  59. 59.

    Dhaka A, Costa RM, Hu H, Irvin DK, Patel A, Kornblum HI, Silva AJ, O'Dell TJ, Colicelli J. The RAS effector RIN1 modulates the formation of aversive memories. J Neurosci. 2003;23:748–57.

  60. 60.

    Pan BX, Vautier F, Ito W, Bolshakov VY, Morozov A. Enhanced cortico-amygdala efficacy and suppressed fear in absence of Rap1. J Neurosci. 2008;28:2089–98.

  61. 61.

    Sabban EL, Serova LI, Newman E, Aisenberg N, Akirav I. Changes in gene expression in the locus Coeruleus-amygdala circuitry in inhibitory avoidance PTSD model. Cell Mol Neurobiol. 2018;38:273–80.

  62. 62.

    Louvart H, Maccari S, Lesage J, Leonhardt M, Dickes-Coopman A, Darnaudery M. Effects of a single footshock followed by situational reminders on HPA axis and behaviour in the aversive context in male and female rats. Psychoneuroendocrinology. 2006;31:92–9.

  63. 63.

    Campbell-Smith EJ, Holmes NM, Lingawi NW, Panayi MC, Westbrook RF. Oxytocin signaling in basolateral and central amygdala nuclei differentially regulates the acquisition, expression, and extinction of context-conditioned fear in rats. Learn Mem. 2015;22:247–57.

  64. 64.

    Viviani D, Charlet A, van den Burg E, Robinet C, Hurni N, Abatis M, Magara F, Stoop R. Oxytocin selectively gates fear responses through distinct outputs from the central amygdala. Science. 2011;333:104–7.

  65. 65.

    Winter J, Jurek B. The interplay between oxytocin and the CRF system: regulation of the stress response. Cell Tissue Res. 2019;375:85–91.

  66. 66.

    Shin LM, Lasko NB, Macklin ML, Karpf RD, Milad MR, Orr SP, Goetz JM, Fischman AJ, Rauch SL, Pitman RK. Resting metabolic activity in the cingulate cortex and vulnerability to posttraumatic stress disorder. Arch Gen Psychiatry. 2009;66:1099–107.

  67. 67.

    Hamner MB, Lorberbaum JP, George MS. Potential role of the anterior cingulate cortex in PTSD: review and hypothesis. Depress Anxiety. 1999;9:1–14.

  68. 68.

    Burgdorf J, Kroes RA, Zhang XL, Gross AL, Schmidt M, Weiss C, Disterhoft JF, Burch RM, Stanton PK, Moskal JR. Rapastinel (GLYX-13) has therapeutic potential for the treatment of post-traumatic stress disorder: characterization of a NMDA receptor-mediated metaplasticity process in the medial prefrontal cortex of rats. Behav Brain Res. 2015;294:177–85.

  69. 69.

    Ponomarev I, Rau V, Eger EI, Harris RA, Fanselow MS. Amygdala transcriptome and cellular mechanisms underlying stress-enhanced fear learning in a rat model of posttraumatic stress disorder. Neuropsychopharmacology. 2010;35:1402–11.

  70. 70.

    Sillivan SE, Joseph NF, Jamieson S, King ML, Chevere-Torres I, Fuentes I, Shumyatsky GP, Brantley AF, Rumbaugh G, Miller CA. Susceptibility and resilience to posttraumatic stress disorder-like behaviors in inbred mice. Biol Psychiatry. 2017;82:924–33.

  71. 71.

    Muhie S, Gautam A, Meyerhoff J, Chakraborty N, Hammamieh R, Jett M. Brain transcriptome profiles in mouse model simulating features of post-traumatic stress disorder. Mol Brain. 2015;8:14.

  72. 72.

    Krzyzewska IM, Ensink JBM, Nawijn L, Mul AN, Koch SB, Venema A, Shankar V, Frijling JL, Veltman DJ, Lindauer RJL, et al. Genetic variant in CACNA1C is associated with PTSD in traumatized police officers. Eur J Hum Genet. 2018;26:247–57.

  73. 73.

    Castro-Vale I, van Rossum EF, Machado JC, Mota-Cardoso R, Carvalho D. Genetics of glucocorticoid regulation and posttraumatic stress disorder--what do we know? Neurosci Biobehav Rev. 2016;63:143–57.

  74. 74.

    Wolf EJ, Mitchell KS, Logue MW, Baldwin CT, Reardon AF, Humphries DE, Miller MW. Corticotropin releasing hormone receptor 2 (CRHR-2) gene is associated with decreased risk and severity of posttraumatic stress disorder in women. Depress Anxiety. 2013;30:1161–9.

  75. 75.

    Holmes SE, Girgenti MJ, Davis MT, Pietrzak RH, DellaGioia N, Nabulsi N, Matuskey D, Southwick S, Duman RS, Carson RE, et al. Altered metabotropic glutamate receptor 5 markers in PTSD: in vivo and postmortem evidence. Proc Natl Acad Sci U S A. 2017;114:8390–5.

  76. 76.

    Young KA, Thompson PM, Cruz DA, Williamson DE, Selemon LD. BA11 FKBP5 expression levels correlate with dendritic spine density in postmortem PTSD and controls. Neurobiol Stress. 2015;2:67–72.

  77. 77.

    Sippel LM, Han S, Watkins LE, Harpaz-Rotem I, Southwick SM, Krystal JH, Olff M, Sherva R, Farrer LA, Kranzler HR, et al. Oxytocin receptor gene polymorphisms, attachment, and PTSD: results from the National Health and resilience in veterans study. J Psychiatr Res. 2017;94:139–47.

  78. 78.

    Logue MW, Baldwin C, Guffanti G, Melista E, Wolf EJ, Reardon AF, Uddin M, Wildman D, Galea S, Koenen KC, et al. A genome-wide association study of post-traumatic stress disorder identifies the retinoid-related orphan receptor alpha (RORA) gene as a significant risk locus. Mol Psychiatry. 2013;18:937–42.

  79. 79.

    Licznerski P, Duric V, Banasr M, Alavian KN, Ota KT, Kang HJ, Jonas EA, Ursano R, Krystal JH, Duman RS, et al. Decreased SGK1 expression and function contributes to behavioral deficits induced by traumatic stress. PLoS Biol. 2015;13:e1002282.

Download references

Acknowledgements

We are thankful to Ms. Laura Tucker for maintenance of animal testing room, and to Dr. David Jacobowitz for assistance to brain dissection.

Funding

This research was supported by NHLBI/USU Collaborative Health Initiative Research Program (CHIRP) grant (308431-9.00-64532).

Availability of data and materials

The datasets during and/or analysed during the current study available from the corresponding author on reasonable request.

Author information

MT and YZ planned and designed the studies. MT performed behavioral tests, and HL prepared samples, and XZ, JS, MW, and CD ran RNA sequencing and analyzed data. MT and XZ performed functional annotation analyses. MT, XZ and YZ drew figures and tables. MT and YZ wrote manuscript. All authors read and approved the final manuscript.

Correspondence to Yumin Zhang.

Ethics declarations

Ethics approval

All the experimental procedures involving animals were conducted in accordance with the NIH guidelines and approved ethically by the Uniformed Services University Animal Care and Use Committee.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional information

The original version of this article was revised: the name of one of the authors had been presented incorrectly. It should be Clifton L. Dalgard, not Clifton Dalgard.

Additional files

Additional file 1:

Figure S1. RNAseq analysis in Anterior Cingulate Cortex (ACC) of stressed mice. Total RNA was isolated from ACC by punching brain slice and applied for RNA sequence. Heatmaps of DEGs were shown for 2 weeks (A) and 5 weeks (B) PS. Green bars and orange bars indicated control and stressed animals, respectively. Venn diagram of overlapped DEGs between the two time points was shown in (C). (PPTX 496 kb)

Additional file 2:

Figure S2. RNAseq analysis in Amygdala (AMY) of stressed mice. Total RNA was isolated from AMY by punching brain slice and applied for RNA sequence. Heatmaps of DEGs were shown for 2 weeks (A) and 5 weeks (B) PS. Green bars and orange bars indicated control and stressed animals, respectively. Venn diagram between the two conditions was shown in (C). (PPTX 491 kb)

Additional file 3:

Table S1. Enriched molecular functions of regulated genes in AMY and ACC at 2 and 5 weeks post stress. (PPTX 44 kb)

Additional file 4:

Table S2. Top 20 enriched biological processes of regulated genes in AMY and ACC at 2 and 5 weeks post stress. (PPTX 47 kb)

Additional file 5:

Figure S3. Heatmap of enriched GOs related to transcriptional regulation. Color index represents level of significance (p-values). (PPTX 54 kb)

Additional file 6:

Figure S4. Heatmap of enriched GOs involved in neuronal development, neurogenesis, and differentiation, and non-neuronal cell differentiation. Color index represents level of significance (p-values). (PPTX 53 kb)

Additional file 7:

Table S3. Enriched pathways of regulated genes in AMY and ACC at 2 and 5 weeks post stress. (PPTX 45 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Keywords

  • Synaptic remodeling
  • Neuroendocrine
  • Long-term regulation
  • Anxiety
  • Avoidance
  • Hyperarousal
  • Sleep
  • RNA-seq
  • Amygdala
  • Anterior cingulate cortex