- 1State Key Laboratory of Chemical Oncogenomics, Guangdong Provincial Key Laboratory of Chemical Genomics, Peking University Shenzhen Graduate School, Shenzhen, China
- 2Department of Addiction Medicine, Shenzhen Engineering Research Center for Precision Psychiatric Technology, Shenzhen Clinical Research Center for Mental Disorders, Shenzhen Kangning Hospital, Shenzhen Mental Health Center, Shenzhen, China
- 3Henan Key Laboratory of Medical Tissue Regeneration, Xinxiang Medical University, Xinxiang, China
- 4National Institute on Drug Dependence and Beijing Key Laboratory of Drug Dependence, Peking University, Beijing, China
- 5Affiliated Mental Health Center, Southern University of Science and Technology, Shenzhen, China
- 6Clinical College of Mental Health, ShenZhen University Health Science Center, Shenzhen, China
- 7School of Mental Health, Jining Medical University, Jining, China
- 8School of Mental Health, Anhui Medical University, Hefei, China
The experience of traumatic stress can engender lasting memories associated with the trauma, often resulting in post-traumatic stress disorder (PTSD). However, only a minority of individuals develop PTSD symptoms upon exposure. The neurobiological mechanisms underlying the pathology of PTSD are poorly understood. Utilizing a rat model of PTSD, the Single Prolonged Stress (SPS) paradigm, we were able to differentiate between resilient and susceptible individuals. Fourteen days after the SPS exposure, we conducted the behavioral analyses using Elevated Plus Maze (EPM) and Open Field (OF) tests to identify male rats as trauma resilient or susceptible. We focused on the microRNA (miRNA) profiles of the infralimbic (IL) and prelimbic (PL) cortical regions, known to be crucial in regulating the stress response. Our investigation of stressed rats exposed to the SPS procedure yielded divergent response, and differential expression microRNAs (DEmiRs) analysis indicated significant differences in the IL and PL transcriptional response. In the IL cortex, the GO analysis revealed enriched GO terms in the resilient versus control comparison, specifically related to mitogen-activated protein kinase and MAP kinase signaling pathways for their molecular functions as well as cytosol and nucleoplasm for the biological process. In the susceptible versus resilient comparison, the changes in molecular functions were only manifested in the functions of regulation of transcription involved in the G1/S transition of the mitotic cell cycle and skeletal muscle satellite cell activation. However, no enriched GO terms were found in the susceptible versus control comparison. In the PL cortex, results indicated that the DEmiRs were enriched exclusively in the cellular component level of the endoplasmic reticulum lumen in the comparison between resilient and control rats. Overall, our study utilized an animal model of PTSD to investigate the potential correlation between stress-induced behavioral dysfunction and variations in miRNA expression. The aforementioned discoveries have the potential to pave the way for novel therapeutic approaches for PTSD, which could involve the targeted regulation of transcriptome expression.
Introduction
The occurrence of a life-threatening event can give rise to the development of PTSD which is distinguished by intrusive memories, avoidance of stimuli that trigger recollections of the event, recurring nightmares, and persistent hyperarousal (1). The symptoms of PTSD may not become apparent immediately after the trauma and may exhibit a delayed onset. Some symptoms may show partial remission and relapse, while others may persist and worsen over time. Despite most people undergo at least one traumatic experience in their lives, only a small fraction, constituting 10–20% of the population, develop enduring manifestations of PTSD (2). Conversely, some other individuals display resilience and do not develop lasting psychological consequences (3). Resilience, which refers to the ability to adapt effectively in the face of significant adversity, is often the prevailing response pattern observed in individuals who have experienced trauma (4, 5). Understanding the neural mechanisms that differentiate individuals from those vulnerable to stress-related disorders is crucial in the prevention of psychiatric disorders and the development of innovative treatment approaches.
Studies conducted on both humans and animals have suggested that the abnormalities in neurocircuits that are linked to emotions and fear, specifically the amygdala, prefrontal cortex (PFC), and hippocampus, may contribute to the development and maintenance of symptoms of PTSD (6). Moreover, aberrations in the medial prefrontal cortex (mPFC) function have been implicated in the perpetuation of fearfulness in anxiety disorders, such as PTSD (7–9). Advanced high-precision neuroimaging and optogenetic methods have confirmed that PFC is an essential brain structure in the regulation of PTSD (10). Studies of rodents have shown that SPS leads to decreased neuronal activation in the IL, which could be a factor in the associated fear extinction deficits (11–13). In contrast, optogenetic excitation of the mPFC has been linked to increased stress resilience, and successful stress coping is associated with heightened mPFC activation following social defeat stress (14). Overall, it becomes evident that the PFC plays a pivotal role in regulating an individual’s susceptibility or resilience to stress. In light of the aforementioned considerations, more analysis to the PFC subregions, including the IL and PL cortex (15), should be further investigated.
The factors that contribute to the development of PTSD in some individuals following exposure to traumatic events remain largely unknown. Nevertheless, it is likely that alterations in the molecular mechanisms responsible for regulating stress vulnerability and adaptability are implicated in this process (16). Epigenetic modifications are primarily regulated through DNA methylation, histone modifications and the control of non-coding RNAs, and any disruption of these modifications can be detrimental to normal neuronal development and functioning, resulting in severe neuropsychiatric disorders (17, 18). Specifically, MicroRNAs (miRNAs), a class of small non-coding RNAs, have emerged as potential biomarkers for various diseases, including psychiatric disorders (19, 20). Acting as essential regulators of gene expression, miRNA primarily repress the expression of their target mRNAs (21, 22). Numerous studies have suggested that miRNAs are involved in the molecular mechanisms of fear memory (23), and stress reactions associated with PTSD (24, 25). In animal studies, altered miRNA expression in the PFC (26) and amygdala (27) of rodents exposed to traumatic stress, exhibited PTSD-like symptoms. Moreover, the circulating miRNAs were examined as potential diagnostic biomarkers in individuals with PTSD and corresponding controls (28, 29). These results indicate that the aberrant miRNA expression observed in PTSD may be correlated with immune responses and neuronal activation. The examination of alterations in gene expression in the key regions of PTSD, particular subregions of the PFC, will provide insights into the fundamental neurobiological processes, ultimately leading to individualized treatment that prevent susceptibility or promote resilience to stress. In doing so, we utilized a transcriptome analysis of differential miRNA expression with subsequent gene ontology-based pathway analysis. Our primary objective is to investigate the transcriptome associated with susceptibility and resilience phenotypes in IL and PL cortex, aiming to comprehend the underlying molecular mechanisms and identify innovative therapeutic targets and biomarkers.
Methods
Animals
Male Sprague–Dawley (SD) rats, approximately 3 months old, were procured from Beijing Vital River Laboratory Animal Technology Co., Ltd. These animals were housed in a controlled environment with a 12 h light/dark cycle and provided with ad libitum access to food and water. They were housed in groups of four per cage for 1 week prior to the commencement of the experiment. The light/dark cycle was maintained at 12 h each. Subsequently, some of the rats were randomly subjected to the SPS protocol to induce PTSD-like behavior, while the remaining rats were kept as controls in separate cages. The behavioral experiments were conducted during the dark phase of the cycle. All animal-related experiments were conducted in accordance with the National Institutes of Health Guide for the Care and Use of Laboratory Animals and the Biomedical Ethics Committee of Peking University’s regulations on animal use and protection. The Biomedical Ethics Committee of Peking University also approved the protocol for animal use and protection.
Experimental design
The experimental design is shown in Figure 1A. Thirty rats (260–280 g) were acclimatized for 7 days, then twenty underwent the SPS stressors, and the other ten rats as control. After undisturbed for 14 days, they were tested on OF and EPM. The animals were sacrificed after the behavioral test, the PL and IL cortex were dissected and stored at −80°C until used for RNA sequencing analysis.
Figure 1. Experimental design and behavioral test. (A) Overview of experimental design. (B) Effect of SPS on anxiety index in control and SPS group. (C) Effect of SPS on anxiety index in control, susceptible, and resilient group. (D) Time spent in open arms. (E) Number of open arm entries. (F) Time spent in central area. (G) Distance spent in central area.
Single prolonged stress
An animal model for PTSD-like symptoms was constructed by the SPS procedure as previously described (30). The SPS rat model is widely regarded as one of the most reliable models for studying PTSD. The procedure included a 2 h restraint stress in a plastic restrainer, followed by a 20 min forced swim at a depth of 40 cm and a temperature of 25°C. After a 15 min rest, the rats were anesthetized using diethyl to induce unconsciousness.
Behavioral test
The OF and EPM test was based on our previous studies (30) and was conducted in a procedure room under dim light conditions. To eradicate the previous animal’s odor, the apparatus underwent a cleaning procedure utilizing a towel drenched in a 75% ethanol solution. These two tests were recorded for 5 min per rat and then analyzed with LabMaze v3.0 (Zhongshi Technology, China). In the OF test, the total travel distance and the time spent in the central area were counted to measure the rats’ behavior. In the EPM test, the number of open/closed arm entries and the time spent in the open/closed arms were recorded. Vulnerable and resilient subpopulations of rats emerge after SPS. According to the previous study (31), individual animals were retrospectively classified as having behavioral responses according to the anxiety index (AI), which is calculated as 1 − [(time spent in open arm/total time on the maze)/2 + (number of entries to the open arms/total exploration on the maze)/2]. For classification criteria, one standard deviation below the mean of the control group’s anxiety index as resilient individuals, and one standard deviation above the mean of the the control group’s anxiety index as susceptible individuals. After SPS procedure, 6 rats reached susceptible criteria and 5 rats reached resilient criteria. A subset of 12 animals were selected for RNA sequencing based on the anxiety index: 4 out of 5 resilient individuals, 4 of 6 susceptible individuals, and 4 out of 10 control animals were randomly chosen.
RNA sequencing and analysis
The animals were sacrificed after the EPM and OF tests, and samples of IL and PL cortex were dissected from each animal. Then, the tissue samples were sent to BGI Co., LTD (Shenzhen, China), and the BGISEQ-500 platform was applied to perform RNA sequencing. DEmiRs were determined based on |log2(fold change) | ≥ 0.5 and q value (Adjusted p-value) < = 0.05. To gain insight into the biological roles of the differentially expressed miRNAs, we used the gene ontology (GO)-term and the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses. All analysis of RNA sequencing data was performed by DR. TOM online system provided by BGI.1
Statistical analysis
The behavioral data were expressed as mean ± SEM. The data were analyzed using an analysis of t-test or one-way ANOVA followed by Bonferroni’s multiple comparison test (see Results). Fisher’s exact test was used to compare the distribution of each GO classification or KEGG pathway in the target miRNA set, and the q value (Adjusted p-value) < = 0.05. Statistical significance was set at p < 0.05. *, **, and *** depicting p < 0.05, p < 0.01, and p < 0.001, respectively.
Results
Divergent responses on EPM and OF after SPS
Fourteen days after the SPS procedure, animals were tested on OF and EPM. Results analysis showed that the anxiety index in the SPS group was significantly higher compared to the control group (t = 3.10, p < 0.001, Figure 1B). To better understand the susceptibility and resilience to stress, SPS animals were classified into two groups, namely susceptible and resilient, based on their performance in the anxiety index assessment. These groups differed significantly on anxiety index (F2, 18 = 54.25, p < 0.001, Figure 1C), specifically, the susceptible group exhibited a statistically significant difference in anxiety index as compared to both the resilient and control groups (p < 0.001). However, no significant difference was observed between the resilient and control groups. In the EPM test, these groups differed significantly on time (F2, 18 = 27.19, p < 0.001, Figure 1D) and the number of entries (F2, 18 = 18.84, p < 0.001, Figure 1E) into the open arms, but the significant difference only appear in susceptible group when compare to control and resilient group (p < 0.001). In the OF test, the susceptible group displayed more anxiety behaviors, and tended to spend less time (F2, 18 = 7.421, p < 0.001, Figure 1F) and trace on central area (F2, 18 = 6.414, p < 0.001, Figure 1G), in contrast, no significant changes were observed in resilient group and control group.
Regional-specific changes in miRNAs expression
We performed small RNA-seq on the IL and PL cortex and compared the miRNA expression between groups. Differential expression analysis of each subset reveals significant differences in both genes being up- and downregulated in IL and PL cortex. For DEmiRs in the IL cortex, we identified 9 DEmiRs in the resilient versus control comparison, 6 DEmiRs in the susceptible versus control comparison, and 2 DEmiRs in resilient versus susceptible comparison (Figure 2A). For DEmiRs in PL cortex, there are 6 DEmiRs in the resilient versus control comparison, 10 DEmiRs in the susceptible versus control comparison, and 8 DEmiRs in resilient versus susceptible comparison (Figure 2B). Venn diagram demonstrating shared DEmiRs were plotted, comparison of resilient/control and susceptible/control profiles showed only 2 common DEmiRs in IL cortex (Figure 2C), and no common DEmiRs in PL cortex. However, between resilient/control and resilient/susceptible, we only found 1 overlapping DEmiR, and we also found 2 overlapping DEmiRs between resilient/control and resilient/susceptible in PL cortex (Figure 2D). Only one gene overlapped between IL and PL cortex in susceptible versus control comparison (Figure 2E), and two genes overlapped between IL and PL cortex in resilient versus control comparison (Figure 2F).
Figure 2. Differentially expressed microRNA (DEmiRs). (A) DEmiRs within the IL cortex. (B) DEmiRs within the PL cortex. (C) Venn diagram of DEmiRs in IL cortex. (D) Venn diagram of DEmiRs in PL cortex. (E) Venn diagram of susceptible versus control comparison between IL and PL cortex. (F) Venn diagram of resilient versus control comparison between IL and PL cortex.
Gene ontology analysis of differentially expressed miRNAs
GO enrichment analysis was used for gene annotation and function description, containing three sub-ontologies of defined terms that describe gene product attributes, such as biological processes (BP), cellular components (CC), and molecular functions (MF). In IL cortex, the gene enrichment of miRNA and coexpressed mRNAs in susceptible versus control comparison did not significantly participate in any signaling pathways (Figures 3A–C). GO analysis of resilient versus control comparison showed that changes in the cellular and molecular function were enriched in the mitogen-activated protein kinase (q value = 0.04) and MAP kinase signaling pathways (q value = 0.04) (Figures 3D,E). And the changes in the biological process of DEmiRs were enriched in the cytosol (q value<0.01) and nucleoplasm (q value = 0.01) (Figure 3F). While in susceptible versus resilient comparison, the changes in biological functions except cellular and molecular were only manifested in the functions of regulation of transcription involved in G1/S transition of mitotic cell (q value = 0.02) and skeletal muscle satellite cell activation (q value = 0.02) (Figures 3G–I). In susceptible vs control, no important cellular, molecular and biological functions were found in PL cortex (Figures 4A–C). There were no important GO enrichments at molecular and biological levels (Figures 4E,F). Similarly, the GO enrichments were not significant at cellular, molecular and biological levels in susceptible vs resilient comparison (Figures 4G,H).
Figure 3. Gene Ontology enrichment analysis for DEmiR in IL cortex. (A–C) GO terms in susceptible versus control comparison. (D–F) GO terms in resilient versus control comparison. (G–I) GO terms in susceptible versus resilient comparison.
Figure 4. Gene Ontology enrichment analysis for DEmiR in PL cortex. (A–C) GO terms in susceptible versus control comparison. (D–F) GO terms in resilient versus control comparison. (G–I) GO terms in susceptible versus resilient comparison.
KEEG ontology analysis of differentially expressed miRNAs
Based on KEEG classification, KEGG analysis in the IL cortex mainly involved Axon regeneration in susceptible versus resilient comparison (q value = 0.028, Figure 5C). However, no altered pathways were obtained in susceptible versus control and resilient versus control comparisons (Figures 5A, B). As for PL cortex, KEGG analysis suggested that DEmiRs did not involve any pathway between group comparison (Figures 5D–F).
Figure 5. KEEG pathway analysis of DEmiRs in different comparisons. (A–C) KEGG pathway enrichment analysis of DEmiRs in IL cortex was performed in susceptible versus control, resilient versus control, and susceptible versus resilient comparisons. (D–F) KEGG pathway enrichment analysis of DEmiRs in PL cortex was performed in susceptible versus control, resilient versus control, and susceptible versus resilient comparisons.
Overall, our research employed an animal model of PTSD to explore the association between stress-induced behavioral issues and alterations in miRNA expression. Results analysis of rats exposed to the SPS procedure revealed disparate reactions, and DEmiRs analysis revealed considerable differences in the IL and PL transcriptional response.
Discussion
In this study, the SPS procedure was utilized to induce different stress response phenotypes, as determined by anxiety index. We profiled the miRNAs in the IL and PL cortex of resilient and susceptible rats and subsequently explored the biological alterations of the DEmiRs by enrichment analysis.
SPS-exposed rats exhibited more severe anxiety behaviors than control rats after a long duration. We also found divergent responses among the SPS-exposed rats, resulting in two distinct phenotypes: susceptible and resilient. Animal models of PTSD have proven to be highly effective in assessing the underlying biochemical mechanisms that contribute to the development of long-lasting impairments. These models also facilitate the identification of biomarkers and the testing of potential therapeutic interventions (32, 33). The behavioral impairments observed in SPS rats resemble the key symptoms associated with PTSD, indicating the reliability of the model (34–36). Moreover, previous studies also showed divergent responses to the SPS in rats (31, 37, 38), indicating the ability of SPS to induce behavioral alterations that reproduce PTSD characteristics only in a subset of exposed individuals. According to the anxiety index, approximately 30% of the SPS-exposed rats can be considered as susceptible, while other 25% of rats as resilient. The current findings further support the reliability of the SPS as an animal model for PTSD, as it mimics the heterogeneous response observed in individuals with PTSD. Individuals who are prone to stress often encounter difficulties in adjusting to stressors and may manifest maladaptive reactions over time. On the other hand, resilient individuals possess the ability to perceive adversity as relatively harmless and can establish adaptive physiological and psychological coping strategies (39). Each individual may perceive exposure to stressful events differently, and the resulting effects can be long-lasting, contingent upon the degree of stress resilience or vulnerability of each individual. Randomly assigning rats to different groups could reduce individual differences, and detecting anxiety levels prior to the SPS procedure would help to identify rats prone to stress. In animal research, individual behavioral profiling allows investigating neurobiological alterations related to susceptibility or resilience in animal models of PTSD and is utilized here to examine the activation of differential transcriptome expression of the stress system associated with trauma susceptibility or resilience.
Our findings demonstrated that the IL and PL cortex have different miRNA expression patterns, suggesting differential contributions within subregions to the stress response. Interestingly, transcriptional regulation of each subregion may have an influential, expression-independent role in PTSD-like behavior, even though only a few modified transcripts are differentially expressed. The mPFC is known to be involved in behavioral and executive functions, which are impacted by stress in both rodents and humans and are disrupted in individuals with stress-related disorders like PTSD and Major Depressive Disorder (40, 41). Human fMRI studies have revealed that increased prefrontal activity may serve as a resilience factor against PTSD, indicating that improved top-down control of subcortical regions promotes more effective emotion regulation (42). Thus, investigating the mPFC, particularly when its impact on vulnerability to stress, is crucial for the discovery of novel biomarkers that can distinguish between resilience and vulnerability to stress, offering valuable insights into the prevention and treatment of stress-related psychiatric disorders. In the present study, we proposed that the subregions of the mPFC, specifically the PL and IL cortex, would be impacted differently by traumatic stress after exposing rats to the SPS procedure. It has been observed that SPS causes a hypoactivity of the ventromedial prefrontal cortex (vmPFC) and hyperactivity of the amygdala, leading to a lack of top-down control, which is thought to be the source of the abnormal fear responses observed in SPS, and is similarly seen in PTSD patients (12, 43). Evidence also suggests that SPS-induced extinction impairments can be reversed with chemogenetic activation of the IL, demonstrating the crucial role of IL cortex in stress-related behavioral deficits (44). Also, the IL is a vital component in the processes of extinction, extinction recall, and reinstatement of conditioned fear (45), and damage or inactivation of this structure leads to extinction memory deficits that are similar to those observed in individuals with PTSD (46). Moreover, the IL region has direct and indirect projections to the intercalated neurons which suppress activation of the CeA, while the PL region sends signals to the BLA to arouse the CeA (47), which suggests differential contributions of different subregions of the mPFC. Notably, previous findings demonstrated enhancement of extinction of fear memory by chemogenetic activation of mPFC excitatory neurons only appeared in IL but not the PL region (44). Moreover, SPS was found to decrease neural activity in the IL but not PL, as observed in manganese-enhanced magnetic resonance imaging (12). Further investigation is needed to elucidate the roles of subregions in PTSD development. Moreover, gender is a fundamental factor to consider when examining PTSD, as women are two times more likely to suffer from it than men (48), consequently, more research is required in the future to explore the gender factor in PTSD.
The present study identified core PFC transcriptional signatures that are associated with behavioral phenotypes, including resilience and susceptibility. Our analysis of DEmiRs data suggests significant distinctions exist between the SPS-exposed groups in IL and PL cortex. Specially, we identified 17 DEmiRs in IL cortex and 24 DEmiRs in PL cortex. Interestingly, there was limited overlap between the IL and PL cortex comparisons, with only one gene (miR-205) overlapping between IL and PL cortex in susceptible versus control comparison, and two genes (miR-133c, miR-3,587) overlapping between IL and PL cortex in resilient versus control comparison. Previous studies indicated that miR-142, miR-124, and miR-153 regulation in the hippocampus may be important therapeutic targets for PTSD and could alleviate SPS-exposed PTSD-like behaviors (24, 25, 49). However, there is currently limited literature on PFC transcriptome changes after SPS. In the IL cortex, the GO analysis revealed enriched GO terms in the resilient versus control, including mitogen-activated protein kinase, MAP kinase signaling pathways for their molecular functions, as well as cytosol and nucleoplasm for the biological process. In the susceptible versus resilient comparison, the changes in molecular functions were primarily related to the functions of regulation of transcription involved in G1/S transition of mitotic cell cycle and skeletal muscle satellite cell activation. Notably, no enriched GO terms were found in susceptible versus control comparison. In the PL cortex, the GO enrichment analysis conducted on the target genes indicated that the DEmiRs were enriched exclusively in the cellular component level of endoplasmic reticulum lumen in the comparison between resilient and control rats. As far as we know, this study is the first to investigate the transcriptome expression profile in the mPFC after SPS in rats. A previous study examining the immediate changes to the transcriptional architecture in the PFC in response to footshock stress in rats showed no substantial changes at the single gene level immediately after the stress session. However, gene set enrichment analysis demonstrated modifications in neuronal pathways linked to glia maturation, glia–neuron communication, and synaptic activity (50). Other brain areas, including nucleus accumbens (NAc) and locus coeruleus (LC), have been investigated in rats using the SPS model, revealing the DEGs in the LC were found to be associated with morphological changes, including the regulation of actin cytoskeleton, regulation of cell size, growth factor activity, brain development and memory, and few DEGs appeared in NAc (31). However, this study only highlights different transcriptional profiles within the LC and NAc between susceptible and resilient groups, without control group.
Conclusion
Overall, we observed that only a subset of rats exposed to a traumatic event exhibited symptoms similar to those of PTSD. Our finding has facilitated the identification of alterations in the transcriptome of the mPFC, providing novel insights into the underlying mechanisms that may account for the various phenotypic effects induced by stress. The investigation of gene expression alterations in the IL and PL and the pathway that underlie these changes will enhance our understanding of the core neurobiological processes of PTSD. Although multiple gene mutations likely contribute to PTSD in different individuals, there may be converging neurobiological gene expression changes that contribute to the divergent symptoms of PTSD.
Data availability statement
The data presented in this study are deposited in the NCBI SRA database with accession number PRJNA1001887 and are accessible at the following link: https://www.ncbi.nlm.nih.gov/sra/PRJNA1001887.
Ethics statement
The animal study was approved by the Biomedical Ethics Committee of Peking University. The study was conducted in accordance with the local legislation and institutional requirements.
Author contributions
XJ, MY, and YX were instrumental in the conception and design of the work. GH, JI, and DS were involved in the acquisition of the data. GH and JI were responsible for analyzing and interpreting the data, as well as writing the article. In the course of a rigorous review process, XJ, MY, and YX have collaborated to scrutinize the article and have committed to being answerable for all elements of the project. All authors contributed to the article and approved the submitted version.
Funding
This work was supported by grants from the Science and Technology Planning Project of Shenzhen Municipality (KCXFZ20211020164543007, 20210617155253001), the Guangdong Basic and Applied Basic Research Foundation (2021A1515012141), the Shenzhen Fund for Guangdong Provincial High-level Clinical Key Specialties (SZGSP013), the Fund of Development and Reform Commission of Shenzhen Municipality (XMHT20220104028), and the National Natural Science Foundation of China (82071498, 82101492 and 32250410298).
Acknowledgments
The authors thank all participants for their contribution.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Footnotes
References
1. Shalev, A , Liberzon, I , and Marmar, C . Post-traumatic stress disorder. N Engl J Med. (2017) 376:2459–69. doi: 10.1056/NEJMra1612499
2. Kessler, RC , Chiu, WT , Demler, O , Merikangas, KR , and Walters, EE . Prevalence, severity, and comorbidity of 12-month DSM-IV disorders in the National Comorbidity Survey Replication. Arch Gen Psychiatry. (2005) 62:617–27. doi: 10.1001/archpsyc.62.6.617
3. Cathomas, F , Murrough, JW , Nestler, EJ , Han, M-H , and Russo, SJ . Neurobiology of resilience: interface between mind and body. Biol Psychiatry. (2019) 86:410–20. doi: 10.1016/j.biopsych.2019.04.011
4. Horn, SR , Charney, DS , and Feder, A . Understanding resilience: new approaches for preventing and treating PTSD. Exp Neurol. (2016) 284:119–32. doi: 10.1016/j.expneurol.2016.07.002
5. Russo, SJ , Murrough, JW , Han, M-H , Charney, DS , and Nestler, EJ . Neurobiology of resilience. Nat Neurosci. (2012) 15:1475–84. doi: 10.1038/nn.3234
6. Hare, BD , and Duman, RS . Prefrontal cortex circuits in depression and anxiety: contribution of discrete neuronal populations and target regions. Mol Psychiatry. (2020) 25:2742–58. doi: 10.1038/s41380-020-0685-9
7. Sotres-Bayon, F , Bush, DEA , and LeDoux, JE . Emotional perseveration: an update on prefrontal-amygdala interactions in fear extinction. Learn Mem. (2004) 11:525–35. doi: 10.1101/lm.79504
8. Alexandra Kredlow, M , Fenster, RJ , Laurent, ES , Ressler, KJ , and Phelps, EA . Prefrontal cortex, amygdala, and threat processing: implications for PTSD. Neuropsychopharmacology. (2022) 47:247–59. doi: 10.1038/s41386-021-01155-7
9. Fenster, RJ , Lebois, LAM , Ressler, KJ , and Suh, J . Brain circuit dysfunction in post-traumatic stress disorder: from mouse to man. Nat Rev Neurosci. (2018) 19:535–51. doi: 10.1038/s41583-018-0039-7
10. Yabuki, Y , and Fukunaga, K . Clinical therapeutic strategy and neuronal mechanism underlying post-traumatic stress disorder (PTSD). Int J Mol Sci. (2019) 20:3614. doi: 10.3390/ijms20153614
11. Nawreen, N , Baccei, ML , and Herman, JP . Single prolonged stress reduces intrinsic excitability and excitatory synaptic drive onto pyramidal neurons in the Infralimbic prefrontal cortex of adult male rats. Front Cell Neurosci. (2021) 15:705660. doi: 10.3389/fncel.2021.705660
12. Piggott, VM , Bosse, KE , Lisieski, MJ , Strader, JA , Stanley, JA , Conti, AC, et al. Single-prolonged stress impairs prefrontal cortex control of amygdala and striatum in rats. Front Behav Neurosci. (2019) 13:18. doi: 10.3389/fnbeh.2019.00018
13. Tao, Y , Cai, C-Y , Xian, J-Y , Kou, X-L , Lin, Y-H , Qin, C, et al. Projections from Infralimbic cortex to paraventricular thalamus mediate fear extinction retrieval. Neurosci Bull. (2021) 37:229–41. doi: 10.1007/s12264-020-00603-6
14. Covington, HE , Lobo, MK , Maze, I , Vialou, V , Hyman, JM , Zaman, S, et al. Antidepressant effect of optogenetic stimulation of the medial prefrontal cortex. J Neurosci. (2010) 30:16082–90. doi: 10.1523/JNEUROSCI.1731-10.2010
15. Heidbreder, CA , and Groenewegen, HJ . The medial prefrontal cortex in the rat: evidence for a dorso-ventral distinction based upon functional and anatomical characteristics. Neurosci Biobehav Rev. (2003) 27:555–79. doi: 10.1016/j.neubiorev.2003.09.003
16. Girgenti, MJ , Hare, BD , Ghosal, S , and Duman, RS . Molecular and cellular effects of traumatic stress: implications for PTSD. Curr Psychiatry Rep. (2017) 19:85. doi: 10.1007/s11920-017-0841-3
17. Kular, L , and Kular, S . Epigenetics applied to psychiatry: clinical opportunities and future challenges. Psychiatry Clin Neurosci. (2018) 72:195–211. doi: 10.1111/pcn.12634
18. Kuehner, JN , Bruggeman, EC , Wen, Z , and Yao, B . Epigenetic regulations in neuropsychiatric disorders. Front Genet. (2019) 10:268. doi: 10.3389/fgene.2019.00268
19. Wang, J , Chen, J , and Sen, S . MicroRNA as biomarkers and diagnostics. J Cell Physiol. (2016) 231:25–30. doi: 10.1002/jcp.25056
20. Girgenti, MJ , and Duman, RS . Transcriptome alterations in posttraumatic stress disorder. Biol Psychiatry. (2018) 83:840–8. doi: 10.1016/j.biopsych.2017.09.023
21. Saliminejad, K , Khorram Khorshid, HR , Soleymani Fard, S , and Ghaffari, SH . An overview of microRNAs: biology, functions, therapeutics, and analysis methods. J Cell Physiol. (2019) 234:5451–65. doi: 10.1002/jcp.27486
22. Bushati, N , and Cohen, SM . microRNA functions. Annu Rev Cell Dev Biol. (2007) 23:175–205. doi: 10.1146/annurev.cellbio.23.090506.123406
23. Dias, BG , Goodman, JV , Ahluwalia, R , Easton, AE , Andero, R , and Ressler, KJ . Amygdala-dependent fear memory consolidation via miR-34a and notch signaling. Neuron. (2014) 83:906–18. doi: 10.1016/j.neuron.2014.07.019
24. Chen, Y-L , Tong, L , Chen, Y , Fu, C-H , Peng, J-B , and Ji, L-L . MiR-153 downregulation alleviates PTSD-like behaviors and reduces cell apoptosis by upregulating the Sigma-1 receptor in the hippocampus of rats exposed to single-prolonged stress. Exp Neurol. (2022) 352:114034. doi: 10.1016/j.expneurol.2022.114034
25. Chen, Y , An, Q , Yang, S-T , Chen, Y-L , Tong, L , and Ji, L-L . MicroRNA-124 attenuates PTSD-like behaviors and reduces the level of inflammatory cytokines by downregulating the expression of TRAF6 in the hippocampus of rats following single-prolonged stress. Exp Neurol. (2022) 356:114154. doi: 10.1016/j.expneurol.2022.114154
26. Schmidt, U , Herrmann, L , Hagl, K , Novak, B , Huber, C , Holsboer, F, et al. Therapeutic action of fluoxetine is associated with a reduction in prefrontal cortical miR-1971 expression levels in a mouse model of posttraumatic stress disorder. Front Psych. (2013) 4:66. doi: 10.3389/fpsyt.2013.00066
27. Balakathiresan, NS , Chandran, R , Bhomia, M , Jia, M , Li, H , and Maheshwari, RK . Serum and amygdala microRNA signatures of posttraumatic stress: fear correlation and biomarker potential. J Psychiatr Res. (2014) 57:65–73. doi: 10.1016/j.jpsychires.2014.05.020
28. Martin, CG , Kim, H , Yun, S , Livingston, W , Fetta, J , Mysliwiec, V, et al. Circulating miRNA associated with posttraumatic stress disorder in a cohort of military combat veterans. Psychiatry Res. (2017) 251:261–5. doi: 10.1016/j.psychres.2017.01.081
29. Zhou, J , Nagarkatti, P , Zhong, Y , Ginsberg, JP , Singh, NP , Zhang, J, et al. Dysregulation in microRNA expression is associated with alterations in immune functions in combat veterans with post-traumatic stress disorder. PLoS One. (2014) 9:e94075. doi: 10.1371/journal.pone.0094075
30. Fang, Q , Li, Z , Huang, G-D , Zhang, H-H , Chen, Y-Y , Zhang, L-B, et al. Traumatic stress produces distinct activations of GABAergic and glutamatergic neurons in amygdala. Front Neurosci. (2018) 12:387. doi: 10.3389/fnins.2018.00387
31. Nahvi, RJ , Tanelian, A , Nwokafor, C , Godino, A , Parise, E , Estill, M, et al. Transcriptome profiles associated with resilience and susceptibility to single prolonged stress in the locus coeruleus and nucleus accumbens in male Sprague-dawley rats. Behav Brain Res. (2023) 439:114162. doi: 10.1016/j.bbr.2022.114162
32. Richter-Levin, G , Stork, O , and Schmidt, MV . Animal models of PTSD: a challenge to be met. Mol Psychiatry. (2019) 24:1135–56. doi: 10.1038/s41380-018-0272-5
33. Daskalakis, NP , Yehuda, R , and Diamond, DM . Animal models in translational studies of PTSD. Psychoneuroendocrinology. (2013) 38:1895–911. doi: 10.1016/j.psyneuen.2013.06.006
34. Souza, RR , Noble, LJ , and McIntyre, CK . Using the single prolonged stress model to examine the pathophysiology of PTSD. Front Pharmacol. (2017) 8:615. doi: 10.3389/fphar.2017.00615
35. Liberzon, I , Krstov, M , and Young, EA . Stress-restress: effects on ACTH and fast feedback. Psychoneuroendocrinology. (1997) 22:443–53. doi: 10.1016/S0306-4530(97)00044-9
36. Török, B , Sipos, E , Pivac, N , and Zelena, D . Modelling posttraumatic stress disorders in animals. Prog Neuro-Psychopharmacol Biol Psychiatry. (2019) 90:117–33. doi: 10.1016/j.pnpbp.2018.11.013
37. Le Dorze, C , and Gisquet-Verrier, P . Sensitivity to trauma-associated cues is restricted to vulnerable traumatized rats and reinstated after extinction by yohimbine. Behav Brain Res. (2016) 313:120–34. doi: 10.1016/j.bbr.2016.07.006
38. Tanelian, A , Nankova, B , Miari, M , Nahvi, RJ , and Sabban, EL . Resilience or susceptibility to traumatic stress: potential influence of the microbiome. Neurobiol Stress. (2022) 19:100461. doi: 10.1016/j.ynstr.2022.100461
39. Del Giudice, M , Ellis, BJ , and Shirtcliff, EA . The adaptive calibration model of stress responsivity. Neurosci Biobehav Rev. (2011) 35:1562–92. doi: 10.1016/j.neubiorev.2010.11.007
40. Drevets, WC , Price, JL , and Furey, ML . Brain structural and functional abnormalities in mood disorders: implications for neurocircuitry models of depression. Brain Struct Funct. (2008) 213:93–118. doi: 10.1007/s00429-008-0189-x
41. Arnsten, AFT , Raskind, MA , Taylor, FB , and Connor, DF . The effects of stress exposure on prefrontal cortex: translating basic research into successful treatments for post-traumatic stress disorder. Neurobiol Stress. (2015) 1:89–99. doi: 10.1016/j.ynstr.2014.10.002
42. Chen, F , Ke, J , Qi, R , Xu, Q , Zhong, Y , Liu, T, et al. Increased inhibition of the amygdala by the mPFC may reflect a resilience factor in post-traumatic stress disorder: a resting-state fMRI granger causality analysis. Front Psych. (2018) 9:516. doi: 10.3389/fpsyt.2018.00516
43. Milad, MR , Rauch, SL , Pitman, RK , and Quirk, GJ . Fear extinction in rats: implications for human brain imaging and anxiety disorders. Biol Psychol. (2006) 73:61–71. doi: 10.1016/j.biopsycho.2006.01.008
44. Omura, J , Fuchikami, M , Araki, M , Miyagi, T , Okamoto, Y , and Morinobu, S . Chemogenetic activation of the mPFC alleviates impaired fear memory extinction in an animal model of PTSD. Prog Neuro-Psychopharmacol Biol Psychiatry. (2021) 108:110090. doi: 10.1016/j.pnpbp.2020.110090
45. Sierra-Mercado, D , Padilla-Coreano, N , and Quirk, GJ . Dissociable roles of prelimbic and infralimbic cortices, ventral hippocampus, and basolateral amygdala in the expression and extinction of conditioned fear. Neuropsychopharmacology. (2011) 36:529–38. doi: 10.1038/npp.2010.184
46. Norrholm, SD , Jovanovic, T , Olin, IW , Sands, LA , Karapanou, I , Bradley, B, et al. Fear extinction in traumatized civilians with posttraumatic stress disorder: relation to symptom severity. Biol Psychiatry. (2011) 69:556–63. doi: 10.1016/j.biopsych.2010.09.013
47. Kim, GS , Smith, AK , Nievergelt, CM , and Uddin, M . Neuroepigenetics of post-traumatic stress disorder. Prog Mol Biol Transl Sci. (2018) 158:227–53. doi: 10.1016/bs.pmbts.2018.04.001
48. Breslau, N , Peterson, EL , Poisson, LM , Schultz, LR , and Lucia, VC . Estimating post-traumatic stress disorder in the community: lifetime perspective and the impact of typical traumatic events. Psychol Med. (2004) 34:889–98. doi: 10.1017/S0033291703001612
49. Nie, P-Y , Tong, L , Li, M-D , Fu, C-H , Peng, J-B , and Ji, L-L . miR-142 downregulation alleviates rat PTSD-like behaviors, reduces the level of inflammatory cytokine expression and apoptosis in hippocampus, and upregulates the expression of fragile X mental retardation protein. J Neuroinflammation. (2021) 18:17. doi: 10.1186/s12974-020-02064-0
Keywords: transcriptome, stress, susceptibility, resilience, prefrontal cortex
Citation: Huang G, Iqbal J, Shen D, Xue Y-x, Yang M and Jia X (2023) MicroRNA expression profiles of stress susceptibility and resilience in the prelimbic and infralimbic cortex of rats after single prolonged stress. Front. Psychiatry. 14:1247714. doi: 10.3389/fpsyt.2023.1247714
Edited by:
Fushun Wang, Nanjing University of Chinese Medicine, ChinaReviewed by:
Wenhua Zhang, Nanchang University, ChinaZili You, University of Electronic Science and Technology of China, China
Copyright © 2023 Huang, Iqbal, Shen, Xue, Yang and Jia. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xiaojian Jia, amlheGlhb2ppYW5Ab3V0bG9vay5jb20=; Mei Yang, eW04MzQyQDE2My5jb20=; Yan-xue Xue, eWFueHVleHVlQGJqbXUuZWR1LmNu
†These authors have contributed equally to this work and share first authorship