- Micro report
- Open Access
Identification of stress resilience module by weighted gene co-expression network analysis in Fkbp5-deficient mice
Molecular Brain volume 12, Article number: 99 (2019)
FKBP5 encodes the FK506 binding protein 5, a glucocorticoid receptor (GR) binding protein known to play an important role in the physiological stress response. However, results from previous studies examining the association between common variants of FKBP5 and stress have been inconsistent. To investigate whether the loss of FKBP5 affects the stress response, we examined the behavior of mice following the induction of chronic restraint stress between homozygous wild-type and Fkbp5 knock-out mice. After 21 days of exposure to restraint stress, WT mice showed anhedonia, a core symptom of depression, which could be measured by a sucrose preference test. However, Fkbp5-deficient mice did not exhibit significant depressive-like behavior compared to the WT after exposure to chronic restraint stress. To investigate the molecular mechanism underlying stress resilience, we performed RNA sequencing analysis. The differentially expressed gene (DEG) analysis showed that chronic stress induced changes in various biological processes involved in cell-cell adhesion and inflammatory response. Weighted gene co-expression network analysis identified 60 characteristic modules that correlated with stress or the FKBP5 genotype. Among them, M55 showed a gene expression pattern consistent with behavioral changes after stress exposure, and the gene ontology analysis revealed that this was involved in nervous system development, gland morphogenesis, and inflammatory response. These results suggest that FKBP5 may be a crucial factor for the stress response, and that transcriptomic data can provide insight into stress-related pathophysiology.
Depression is one of the most common mental disorders affecting people of all ages, and can arise from a variety of causes, including genetic susceptibility, endocrine dysregulation, and stresses in life . When exposed to acute and temporary stress, while the body can protect itself from stress, chronic stress can disturb the function of the brain system. Chronic accumulation of stress leads to abnormal and excessive cortisol secretion in the hypothalamic-pituitary-adrenal axis, which affects a variety of physical activities, including brain function, leading to mental disorders such as depression or post-traumatic stress disorder (PTSD) . Proper regulation through negative feedback of glucocorticoid receptors (GR) is important for the stress response, and long-term or excessive activation of this system is associated with the development of depression or anxiety disorders . The FK506-binding protein 51 (FKBP5) is a co-chaperone of Hsp90 in the GR molecular complex and is a key modulator of GR sensitivity [4, 5]. Although FKBP5 is an important factor that is responsible for coping behavior as well as neuroendocrine responses , results from previous studies investigating the association between Fkbp5 gene variants and stress remain controversial. In order to investigate whether genetic FKBP5 variants affect behavior in response to chronic restraint stress exposure, we examined the behavior of mice following the induction of chronic restraint stress in homozygous wild-type (WT) and knock-out (KO) mice. After 21 days of exposure to restraint stress, while WT mice showed anhedonia in the sucrose preference test, Fkbp5-deficient mice did not exhibit significant depressive-like behavior compared to the WT (Fig. 1)a and b.
Recent studies performed in humans and rodents have suggested that long-term stress and pathological anxiety leads to structural degeneration and functional alteration of the frontal cortex, and increases the risk of mental disorders [7, 8]. In addition, it has been suggested that the medial prefrontal cortex (mPFC), a region controlling higher brain function including cognition and emotion, is a primary target of stress hormones [9, 10]. However, little is known about the molecular mechanisms in the mPFC involved in stress-associated psychiatric disorders. Transcriptome profiling has helped provide an unbiased insight into the pathophysiological mechanism underlying complicated brain disorders . Therefore, using RNA sequencing (RNAseq) analysis, we investigated the dynamic transcriptomic changes that occur after stress in the mPFC of Fkbp5-deficient mice.
Multidimensional scaling analysis showed a clear separation between the stressed and the non-stressed control (CT) mice. There was also a slight overlap between the CRS-exposed WT mice (WT_ST) and CRS-exposed Fkbp5 KO mice (KO_ST) groups, which were clustered according to their genotype (Fig. 1)c. To identify the genotypes and genes affected by stress, we analyzed differentially expressed genes (DEG) between each group (Additional file 1). Of the 24,532 mRNA genes profiled, 224 (0.91%; CT vs WT_ST), 258 (1.05%; CT vs KO_ST), and 135 (0.55%; WT_ST vs KO_ST) genes were dysregulated in the mPFC following chronic restraint stress, and the percentage of DEG induced by stress was higher than the percentage of DEG by genotype (Additional file 2: Figure S1, Additional file 3: Table S1-S3). Gene ontology (GO) enrichment analysis of these expression profiles showed functional categories that were potentially dysregulated, including lipid metabolic process (Benjamini adjusted p = 9.63 × 10− 4), regulation of immune response (Benjamini adjusted p = 7.24 × 10− 4), cell adhesion (Benjamini adjusted p = 1.52 × 10− 3), regulation of cell differentiation (Benjamini adjusted p = 1.17 × 10− 3), and neurogenesis (Benjamini adjusted p = 7.66 × 10− 2) (Additional file 4: Table S4-S6).
To compare the differences in expression observed in the multidimensional data set with the pattern of stress-response behavior, we performed a weighted gene co-expression network analysis (WGCNA) . Through WGCNA, we identified 60 modules of co-expressed genes following chronic restraint stress in both WT and KO homozygous genotypes (Additional file 5: Figure S2, Additional file 6: Table S7). Among the 60 modules, we identified characteristic modules which showed a significant correlation to the genotype (M3, M6, M25, M33, M39, M44 and M57) and to the stress exposure (M3, M10, M18, M21, M31, M34, M50 and M55) (Fig. 1, Additional file 7: Figure S3)d, a-c. Interestingly, one of the modules that negatively correlated with stress, M55, had a pattern similar to the stress resilience behavior of Fkbp5-deficient mice. This was down-regulated, consistent with depression-like behavior in the WT_ST group, and was restored in the KO_ST group (Fig. 1)d and e. GO enrichment analysis of M55 revealed the biological functions potentially involved, including gland morphogenesis (p = 1.25 × 10− 2), activation of immune response (p = 2.25 × 10− 2), and nervous system development (p = 7.65 × 10− 3) (Fig. 1, Additional file 8: Table S8)f and g. Hub genes, with the highest degree of connectivity within a module of the M55 include Mmp25. This gene has been functionally implicated in the regulation of immune response through NF-B signaling  and has been linked to neuropsychiatric disorders including PTSD .
In this study, we compared brain transcriptome altered by chronic stress in the mPFC between Fkbp5-deficient and wild-type mice by RNAseq analysis. In addition to the DEG analysis, by employing WGCNA, we identified a distinct co-expression network module associated with stress resilience caused by Fkbp5 knock-out, and characterized the biological processes affected by this module, leading to this unique behavior. Our systematic transcriptome analysis demonstrated that aberration in the development of the neuroendocrine system, and regulation of the immune response may underlie the stress resilient behavior observed in the Fkbp5 deficient mice. This is the first study, to the best of our knowledge, to identify the stress resilience associated genes through gene co-expression network analysis in Fkbp5 deficient mice. Altogether, we confirmed that FKBP5 may be an important component of the stress response, suggesting that identification of the module associated with the stress response can provide a treatment strategy and therapeutic target to attenuate the depressive symptoms caused by stress.
Availability of data and materials
The datasets supporting the conclusions of this article are included in this published article and its additional files. RNAseq raw data reported in this paper are submitted to the GEO repository under the accession number GSE138240.
Benjamini-Hochberg adjusted P-value
Chronic restraint stress
Differentially expressed genes
CRS-exposed Fkbp5 KO mice
Medial prefrontal cortex
Post-traumatic stress disorder
Weighted gene co-expression network analysis
CRS-exposed WT mice
Jabbi M, Korf J, Kema IP, Hartman C, van der Pompe G, Minderaa RB, et al. Convergent genetic modulation of the endocrine stress response involves polymorphic variations of 5-HTT. COMT MAOA Mol Psychiatry. 2007;12:483.
de Kloet CS, Vermetten E, Geuze E, Kavelaars A, Heijnen CJ, Westenberg HG. Assessment of HPA-axis function in posttraumatic stress disorder: pharmacological and non-pharmacological challenge tests, a review. J Psychiatr Res. 2006;40:550.
Lucassen PJ, Pruessner J, Sousa N, Almeida OF, Van Dam AM, Rajkowska G, et al. Neuropathology of stress. Acta Neuropathol. 2014;127:109.
Wiederrecht G, Hung S, Chan HK, Marcy A, Martin M, Calaycay J, et al. Characterization of high molecular weight FK-506 binding activities reveals a novel FK-506-binding protein as well as a protein complex. J Biol Chem. 1992;267:21753.
Zannas AS, Wiechmann T, Gassen NC, Binder EB. Gene-stress-epigenetic regulation of FKBP5: clinical and translational implications. Neuropsychopharmacol. 2016;41:261.
Touma C, Gassen NC, Herrmann L, Cheung-Flynn J, Bull DR, Ionescu IA, et al. FK506 binding protein 5 shapes stress responsiveness: modulation of neuroendocrine reactivity and coping behavior. Biol Psychiatry. 2011;70:928.
Korgaonkar MS, Fornito A, Williams LM, Grieve SM. Abnormal structural networks characterize major depressive disorder: a connectome analysis. Biol Psychiatry. 2014;76:567.
Laine MA, Sokolowska E, Dudek M, Callan SA, Hyytia P, Hovatta I. Brain activation induced by chronic psychosocial stress in mice. Sci Rep. 2017;7:15061.
Arnsten AF. Stress signalling pathways that impair prefrontal cortex structure and function. Nat Rev Neurosci. 2009;10:410.
Broadbent DE. Decis Stress. London: Academic; 1971.
Girgenti MJ, Duman RS. Transcriptome alterations in posttraumatic stress disorder. Biol Psychiatry. 2018;83:840.
Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. 2008;9:559.
Soria-Valles C, Gutierrez-Fernandez A, Osorio FG, Carrero D, Ferrando AA, Colado E, et al. MMP-25 Metalloprotease regulates innate immune response through NF-kappaB signaling. J Immunol. 2016;197:296.
Bam M, Yang X, Zumbrun EE, Zhong Y, Zhou J, Ginsberg JP, et al. Dysregulated immune system networks in war veterans with PTSD is an outcome of altered miRNA expression and DNA methylation. Sci Rep. 2016;6:31209.
This research was supported by the Chung-Ang University Scholarship Grants in 2016 and the Brain Research Program of the National Research Foundation of Korea (NRF) funded by the Ministry of Science and ICT (NRF-2018M3C7A1024150).
The animal experiments were conducted in accordance with the Chung-Ang University Council Directive for the use and care of laboratory animals, and National Institutes of Health animal care guidelines. The Institutional Animal Care and Use Committee of Chung-Ang University approved all protocols described in this study.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Materials and Methods.
Heatmap representing the expression profiles of the DEGs in the three groups of CT, WT_ST, and KO_ST.
List of DEG between WT ST and CT. Table S2. List of DEG between KO ST and CT. Table S3. List of DEG between KO ST and WT ST.
GO analysis for DEGs between WT ST and CT. Table S5. GO analysis for DEGs between KO ST and CT. Table S6. GO analysis for DEGs between KO ST and WT ST.
Heatmap of the correlation of WGCNA modules with traits. The correlation between each module eigengene and sample trait was calculated. Values in the figure indicate the correlation coefficient between modules and traits. Values in brackets are the p-values for the association test.
List of genes within each WGCNA module (kME > 0.7).
The interleaved scatter plots of module which have a significant correlation with genotype and stress (a-b) Modules with positively (a) and negatively (b) correlated to the genotype. (c) Modules with positively correlated to the stress. Data represent mean ± SEM. One-way ANOVA; Fisher’s LSD: *p ≤ 0.05, **p < 0.01, ***p < 0.001.
GO analysis for genes in M55.
About this article
Cite this article
Kwon, J., Kim, Y.J., Choi, K. et al. Identification of stress resilience module by weighted gene co-expression network analysis in Fkbp5-deficient mice. Mol Brain 12, 99 (2019). https://doi.org/10.1186/s13041-019-0521-9
- Chronic stress