Cheng Qin
Yujia Wang
Yang Zhang
Yan Zhu
Yabin Wang
Feng Cao- 1Department of Cardiology, National Clinical Research Center for Geriatric Diseases and Second Medical Center of Chinese PLA General Hospital, Beijing, China
- 2Nankai University School of Medicine, Nankai University, Tianjin, China
- 3Beijing Key Laboratory of Research on Aging and Related Diseases, Beijing, China
Background: Growing evidence has supported that chronic psychological stress would cause heart damage, However the mechanisms involved are not clear and effective interventions are insufficient. Cannabinoid type 2 receptor (CB2R) can be a potential treatment for cardiac injury. This study is aimed to investigate the protective mechanism of CB2R agonist against chronic psychological stress-induced cardiac injury.
Methods: A mouse chronic psychological stress model was constructed based on a chronic unpredictable stress pattern. Mice were performed a three-week psychological stress procedure, and cardiac tissues of them were collected for whole-transcriptome sequencing. Overlap analysis was performed on differentially expressed mRNAs (DE-mRNAs) and ER stress-related genes (ERSRGs), and bioinformatic methods were used to predict the ceRNA networks and conduct pathway analysis. The expressions of the DE-ERSRGs were validated by RT-qPCR.
Results: In the comparison of DE mRNA in Case group, Control group and Treatment group, three groups of ceRNA networks and ceRNA (circ) networks were constructed. The DE-mRNAs were mainly enriched in chromatid-relevant terms and Hematopoietic cell lineage pathway. Additionally, 13 DE-ERSRGs were obtained by the overlap analysis, which were utilized to establish a ceRNA network with 15 nodes and 14 edges and a ceRNA (circ) network with 23 nodes and 28 edges. Furthermore, four DE-ERSRGs (Cdkn1a, Atf3, Fkbp5, Gabarapl1) in the networks were key, which were mainly enriched in response to extracellular stimulus, response to nutrient levels, cellular response to external stimulus, and FoxO signaling pathway. Finally, the RT-qPCR results showed almost consistent expression patterns of 13 DE-ERSRGs between the transcriptome and tissue samples.
Conclusion: The findings of this study provide novel insights into the molecular mechanisms of chronic psychological stress-induced cardiac diseases and reveal novel targets for the cardioprotective effects of CB2R agonists.
1 Introduction
Psychological stress is a maladaptation in response to exposure to an environmental threat or severe external stimuli. Traumatic events, such as attacks, natural disasters, or military combat, can result in psychological stress, which is detrimental to the health of an individual (Glaser and Kiecolt-Glaser, 2005). A potential link might exist between psychological stress and cardiovascular diseases (CVD), a topic that has received increasing attention from scholars (Gold et al., 2020). Findings from multiple clinical studies have shown that psychological stress can increase the risk of adverse cardiovascular events, such as myocardial ischemia (Ebrahimi et al., 2021; Vaccarino et al., 2021; Shah et al., 2022), sudden death (Meisel et al., 1991; Leor et al., 1996), arrhythmia (Steinberg et al., 2004; Hourani et al., 2020), stroke (Jordan et al., 2013; Chen et al., 2015), and the aggravation of heart failure (Roy et al., 2015).
However, despite the epidemiological evidence, the pathophysiological mechanisms underlying these associations remain incomplete. Findings from recent studies suggest that autonomic dysfunction and inflammatory response promote plaque vulnerability, which is the cause of cardiovascular events induced by acute stress (Hinterdobler et al., 2021). However, the long-term effects of chronic stress on the cardiovascular system are not completely understood. Therefore, specific and effective interventions are urgently needed. Moreover, whether healthy individuals exposed to stress elicit similar pathophysiological responses as individuals with pre-existing cardiovascular diseases remains doubtful.
The endoplasmic reticulum (ER) is the primary site for protein synthesis, folding, and transport (Wang and Kaufman, 2016). The exposure of cells to oxidative stress, hypoxia, nutrient deprivation, and pathological infection will cause the accumulation of misfolded or unfolded proteins in the ER and initiate the unfolded protein response (UPR) (Ochoa et al., 2018). If stimuli persist, the UPR will lead to apoptosis and autophagy. ATF6, PERK, and IRE1 are the representative molecules of UPR signaling (Ron and Walter, 2007). Under usual conditions, these molecules bind to their partner molecule glucose-regulated protein 78 (GRP78). Under ER stress, unfolded or misfolded proteins recruit GRP78/BiP to dissociate from the complex composed of IRE1, ATF6, and PERK. This process marks the initiation of ER stress (Bertolotti et al., 2000; Frakes and Dillin, 2017). The expression of the abovementioned ER stress-related genes (ESRGs) was found to be related to the occurrence and development of various cardiovascular diseases, including myocardial ischemia (Bi et al., 2018; Zhang et al., 2019), ventricular remodeling (Blackwood et al., 2019), atherosclerosis (Tang et al., 2019), and vascular calcification (Fu et al., 2019). Likewise, the expression of ESRGs was found to be upregulated in the brain tissues of mice subjected to psychological stress (Gold et al., 2013; Mondal et al., 2015; Mao et al., 2019). Thus far, small molecule drugs targeting ESRGs have been developed, but they have not been adopted widely in clinical practice owing to their relative toxicity and lack of specificity (Atkins et al., 2013; Almanza et al., 2019).
Cannabinoid receptor type 2 (CB2R) is primarily localized in immune cells as well as in endothelial cells and cardiac tissues. It plays an important role in the regulation of the immune response and inflammation. CB2R agonists can exhibit a cardioprotective role in myocardial infarction, cardiomyopathy, arrhythmia, and other CVDs by reducing reactive oxygen species (ROS) production, inhibiting necroptosis and apoptosis, promoting efferocytosis, upregulating endothelial NO synthase, and reducing macrophage infiltration (Molica et al., 2012; Li et al., 2016; Yu et al., 2019; Liu X. et al., 2020). Owing to the anti-oxidative stress effect of CB2R agonists and the regulation of autophagy, CB2R agonists may affect the cellular process underlying ER stress.
Here, we aimed to elucidate the mechanism underlying cardiac injury induced by mental stress and the therapeutic targets of CB2R agonists. We constructed a mouse model induced by psychological stress (case group), stressed mice treated with CB2R agonists (treatment group), and normal mice (control group) and collected samples for whole transcriptome sequencing analysis to analyze the differences in gene expression profiles between the different groups.
2 Methods
2.1 Animal models
All mice were obtained from SPF Biotechnology and allowed 1 week of acclimation to the specific-pathogen-free animal facilities before the experiments were conducted. Six-week-old male mice were used in this study and were randomly divided into three groups: control group (three normal mice were maintained with free access to food and water on a 12 h light/dark cycle), case group (four mice with psychological stress-induced heart disease), and treatment group (four mice with psychological stress-induced heart disease treated with the CB2R agonist JWH133). JWH133 (ApexBio, B7941) was injected twice daily intraperitoneally at 10 mg/kg from the beginning of the stress induction procedure for 21 days. The stress induction process is shown inSupplementary Table S1. Behavioral tests were performed 24 h after stress induction.
2.2 Data sources
After stress induction, the hearts of the mice were removed, and the tissue was processed for whole-transcriptome sequencing analysis. After RNA extraction, purification, and library construction, these libraries were double-ended (paired-end, PE) sequenced by next-generation sequencing (NGS) using the Illumina HiSeq sequencing platform. The sequencing depth for the transcriptome-wide analysis was showed in Supplementary Figure S1.
In addition, 787 ERS-related genes with a relevance score of ≥7 were extracted from a report by Zhang et al. (2021). Since these genes were from a human gene set, the biomaRt package in R software was used to identify the homologous genes of mice. We identified 828 ERS-related genes in mice (Supplementary Table S2). GSE210252 was downloaded from the Gene Expression Omnibus (GEO) dataset (https://www.ncbi.nlm.nih.gov/), which contained four pairs samples from mice received restraint stress or no stress.
2.3 Behavioral test
2.3.1 Open-field test
The open-field test can be used to observe voluntary movement and assess the severity of stress. The open-field system comprised an open-field box (40 × 40 × 40 cm3), an imaging system, and analysis software. During the experiment, the mice were placed in a specific position in the open-field, the camera system was used to monitor the activities of the animals in the open-field, and an analysis software was used to track the trajectory of the animals and collect data on their length of stay in the central area. The lesser time the mice remained in the central area, the more severe stress they experienced (Yan et al., 2022).
2.3.2 Elevated plus maze (EPM)
EPM is an experimental tool used to evaluate the stress response of rodents. In this study, EPM was a cross-shaped platform composed of two open arms and two closed arms (arms 35 cm long × 5 cm wide) positioned approximately 55 cm above the ground. A monitoring camera for animal behavior experiments was installed above the EPM. The residence time in the open arm (open arm time, OT) was recorded using the Labmaze 3.0 animal behavior analysis software (Zhongshidichuang Sci. and Tech). For the EPM test, animals were placed in the center of the maze facing an open arm (Cho et al., 2022).
2.3.3 Forced swimming test
The forced swimming test was conducted in a cylindrical water tank (50 cm in height and 10 cm in diameter). The water temperature was 23 ∼ 25 °C, and the water had a depth that was sufficient to position the tail of the animal at a certain distance from the bottom. The swimming time was 6 min, and the immobility time in the last 4 min of the swimming duration was recorded. Moreover, the standard of immobility was that the animal stopped struggling in the water and started floating, with only slight limb movement to keep the head afloat (Gu et al., 2021). A longer immobility time represented more severe stress.
2.3.4 Tail suspension test
This experiment was conducted in a dark and quiet environment. The tail of the mouse was fixed and hung in a 25 cm × 25 cm × 30 cm box, and the head of the mouse was positioned approximately 5 cm from the bottom. Furthermore, the activity of the mouse was recorded within 6 min. The immobility time in the last 4 min was recorded. Small movements of the forelimbs rather than the hindlimbs were considered to represent immobility (Gu et al., 2021).
2.4 Differential analysis
The Limma package (version 3.44.3) was used to screen differentially expressed mRNAs (DE-mRNAs), miRNAs (DE-miRNAs), lncRNAs (DE-lncRNAs), and circRNAs (DE-circRNAs) in the three comparison groups (Control vs. Case; Case vs. Treatment; Control vs. Treatment). In detail, the count values obtained by comparing the mRNA and lncRNA data were converted to the fragments per kilobase per million mapped reads (FPKM) values. The count values used by miRNA and circRNA were converted to the counts per million mapped (CPM) value, and the selection standard was set at p < .05 and |Log2FC|>1. Moreover, the distribution of DE-mRNAs, DE-miRNAs, DE-circRNAs, and DE-lncRNAs in the three comparison groups was demonstrated at the chromosomal level using Circos in R package.
2.5 Construction of ceRNA networks
To predict target miRNAs that could potentially bind to the DE-mRNAs in the three comparison groups, mirwalk 3.0 was used with a threshold score >.95. Simultaneously, only the DE-miRNAs with opposite expression trends to DE-mRNAs were retained. After the target miRNA was obtained, starBase 2.0 was used to predict the miRNA-lncRNA pairs, and only the relationship pairs with opposite differential expression trends were retained. Thus, an mRNA-miRNA-lncRNA binding network was obtained from each comparison group, which was considered the ceRNA network.
Similarly, in each comparison group, mirwalk was used to predict the binding of DE-mRNA and DE-miRNA with a threshold score >.95, and only the pairs with contrasting trends of differential expression were retained. To predict the circRNAs that might be related to the miRNAs after obtaining miRNAs, the circRNA sequences were extracted from assembled gene files, miRNA sequences of mice were retrieved from miRbase, and the miRanda software was used to predict the binding with the selection criteria as binding score = 140, and only the opposite relationship pairs of differential expression trends were retained. Thus, an mRNA-miRNA-circRNA network was constructed for each group, known as the ceRNA (circ) network.
2.6 Screening and function annotation of differentially expressed ER stress-related genes (DE-ERSRGs)
Overlap analysis was performed on the DE-mRNAs from the Case vs. Control and Treatment vs. Case comparison groups and the 828 mouse ER stress-related genes. The intersected genes were considered to be DE-ERSRGs.
2.7 Construction of the DE-ERSRGs-based ceRNA network
Initially, mirwalk was used to predict the target miRNAs that can potentially bind to the DE-ERSRGs with a threshold score >.95. Simultaneously, only miRNAs in the pairs with contrasting trends of differential expression were retained in the Treatment vs. Case group. After obtaining the target miRNA, starBase was applied to predict the miRNA-lncRNA pairs, and the relationship pairs with opposite differential expression trends were retained as well. Thus, an mRNA-miRNA-lncRNA binding network was obtained for the Treatment vs. Case group.
Likewise, mirwalk was used to predict the target miRNAs of DE-ERSRG, and only miRNAs in the miRNAs-DE-ERSRG pairs with opposite trends of differential expressions were retained. After the miRNAs were obtained, the target circRNAs of the miRNAs were predicted using miRanda, and the pairs with the contrasting relationship of differential expression trends were selected. An mRNA-miRNA-circRNA network was established using the data.
2.8 Functional enrichment analyses of DE-ERSRGs in the ceRNA network
Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed with DE-mRNAs, mRNAs in the ceRNA and ceRNA (circ) networks as well as DE-ERSRGs from the three comparison groups individually using clusterProfiler (version4.0.2), with p < .05 as the selection standard.
2.9 Real-time qPCR
Eleven frozen samples were segregated into three groups, with three samples each in the control and case groups and five samples in the treatment group. Then, 50 mg of tissue was extracted from each sample and lysed using TRIzol Reagent (Life Technologies, CA, United States), and total RNA was isolated by following the manufacturer’s instructions. After the concentration and purity of RNA were determined, the RNA was reverse-transcribed to cDNA using the SureScript-First-strand-cDNA-synthesis-kit (Genecopoeia, Guangzhou, China) before qRT-PCR. The qRT-PCR mixture was composed of 3 µL of cDNA, 5 µL of 2x Universal Blue SYBR Green qPCR Master Mix (Servicebio, Wuhan, China), and 1 µL each of forward and reverse primers. PCR was performed in a BIO-RAD CFX96 Touch TM PCR detection system (Bio-Rad Laboratories, Inc., United States) under the following thermal cycling conditions: 40 cycles at 95°C for 60 s, 95°C for 20 s, 55°C for 20 s, and 72°C for 30 s. The 2−△△Ct method was used to calculate gene expression levels, and GraphPad Prism 5 was applied to plot and calculate the statistical significance. The primer sequences used in this study are listed in Supplementary Table S3.
3 Results
3.1 Successful establishment of the chronic psychological stress model
Behavioral tests were performed 24 h after the end of stress induction in the three groups of mice. As shown in Supplementary Figure S2A, in the open-field experiment, mice from the case group spent significantly less time in the central zone than mice from the control and CB2R agonist treatment groups. Additionally, in the elevated plus maze test, the time of stay in the open arm in the case group was significantly lesser than that in the control group, but no difference was observed between the case and treatment groups (Supplementary Figure S2B). Furthermore, in the forced swimming test, the immobility time in the case group was significantly greater than that in the control and treatment groups (Supplementary Figure S2C). Finally, in the tail suspension test, the immobility time in the case group was significantly greater than that in the control group, but no statistically significant differences were observed compared with the treatment group (Supplementary Figure S2D). The above results indicated that the protocol could induce significant stress in the mice.
3.2 Identification of DE-mRNAs, DE-miRNAs, DE-lncRNAs, and DE-circRNAs
In the Control vs. Case group, 701 DE-mRNAs (357 upregulated, 344 downregulated), 80 DE-miRNAs (19 upregulated, 61 downregulated), 530 DE-lncRNAs (268 upregulated, 262 downregulated), and 95 DE-circRNAs (26 upregulated, 69 downregulated) were screened between four case mice and three normal case mice (Figure 1A). In the Treatment vs. Case group, 466 DE-mRNAs (230 upregulated, 236 downregulated), 18 DE-miRNAs (eight upregulated, ten downregulated), 500 DE-lncRNAs (290 upregulated, 210 downregulated), and 83 DE-circRNAs (45 upregulated, 38 downregulated) were obtained (Figure 1B). In the Control vs. Treatment group, 475 DE-mRNAs (231 upregulated, 245 downregulated), 24 DE-miRNAs (four upregulated, 20 downregulated), 349 DE-lncRNAs (196 upregulated, 153 downregulated), and 67 DE-circRNAs (27 upregulated, 40 downregulated) were eventually selected (Figure 1C).
FIGURE 1. Differential expression analysis. Volcano plot for differentially expressed mRNAs (DE-mRNAs), lncRNAs (DE-lncRNAs), miRNAs (DE-miRNAs) and circRNAs (DE-circRNAs) in the three comparison groups. (A) Control vs. Case; (B) Case vs. Treatment; (C) Control vs. Treatment. Blue: downregulation, red: upregulation.
In the Control vs. Case group, the majority of DE-mRNAs were located in the chromosomes, except the Y chromosome (only one downregulated DE-mRNA was located in Y chromosome). The DE-lncRNAs were almost evenly situated in each pair of chromosomes, and DE-miRNAs were situated in chromosomes 2, 3, 4, 7, 11, 12, 15, 17, 19, and X. DE-circRNAs were located in the euchromosomes and chromosome X (Supplementary Figure S3A).
In the Treatment vs. Case group, the DE-mRNAs, DE-lncRNAs, and DE-circRNAs were situated in all chromosomes except the Y chromosome. The DE-miRNAs were located in the 3, 5, 6, 7, 8, 11, 12, 13, 15, and X chromosomes (Supplementary Figure S3B).
In the Treatment vs. Control group, the locations of DE-mRNAs and DE-lncRNAs were similar to those in the Treatment vs. Case group, but the DE-circRNAs were located in the euchromosomes and Y chromosomes, and DE-miRNAs were located in the 1, 3, 5, 6, 7, 8, 10, 11, 12, 14, 15, 19, and X chromosomes (Supplementary Figure S3C).
3.3 Functional enrichment analyses of DE-mRNAs in the three groups
The functional enrichment results of DE-mRNAs in the Control vs. Case group illustrated that 1261 GO terms and 98 KEGG pathways were enriched by 701 DE-mRNAs. For example, chromosome separation, nuclear division, nuclear chromosome segregation, organelle fission, and mitotic nuclear segregation were the major enriched GO BP terms. Chromosomes, centromeric region, spindle, and chromosomal region, among others, were the primarily enriched terms in CC. In MF, microtubule binding, tubulin binding, metal ion transmembrane transporter activity, microtubule motor activity, among others, were the major common terms. Furthermore, systemic lupus erythematosus, neutrophil extracellular trap formation, and alcoholism were the primarily enriched KEGG pathways (Figure 2A).
FIGURE 2. Functional Enrichment Analysis. (A) Bar plots of KEGG pathway and GO enrichment analyses in the control and case groups. (B) Bar plots of KEGG pathway and GO enrichment analyses in the case and treatment groups. (C) Bar plots of KEGG pathway and GO enrichment analyses in the control and case groups. Each GO map presents the top 10 enriched terms in BP, CC, and MF. Each KEGG map shows the top 15 enriched KEGG pathways.
With respect to DE-mRNAs in the Treatment vs. Case group, 928 GO terms and 33 KEGG pathways were enriched by the 466 DE-mRNAs. The common enriched GO terms primarily included regulation of metal ion transport, anatomical structure homeostasis, negative regulation of neuron death, and others in BP; receptor complex, plasma membrane signaling receptor complex, integral component, and others in CC; and carboxylic acid binding, amide binding, kinase phosphatase activity, and others in MF. Additionally, hematopoietic cell lineage, neuroactive ligand-receptor interaction, transcriptional dysregulation in cancer, tryptophan metabolism, and others were the primarily enriched pathways (Figure 2B).
Furthermore, in the Treatment vs. Control group, the 476 DE-mRNAs enriched 952 GO terms. The major BP terms were chromosome segregation, nuclear division, sister chromatid segregation, nuclear chromosome segregation, and others. The major CC terms were chromosome, centromeric region, microtubule, chromosomal region, spindle, and others. The major MF terms were microtubule binding, receptor-ligand activity, tubulin binding, and others. Moreover, 31 KEGG pathways were enriched, and the top fifteen most significant pathways were shown in the chart, which primarily included alcoholism, neutrophil extracellular trap formation, and systemic lupus erythematosus, among others (Figure 2C).
3.4 ceRNA and ceRNA (circ) networks of the three groups
A ceRNA network containing 99 nodes and 127 edges and a ceRNA (circ) network containing 157 nodes and 354 edges were established using DE-mRNAs in the Case vs. Control group (Figure 3A, B). Based on the mRNAs in the Treatment vs. Case group, a ceRNA network with 104 nodes and 112 edges and a ceRNA (circ) network with 293 nodes and 492 edges were established (Supplementary Figures S4A, B). In the Treatment vs. Control group, a ceRNA network with 68 nodes and 66 edges and a ceRNA (circ) network with 226 nodes and 424 edges were established (Supplementary Figures S4C, D).
FIGURE 3. ceRNA network construction. (A) mRNA-miRNA-lncRNA network in control and case groups. (B) mRNA-miRNA-circRNA network in control and case groups. In the mRNA-miRNA-lncRNA (circRNA) network, the hexagon represents miRNA, the rhombus represents lncRNA (circRNA), the oval represents mRNA, and the color represents logFC value. The redder the color, the more upregulated the expression of the gene, and the bluer the color, the more downregulated the expression of the gene.
3.5 Functional enrichment analyses of mRNAs in the networks
Five hundred and sixty-seven GO terms and nine KEGG pathways were enriched by the mRNAs in the ceRNA network of the Case vs. Control group. For instance, the purine nucleoside metabolic process, nucleoside metabolic process, response to xenobiotic stimulus, leading-edge membrane, cell projection membrane, cell leading edge, carboxylic acid binding, metal ion transmembrane transporter activity, and organic acid:sodium symporter activity were the primary enriched terms. The nine KEGG pathways were basal cell carcinoma, biosynthesis of amino acids, GABAergic synapse, glutamatergic synapse, JAK-STAT signaling pathway, nicotinate and nicotinamide metabolism, nucleotide metabolism, purine metabolism, and synaptic vesicle cycle (Figure 4A). Additionally, the mRNAs in the ceRNA (circ) network of the Case vs. Control group enriched 694 GO terms and 17 KEGG pathways, such as response to xenobiotic stimulus, small molecule catabolic process, cell projection membrane, and carboxylic acid binding in GO. Parathyroid hormone synthesis, secretion, and action, glyoxylate and dicarboxylate metabolism, breast cancer, and basal cell carcinoma were the major KEGG pathways (Supplementary Figure S5A).
FIGURE 4. GO enrichment and KEGG pathway analyses of mRNAs enrolled in ceRNA networks. (A) KEGG pathway and GO enrichment analysis of mRNAs in the ceRNA (lncRNA) network in control and case groups. (B) KEGG pathway and GO enrichment analysis of mRNAs in the ceRNA (lncRNA) network in the case and treatment groups. (C) KEGG pathway and GO enrichment analysis of mRNAs in the ceRNA (lncRNA) network in the control and treatment groups. Each GO map shows the top 10 enriched phrases in BP, CC, and MF. On (A, B), the top 9 enriched KEGG pathways are shown, while (C) shows the top 15 enriched KEGG pathways.
With respect to the Case vs. Treatment group, the mRNAs in ceRNA network enriched axonogenesis, response to steroid hormone, epithelial cell development, dendrite membrane, G protein-coupled receptor binding, RNA polymerase II-specific DNA-binding transcription factor binding, and others among 426 GO terms. The nine enriched KEGG pathways primarily included GABAergic synapse, histidine metabolism, and mTOR signaling pathway (Figure 4B). In addition, the mRNAs in the ceRNA (circ) network commonly enriched 684 GO terms and 23 KEGG pathways, such as regulation of metal ion transport, glutamatergic synapse, carbohydrate binding, transcriptional dysregulation in cancer, and FoxO signaling pathway (Supplementary Figure S5B).
Finally, in the Treatment vs. Control group, the mRNAs from the ceRNA network enriched 314 GO terms and 14 KEGG pathways. The enriched BP terms were more significant than the CC and MF terms, such as mitotic cell cycle phase transition, cell cycle phase transition, and positive regulation of fibroblast proliferation. The most significantly enriched pathways included cell cycle, cellular senescence, and human papillomavirus infection, among others (Figure 4C). Besides, the enrichment results of mRNAs in the ceRNA (circ) network were similar in that the enriched BP terms were more significant than those in CC and MF as well, such as regulation of body fluid levels, G2/M transition of mitotic cell cycle, and cell cycle G2/M phase transition. Hematopoietic cell lineage, cytokine-cytokine receptor interaction, Staphylococcus aureus infection, and others were some of the primarily enriched KEGG pathways (Supplementary Figure S5C).
3.6 Identification and functional annotation of DE-ERSRGs
Thirteen intersected mRNAs (Fos, Hp, Serpina3n, Cdkn1a, Atf3, Nos2, Fkbp5, Tfrc, Nox4, Gabarapl1, Herpud1, Casq1, and Ppp1r15a) were identified from the overlap analysis; these were considered DE-ERSRGs (Supplementary Table S4; Figure 5A). The expression levels of DE-ERSRGs in the Case group were significantly different from those in the Control groups (Supplementary Figure S6A). Of note, the expression of Nos2, Tfrc, and Casq1 was higher in the control group (Figure 5B). Furthermore, these DE-ERSRGs were expressed differently between the Treatment and Case groups as well (Supplementary Figure S6B), and the expression levels of Tfrc, Nos2, and Casq1 were significantly higher in the treatment group (Figure 5C).
FIGURE 5. Identification of DE-ERSRGs and functional enrichement analysis. (A) Venn diagram of 13 DE-ERSRGs by overlapping DEGs in control vs. case group, DEGs in case group vs. treatment, with 826 ER stress-related genes. (B) Heat map of expression of 13 DE-ERSRGs in control and case groups. (C) Heat map of expression of 13 DE-ERSRGs in the case and treatment groups. (D) KEGG pathway analysis of 13 DE-ERSRGs. (E) GO enrichment analysis of 13 DE-ERSRGs.
Moreover, the enrichment results indicated that 24 KEGG and 429 GO term pathways were enriched by the 13 DE-ERSRGs. HIF-1 signaling was the most significantly enriched KEGG pathway, followed by Leishmaniasis, Pertussis, and others (Figure 5D). Furthermore, cellular response to extracellular stimulus, cellular response to external stimulus, response to extracellular stimulus, and other terms, were the primary BP terms. Meanwhile, the CC terms included dendrite cytoplasm, high-density lipoprotein particle, neuron projection cytoplasm, and others. Heat shock protein binding, ubiquitin protein ligase binding, kinase regulator activity, and others were the primary terms in MF (Figure 5E).
3.7 Establishment of DE-ERSRG-based ceRNA and ceRNA (circ) networks
The ceRNA network of 13 DE-ERSRGs illustrated that 15 nodes and 14 edges were present in the network, which only contained four of the 13 DE-ERSRGs (Cdkn1a, Atf3, Fkbp5, and Gabarapl1) (Figures 6A, B). In addition, a ceRNA (circ) network with 23 nodes and 28 edges was established based on the 13 DE-ERSRGs, which only contained the same four DE-ERSRGs (Figures 6C, D). Therefore, these four genes were regarded as key DE-ERSRGs. Similarly, according to the online dataset from GEO database (GSE210252), the expression patterns of four key DE-ERSRGs between mice received restraint stress or no stress were further explored and exhibited in Supplementary Figure S7, supporting that Cdkn1a, Atf3, and Gabarapl1 expressed higher in case group than that in control group, which were consistent to the results of our study above that shown in Supplementary Figure S6A.
FIGURE 6. Construction of ceRNA networks for 13 DE-ERSRGs. (A) mRNA-miRNA-lncRNA regulatory network. (B) Sankey diagram of the mRNA-miRNA-lncRNA regulatory. (C) mRNA-miRNA-circRNA regulatory network. (D) Sankey diagram of mRNA-miRNA-circRNA regulatory. In the mRNA-miRNA-lncRNA (circRNA) network, the hexagon represents miRNA, the rhombus represents lncRNA (circRNA), the oval represents mRNA, and the color represents logFC value. The closer the color is to red, the more upregulated the expression is, and the closer the color is to blue, the downregulated the expression is.
3.8 RT-qPCR validation of 13 DE-ERSRGs
To gain deeper insights on the expression levels of the 13 DE-ERSRGs (Fos, Hp, Serpina3n, Cdkn1a, Atf3, Nos2, Fkbp5, Tfrc, Nox4, Gabarapl1, Herpud1, Casq1, and Ppp1r15a), RT-qPCR was performed using eleven tissue samples. The expression levels of the DE-ERSRGs were distinctly different between the control and case samples (all p < .05), which was consistent with the results of previous bioinformatics analysis. Conversely, the expression levels in the case and treatment groups were significantly different, except for the levels of Ppp1r15a and Nos2. Although the expression levels of these two genes were not statistically significant, they were altered, as confirmed by the sequencing results (Figure 7).
FIGURE 7. Validation of the expression of 13 DE-ERSRGs by RT-qPCR. Validation of the expression levels of 13 DE-ERSRGs in mice tissue samples by quantitative RT-PCR. Error bars show mean ± SD. Data are representative of three independent experiments. *p < .05, **p < .01, ***p < .001.
4 Discussion
Numerous studies, including epidemiological and experimental studies in humans and animals, have yielded compelling evidence favoring the strong association between psychosocial factors and cardiovascular morbidity. In current European guidelines for the prevention of CVDs, psychological stress is recognized as a potential contributing factor to both CVD development and progression (Piepoli et al., 2016). Yet, the limited understanding of the molecular mechanisms underpinning chronic stress-induced cardiac injury has led to the lack of effective treatment options. Several studies have shown that psychological stress is a determinant of the prognosis and outcome of cardiovascular events, rather than a cause of CVD (Sara et al., 2021). On one hand, psychological stress can contribute to some adverse lifestyle patterns in humans, such as the abuse of tobacco and alcohol, poor eating habits, or sleep disorders. These confounding factors complicate the assessment of the effects of mental factors on the cardiovascular system in a given context. On the other hand, in animal experiments, subjects suffering from CVDs were simultaneously exposed to psychological stress to investigate the effect of stress on the progression of pre-existing CVDs, but limited attention was paid to the potential molecular processes induced by chronic psychological stress that induce cardiac tissue damage in healthy individuals. Cho et al. evaluated the molecular basis of cardiac damage in healthy persons exposed to psychological stress using a transcriptome sequencing approach, but their results failed to adequately account for the irreversible damage caused to the heart by chronic stress in American soldiers (Cho et al., 2014). This might be related to the methods (social defeat) they used, because the recurrence of the same stressor usually leads to adaptation and failure to reflect the unpredictable, uncontrollable, complex, and variable properties of naturally occurring stressors. In our study, we addressed and resolved some of the limitations of their study. The chronic psychological stress model was slightly modified from the chronic unpredictability stress model by increasing the stimulus intensity to reduce adaptation and better simulate chronic stress. The results suggested that mice that were subjected to the stress stimuli exhibited significant differences in all four behavioral experiments compared to control mice. In this modified model, we observed the abnormal expression of some ERSRGs, and the occurrence of ER stress may be a cause of permanent cardiac damage in healthy individuals suffering from chronic psychological stress.
In recent years, a large body of evidence from research has shown that ER stress is linked to the pathogenesis and progression of multiple CVDs and is an important target for the treatment of diverse CVDs. In this study, we identified 13 DE-ERSRGs. Among them, Fos was previously reported to be a downstream signaling molecule of ER stress in obesity models (Levi et al., 2018), and GRP78 was shown to interact with the transcription factor Fos to promote ER stress (Gopal and Pizzo, 2021). Fu et al. found that the overexpression of sarcoplasmic ER calcium ATPase (SERCA2a) reduced Hp expression while alleviating ER stress, thereby improving cardiac function in animals with heart failure. This suggests that Hp expression may be positively correlated with the ER stress levels in cardiac tissues (Fu et al., 2012). Serpina3 has been detected as a marker of ER stress in astrocytes in immune system diseases (Masvekar et al., 2022), and the role played by its Serpina3n (a protein from the same family) in ER stress has never been described. The role of NOX4, a member of the NADPH oxidase family, in ER stress in CVDs has garnered greater attention, as NOX4 participates in two UPR signaling pathways (Wu et al., 2010). In vascular smooth muscle cells and cardiomyocytes ontogenically related to VSMCs, high NOX4 expression induced by ER stress reduced cell death but promoted cell proliferation, leading to arterial remodeling and myocardial hypertrophy, and exacerbated the deterioration of cardiac function (Mittal et al., 2007; Kuroda et al., 2010). The PP1-Ppp1r15a phosphatase complex dephosphorylates eIF2α, and the anti-hypertensive compound guanabenz disrupts this complex and thereby attenuates ER stress in cardiomyocytes (Neuber et al., 2014). We also detected the upregulation of this complex in cardiac tissues post stress, which was reduced upon treatment with CB2R agonists. Herpud1 is an ER-resident membrane protein that plays a part in ER degradation, and its upregulation can exert a cytoprotective role, which is why we assumed that the upregulation of Herpud1 in this experiment is a compensatory protective mechanism (Navarro-Marquez et al., 2018). Furthermore, three DE-ERSRGs in our results showed a trend of downregulation in response to psychological stress and upregulation after treatment. Among them, Tfrc and Casq1 usually play a protective role and lead to the generation of ER stress when external factors cause their lowering (Jung et al., 2015; Hanna et al., 2021), which is also consistent with our results. Payne et al. suggested that Nos2 could protect against the development of ER stress (Payne et al., 2001), whereas in other studies, increased Nos2 expression has been suggested to exacerbate ER stress (Deslauriers et al., 2011; Zanotto et al., 2017). Whether the tendency of Nos2 downregulation in response to psychological stress results from a compensatory mechanism or some other mechanism remains elusive. Cdkn1a has been reported in ER stress (Lopez et al., 2015; Inoue et al., 2017). In this model, Cdkn1a probably acted as a negative regulator of the cell cycle, allowing cells sufficient time for repair of damaged DNA after external stimuli and various metabolic abnormalities. Moreover, Gabapal1 is involved in the formation of autophagosomal vacuoles (Chakrama et al., 2010). We hypothesized that it may participate in ER stress-induced autophagy. Atf3 is involved in ER stress through multiple pathways. In the most common pathway, Atf3 acts as an Atf4 downstream molecule in PERK-ATF4-ATF3-CHOP-induced apoptosis (Tang et al., 2020). In addition, Atf3 can also act as a transcription factor for NOX4 and upregulate NOX4-induced ROS production (Liu Z. et al., 2020). This connection can also explain the increased Nox4 expression observed in our results. Additionally, Fkbp5 is a well-known stress response gene that plays a central role in psychological stress-induced diseases (Le-Niculescu et al., 2020). The increase in Fkbp5 expression also corroborates the validity of our established model of chronic psychological stress. Targeted drugs of FKBP5 are being developed actively, and CB2R agonists may be one of the promising candidates.
Among the 13 DE-ERSRGs, four key gene could constitute the ceRNA regulative network. They have been reported to be strongly associated with the onset and progression of CVDs. The protein p21 encoded by cdkn1a is considered to be a marker of cardiomyocyte senescence in addition to inhibiting the cell cycle (Bernard et al., 2020). By increasing the levels of p21 intracellularly in cardiac myocytes, phenylalanine causes premature cardiomyocyte senescence (Czibik et al., 2021). Additionally, p21 overexpression in cardiac fibroblasts induces differentiation of primary cardiac fibroblasts into myofibroblasts and promotes cardiac remodeling (Roy et al., 2007). In a pressure overload-induced cardiac remodeling paradigm, cdkn1a knockdown mice showed reduced myocardial hypertrophy and fibrosis and improved cardiac function when compared to wild type mice (Xu and Tang, 2016). It is evident that cdkn1a is crucial in cardiac injury. Acute myocardial infarction (AMI) can induce the release of gabarapl1 from circulating endothelial cells, which can stimulate an increase in monocytes and neutrophils by recruiting NLRP3 inflammasomes, finally causing vascular inflammation during AMI (Zhang et al., 2022). ATF3, which functions as a “hub” for numerous adaptive cellular responses, is important for the progression and development of CVDs. ATF3 overexpression in cardiac tissues can promote myocardial hypertrophy and fibrosis (Okamoto et al., 2001). Elevated ATF3 expression due to different stimuli can damage endothelial cells to promote atherosclerosis (Cai et al., 2000). ATF3 was increased in carotid artery of hypertensive rats, and downregulation of ATF3 was associated with the protective effect of enalapril against hypertension (Liu et al., 2015). However, ATF3 plays a paradoxical role in the development of heart failure (Lin et al., 2014; Brooks et al., 2015). Thus, the role of ATF3 in CVDs is complicated. FKBP5 is a key molecule in the triggering of CVDs by psychological stress. FKBP5 can increase plaque destabilization by promoting local inflammation via the NF-κB pathway, according to previous research (Zannas et al., 2019). High levels of FKBP5 expression were also detected in platelets after AMI (Eicher et al., 2016). FKBP5 single nucleotide polymorphisms (SNPs) are closely related with depression and coronary heart disease co-morbidity (Wang et al., 2020).
Cho et al. assessed the dynamic course of gene expression in the mouse heart in the transcriptional level analysis after 1–10 days of exposure to social defeat (Cho et al., 2014). They identified extracellular matrix remodeling, immune responses (e.g., complement activation), and cell proliferation as the primary molecular-level alterations in cardiac tissues after psychological stress. Despite the differences in the models (social defeat vs. modified CMS), sequencing methods (microarray vs. whole transcriptome sequencing), and durations of stress (10 days vs. 3 weeks), we observed similar changes in gene expression patterns by analyzing the KEGG and GO enrichment results, suggesting that the same pathophysiological changes mentioned above were the basic reactions to psychological stress affecting the heart. In addition, we observed some notable changes, such as metabolic alterations in tissues (including glycine, serine, and threonine metabolism, purine metabolism, nucleotide metabolism, and nicotinate and nicotinamide metabolism.) and ER stress, which was the focus of this study.
Notably, the use of CB2R agonists diminishes the severity of stress in mice, as seen in the behavioral results, although there was no statistically significant difference between the results of the elevated plus maze and tail suspension tests, and only a trend toward improvement was observed. Possibly, ER stress occurs in neurons in regions such as the hippocampus, amygdala, and striatum, which are engaged in the development of psychological stress (Pavlovsky et al., 2013; Jangra et al., 2016; Jangra et al., 2017; Liu et al., 2019). CB2R may attenuate the symptoms of psychological stress by reducing the expression of ERSRGs and other pathways. CB2R agonists have been reported to exert neuroprotective effects to counteract the symptoms of psychological stress through anti-inflammatory functions (Zoppi et al., 2014).
However, there remain some limitations in our study. First, the sample size was small, which could have led to the low concordance between the qRT-PCR and sequencing results. Second, the experimental subjects used were male mice, which might have introduced gender bias. Third, the sequencing of whole heart tissues may have masked specific pathological changes in a particular region. We intend to follow up with experimental validation of the regulatory network of the four key genes identified and the role of CB2R agonists. We will also focus on the impact of metabolite changes in the results. In summary, our findings provide a comprehensive overview of the molecular changes in post-transcriptional regulation underlying chronic psychological stress-induced cardiac injury, with ER stress as one of the important mechanisms. CB2R agonists can alleviate ER stress through a competitive endogenous RNA mechanism and thus can be used as a therapeutic target in chronic psychological stress-induced cardiac injury.
5 Conclusion
Our findings revealed, for the first time, the expression profiles of circRNA, lncRNA, miRNA, and mRNA in cardiac tissues exposed to chronic psychological stress. We identified four key DE-ERSRGs-related genes and established relevant ceRNA networks. Our findings may provide novel insights on the molecular mechanisms underlying the effects of CB2R agonists in the treatment of chronic psychological stress-induced cardiac diseases.
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/geo/, GSE210820; https://www.ncbi.nlm.nih.gov/geo/, GSE210252.
Ethics statement
The animal study was reviewed and approved by Institutional Animal Care and Use Committee of Chinese PLA General Hospital.
Author contributions
FC and CQ conceived and designed the experiments. CQ, YZ (3rd author), and YZ (4th author) performed the experiments. CQ, YuW, and YaW analyzed the data. CQ wrote the manuscript. FC edited the manuscript.
Funding
This study was supported by the National Key Research and Development Projects (2022YFC3602400); the Key Health Care Projects of National Health Commission (2020ZD05); Basic Research Reinforcement Project (2022-JCJQ-ZD-079-00); Project of National Clinical Research Center for Geriatric Diseases (NCRCG-PLAGH-2022007).
Acknowledgments
We thank Prof. Jiyan Zhang at Beijing Institute of Basic Medical Sciences for providing the behavioral equipment.
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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.1095428/full#supplementary-material
References
Almanza, A., Carlesso, A., Chintha, C., Creedican, S., Doultsinos, D., Leuzzi, B., et al. (2019). Endoplasmic reticulum stress signalling - from basic mechanisms to clinical applications. FEBS J. 286 (2), 241–278. doi:10.1111/febs.14608
Atkins, C., Liu, Q., Minthorn, E., Zhang, S. Y., Figueroa, D. J., Moss, K., et al. (2013). Characterization of a novel PERK kinase inhibitor with antitumor and antiangiogenic activity. Cancer Res. 73 (6), 1993–2002. doi:10.1158/0008-5472.CAN-12-3109
Bernard, M., Yang, B., Migneault, F., Turgeon, J., Dieudé, M., Olivier, M. A., et al. (2020). Autophagy drives fibroblast senescence through MTORC2 regulation. Autophagy 16 (11), 2004–2016. doi:10.1080/15548627.2020.1713640
Bertolotti, A., Zhang, Y., Hendershot, L. M., Harding, H. P., and Ron, D. (2000). Dynamic interaction of BiP and ER stress transducers in the unfolded-protein response. Nat. Cell Biol. 2 (6), 326–332. doi:10.1038/35014014
Bi, X., Zhang, G., Wang, X., Nguyen, C., May, H. I., Li, X., et al. (2018). Endoplasmic reticulum chaperone GRP78 protects heart from ischemia/reperfusion injury through akt activation. Circ. Res. 122 (11), 1545–1554. doi:10.1161/CIRCRESAHA.117.312641
Blackwood, E. A., Hofmann, C., Santo Domingo, M., Bilal, A. S., Sarakki, A., Stauffer, W., et al. (2019). ATF6 regulates cardiac hypertrophy by transcriptional induction of the mTORC1 activator, rheb. Circ. Res. 124 (1), 79–93. doi:10.1161/CIRCRESAHA.118.313854
Brooks, A. C., DeMartino, A. M., Brainard, R. E., Brittian, K. R., Bhatnagar, A., and Jones, S. P. (2015). Induction of activating transcription factor 3 limits survival following infarct-induced heart failure in mice. Am. J. Physiol. Heart Circ. Physiol. 309 (8), H1326–H1335. doi:10.1152/ajpheart.00513.2015
Cai, Y., Zhang, C., Nawa, T., Aso, T., Tanaka, M., Oshiro, S., et al. (2000). Homocysteine-responsive ATF3 gene expression in human vascular endothelial cells: Activation of c-jun NH(2)-terminal kinase and promoter response element. Blood 96 (6), 2140–2148. doi:10.1182/blood.v96.6.2140.h8002140_2140_2148
Chakrama, F. Z., Seguin-Py, S., Le Grand, J. N., Fraichard, A., Delage-Mourroux, R., Despouy, G., et al. (2010). GABARAPL1 (GEC1) associates with autophagic vesicles. Autophagy 6 (4), 495–505. doi:10.4161/auto.6.4.11819
Chen, M. H., Pan, T. L., Li, C. T., Lin, W. C., Chen, Y. S., Lee, Y. C., et al. (2015). Risk of stroke among patients with post-traumatic stress disorder: Nationwide longitudinal study. Br. J. Psychiatry 206 (4), 302–307. doi:10.1192/bjp.bp.113.143610
Cho, J. H., Lee, I., Hammamieh, R., Wang, K., Baxter, D., Scherler, K., et al. (2014). Molecular evidence of stress-induced acute heart injury in a mouse model simulating posttraumatic stress disorder. Proc. Natl. Acad. Sci. U. S. A. 111 (8), 3188–3193. doi:10.1073/pnas.1400113111
Cho, W. H., Noh, K., Lee, B. H., Barcelon, E., Jun, S. B., Park, H. Y., et al. (2022). Hippocampal astrocytes modulate anxiety-like behavior. Nat. Commun. 13 (1), 6536. doi:10.1038/s41467-022-34201-z
Czibik, G., Mezdari, Z., Murat Altintas, D., Bréhat, J., Pini, M., d'Humières, T., et al. (2021). Dysregulated phenylalanine catabolism plays a key role in the trajectory of cardiac aging. Circulation 144 (7), 559–574. doi:10.1161/circulationaha.121.054204
Deslauriers, A. M., Afkhami-Goli, A., Paul, A. M., Bhat, R. K., Acharjee, S., Ellestad, K. K., et al. (2011). Neuroinflammation and endoplasmic reticulum stress are coregulated by crocin to prevent demyelination and neurodegeneration. J. Immunol. 187 (9), 4788–4799. doi:10.4049/jimmunol.1004111
Ebrahimi, R., Lynch, K. E., Beckham, J. C., Dennis, P. A., Viernes, B., Tseng, C. H., et al. (2021). Association of posttraumatic stress disorder and incident ischemic heart disease in women veterans. JAMA Cardiol. 6 (6), 642–651. doi:10.1001/jamacardio.2021.0227
Eicher, J. D., Wakabayashi, Y., Vitseva, O., Esa, N., Yang, Y., Zhu, J., et al. (2016). Characterization of the platelet transcriptome by RNA sequencing in patients with acute myocardial infarction. Platelets 27 (3), 230–239. doi:10.3109/09537104.2015.1083543
Frakes, A. E., and Dillin, A. (2017). The UPR(ER): Sensor and coordinator of organismal homeostasis. Mol. Cell 66 (6), 761–771. doi:10.1016/j.molcel.2017.05.031
Fu, Z., Li, F., Jia, L., Su, S., Wang, Y., Cai, Z., et al. (2019). Histone deacetylase 6 reduction promotes aortic valve calcification via an endoplasmic reticulum stress-mediated osteogenic pathway. J. Thorac. Cardiovasc Surg. 158 (2), 408–417. doi:10.1016/j.jtcvs.2018.10.136
Fu, Z. Q., Li, X. Y., Lu, X. C., Mi, Y. F., Liu, T., and Ye, W. H. (2012). Ameliorated stress related proteins are associated with improved cardiac function by sarcoplasmic reticulum calcium ATPase gene transfer in heart failure. J. Geriatr. Cardiol. 9 (3), 269–277. doi:10.3724/SP.J.1263.2012.05299
Glaser, R., and Kiecolt-Glaser, J. K. (2005). Stress-induced immune dysfunction: Implications for health. Nat. Rev. Immunol. 5 (3), 243–251. doi:10.1038/nri1571
Gold, P. W., Licinio, J., and Pavlatou, M. G. (2013). Pathological parainflammation and endoplasmic reticulum stress in depression: Potential translational targets through the CNS insulin, klotho and PPAR-gamma systems. Mol. Psychiatry 18 (2), 154–165. doi:10.1038/mp.2012.167
Gold, S. M., Kohler-Forsberg, O., Moss-Morris, R., Mehnert, A., Miranda, J. J., Bullinger, M., et al. (2020). Comorbid depression in medical diseases. Nat. Rev. Dis. Prim. 6 (1), 69. doi:10.1038/s41572-020-0200-2
Gopal, U., and Pizzo, S. V. (2021). Cell surface GRP78 signaling: An emerging role as a transcriptional modulator in cancer. J. Cell Physiol. 236 (4), 2352–2363. doi:10.1002/jcp.30030
Gu, Y., Ye, T., Tan, P., Tong, L., Ji, J., Gu, Y., et al. (2021). Tolerance-inducing effect and properties of innate immune stimulation on chronic stress-induced behavioral abnormalities in mice. Brain Behav. Immun. 91, 451–471. doi:10.1016/j.bbi.2020.11.002
Hanna, A. D., Lee, C. S., Babcock, L., Wang, H., Recio, J., and Hamilton, S. L. (2021). Pathological mechanisms of vacuolar aggregate myopathy arising from a Casq1 mutation. FASEB J. 35 (5), e21349. doi:10.1096/fj.202001653RR
Hinterdobler, J., Schott, S., Jin, H., Meesmann, A., Steinsiek, A. L., Zimmermann, A. S., et al. (2021). Acute mental stress drives vascular inflammation and promotes plaque destabilization in mouse atherosclerosis. Eur. Heart J. 42 (39), 4077–4088. doi:10.1093/eurheartj/ehab371
Hourani, L. L., Davila, M. I., Morgan, J., Meleth, S., Ramirez, D., Lewis, G., et al. (2020). Mental health, stress, and resilience correlates of heart rate variability among military reservists, guardsmen, and first responders. Physiol. Behav. 214, 112734. doi:10.1016/j.physbeh.2019.112734
Inoue, Y., Kawachi, S., Ohkubo, T., Nagasaka, M., Ito, S., Fukuura, K., et al. (2017). The CDK inhibitor p21 is a novel target gene of ATF4 and contributes to cell survival under ER stress. FEBS Lett. 591 (21), 3682–3691. doi:10.1002/1873-3468.12869
Jangra, A., Dwivedi, S., Sriram, C. S., Gurjar, S. S., Kwatra, M., Sulakhiya, K., et al. (2016). Honokiol abrogates chronic restraint stress-induced cognitive impairment and depressive-like behaviour by blocking endoplasmic reticulum stress in the hippocampus of mice. Eur. J. Pharmacol. 770, 25–32. doi:10.1016/j.ejphar.2015.11.047
Jangra, A., Sriram, C. S., Dwivedi, S., Gurjar, S. S., Hussain, M. I., Borah, P., et al. (2017). Sodium phenylbutyrate and edaravone abrogate chronic restraint stress-induced behavioral deficits: Implication of oxido-nitrosative, endoplasmic reticulum stress cascade, and neuroinflammation. Cell Mol. Neurobiol. 37 (1), 65–81. doi:10.1007/s10571-016-0344-5
Jordan, H. T., Stellman, S. D., Morabia, A., Miller-Archie, S. A., Alper, H., Laskaris, Z., et al. (2013). Cardiovascular disease hospitalizations in relation to exposure to the September 11, 2001 World Trade Center disaster and posttraumatic stress disorder. J. Am. Heart Assoc. 2 (5), e000431. doi:10.1161/JAHA.113.000431
Jung, I. R., Choi, S. E., Jung, J. G., Lee, S. A., Han, S. J., Kim, H. J., et al. (2015). Involvement of iron depletion in palmitate-induced lipotoxicity of beta cells. Mol. Cell Endocrinol. 407, 74–84. doi:10.1016/j.mce.2015.03.007
Kuroda, J., Ago, T., Matsushima, S., Zhai, P., Schneider, M. D., and Sadoshima, J. (2010). NADPH oxidase 4 (Nox4) is a major source of oxidative stress in the failing heart. Proc. Natl. Acad. Sci. U. S. A. 107 (35), 15565–15570. doi:10.1073/pnas.1002178107
Le-Niculescu, H., Roseberry, K., Levey, D. F., Rogers, J., Kosary, K., Prabha, S., et al. (2020). Towards precision medicine for stress disorders: Diagnostic biomarkers and targeted drugs. Mol. Psychiatry 25 (5), 918–938. doi:10.1038/s41380-019-0370-z
Leor, J., Poole, W. K., and Kloner, R. A. (1996). Sudden cardiac death triggered by an earthquake. N. Engl. J. Med. 334 (7), 413–419. doi:10.1056/NEJM199602153340701
Levi, N. J., Wilson, C. W., Redweik, G. A. J., Gray, N. W., Grzybowski, C. W., Lenkey, J. A., et al. (2018). Obesity-related cellular stressors regulate gonadotropin releasing hormone gene expression via c-Fos/AP-1. Mol. Cell Endocrinol. 478, 97–105. doi:10.1016/j.mce.2018.07.011
Li, X., Han, D., Tian, Z., Gao, B., Fan, M., Li, C., et al. (2016). Activation of cannabinoid receptor type II by AM1241 ameliorates myocardial fibrosis via nrf2-mediated inhibition of TGF-β1/smad3 pathway in myocardial infarction mice. Cell Physiol. Biochem. 39 (4), 1521–1536. doi:10.1159/000447855
Lin, H., Li, H. F., Chen, H. H., Lai, P. F., Juan, S. H., Chen, J. J., et al. (2014). Activating transcription factor 3 protects against pressure-overload heart failure via the autophagy molecule Beclin-1 pathway. Mol. Pharmacol. 85 (5), 682–691. doi:10.1124/mol.113.090092
Liu, L., Liu, J., Huang, Z., Yu, X., Zhang, X., Dou, D., et al. (2015). Berberine improves endothelial function by inhibiting endoplasmic reticulum stress in the carotid arteries of spontaneously hypertensive rats. Biochem. Biophys. Res. Commun. 458 (4), 796–801. doi:10.1016/j.bbrc.2015.02.028
Liu, L., Zhao, Z., Lu, L., Liu, J., Sun, J., Wu, X., et al. (2019). Icariin and icaritin ameliorated hippocampus neuroinflammation via inhibiting HMGB1-related pro-inflammatory signals in lipopolysaccharide-induced inflammation model in C57BL/6J mice. Int. Immunopharmacol. 68, 95–105. doi:10.1016/j.intimp.2018.12.055
Liu, X., Zhang, D., Dong, X., Zhu, R., Ye, Y., Li, L., et al. (2020a). Pharmacological activation of CB2 receptor protects against ethanol-induced myocardial injury related to RIP1/RIP3/MLKL-mediated necroptosis. Mol. Cell Biochem. 474 (1-2), 1–14. doi:10.1007/s11010-020-03828-1
Liu, Z., Gu, S., Lu, T., Wu, K., Li, L., Dong, C., et al. (2020b). IFI6 depletion inhibits esophageal squamous cell carcinoma progression through reactive oxygen species accumulation via mitochondrial dysfunction and endoplasmic reticulum stress. J. Exp. Clin. Cancer Res. 39 (1), 144. doi:10.1186/s13046-020-01646-3
Lopez, I., Tournillon, A. S., Nylander, K., and Fahraeus, R. (2015). p53-mediated control of gene expression via mRNA translation during Endoplasmic Reticulum stress. Cell Cycle 14 (21), 3373–3378. doi:10.1080/15384101.2015.1090066
Mao, J., Hu, Y., Ruan, L., Ji, Y., and Lou, Z. (2019). Role of endoplasmic reticulum stress in depression (Review). Mol. Med. Rep. 20 (6), 4774–4780. doi:10.3892/mmr.2019.10789
Masvekar, R., Kosa, P., Barbour, C., Milstein, J. L., and Bielekova, B. (2022). Drug library screen identifies inhibitors of toxic astrogliosis. Mult. Scler. Relat. Disord. 58, 103499. doi:10.1016/j.msard.2022.103499
Meisel, S. R., Kutz, I., Dayan, K. I., Pauzner, H., Chetboun, I., Arbel, Y., et al. (1991). Effect of Iraqi missile war on incidence of acute myocardial infarction and sudden death in Israeli civilians. Lancet 338 (8768), 660–661. doi:10.1016/0140-6736(91)91234-l
Mittal, M., Roth, M., Konig, P., Hofmann, S., Dony, E., Goyal, P., et al. (2007). Hypoxia-dependent regulation of nonphagocytic NADPH oxidase subunit NOX4 in the pulmonary vasculature. Circ. Res. 101 (3), 258–267. doi:10.1161/CIRCRESAHA.107.148015
Molica, F., Matter, C. M., Burger, F., Pelli, G., Lenglet, S., Zimmer, A., et al. (2012). Cannabinoid receptor CB2 protects against balloon-induced neointima formation. Am. J. Physiol. Heart Circ. Physiol. 302 (5), H1064–H1074. doi:10.1152/ajpheart.00444.2011
Mondal, T. K., Emeny, R. T., Gao, D., Ault, J. G., Kasten-Jolly, J., and Lawrence, D. A. (2015). A physical/psychological and biological stress combine to enhance endoplasmic reticulum stress. Toxicol. Appl. Pharmacol. 289 (2), 313–322. doi:10.1016/j.taap.2015.09.013
Navarro-Marquez, M., Torrealba, N., Troncoso, R., Vasquez-Trincado, C., Rodriguez, M., Morales, P. E., et al. (2018). Herpud1 impacts insulin-dependent glucose uptake in skeletal muscle cells by controlling the Ca(2+)-calcineurin-Akt axis. Biochim. Biophys. Acta Mol. Basis Dis. 1864 (5), 1653–1662. doi:10.1016/j.bbadis.2018.02.018
Neuber, C., Uebeler, J., Schulze, T., Sotoud, H., El-Armouche, A., and Eschenhagen, T. (2014). Guanabenz interferes with ER stress and exerts protective effects in cardiac myocytes. PLoS One 9 (6), e98893. doi:10.1371/journal.pone.0098893
Ochoa, C. D., Wu, R. F., and Terada, L. S. (2018). ROS signaling and ER stress in cardiovascular disease. Mol. Asp. Med. 63, 18–29. doi:10.1016/j.mam.2018.03.002
Okamoto, Y., Chaves, A., Chen, J., Kelley, R., Jones, K., Weed, H. G., et al. (2001). Transgenic mice with cardiac-specific expression of activating transcription factor 3, a stress-inducible gene, have conduction abnormalities and contractile dysfunction. Am. J. Pathol. 159 (2), 639–650. doi:10.1016/s0002-9440(10)61735-x
Pavlovsky, A. A., Boehning, D., Li, D., Zhang, Y., Fan, X., and Green, T. A. (2013). Psychological stress, cocaine and natural reward each induce endoplasmic reticulum stress genes in rat brain. Neuroscience 246, 160–169. doi:10.1016/j.neuroscience.2013.04.057
Payne, C. M., Bernstein, H., Bernstein, C., Kunke, K., and Garewal, H. (2001). The specific NOS2 inhibitor, 1400W, sensitizes HepG2 cells to genotoxic, oxidative, xenobiotic, and endoplasmic reticulum stresses. Antioxid. Redox Signal 3 (5), 931–936. doi:10.1089/15230860152665082
Piepoli, M. F., Hoes, A. W., Agewall, S., Albus, C., Brotons, C., Catapano, A. L., et al. (2016). 2016 European Guidelines on cardiovascular disease prevention in clinical practice. Rev. Esp. Cardiol. Engl. Ed. 69 (10), 939. doi:10.1016/j.rec.2016.09.009
Ron, D., and Walter, P. (2007). Signal integration in the endoplasmic reticulum unfolded protein response. Nat. Rev. Mol. Cell Biol. 8 (7), 519–529. doi:10.1038/nrm2199
Roy, S., Khanna, S., Rink, T., Radtke, J., Williams, W. T., Biswas, S., et al. (2007). P21waf1/cip1/sdi1 as a central regulator of inducible smooth muscle actin expression and differentiation of cardiac fibroblasts to myofibroblasts. Mol. Biol. Cell 18 (12), 4837–4846. doi:10.1091/mbc.e07-03-0270
Roy, S. S., Foraker, R. E., Girton, R. A., and Mansfield, A. J. (2015). Posttraumatic stress disorder and incident heart failure among a community-based sample of US veterans. Am. J. Public Health 105 (4), 757–763. doi:10.2105/AJPH.2014.302342
Sara, J. D. S., Lerman, L. O., and Lerman, A. (2021). The endothelium is a key player in the vascular response to acute mental stress. Eur. Heart J. 42 (39), 4089–4091. doi:10.1093/eurheartj/ehab510
Shah, A. J., Weeks, V., Lampert, R., Bremner, J. D., Kutner, M., Raggi, P., et al. (2022). Early life trauma is associated with increased microvolt T-wave alternans during mental stress challenge: A substudy of mental stress ischemia: Prognosis and genetic influences. J. Am. Heart Assoc. 11 (5), e021582. doi:10.1161/JAHA.121.021582
Steinberg, J. S., Arshad, A., Kowalski, M., Kukar, A., Suma, V., Vloka, M., et al. (2004). Increased incidence of life-threatening ventricular arrhythmias in implantable defibrillator patients after the World Trade Center attack. J. Am. Coll. Cardiol. 44 (6), 1261–1264. doi:10.1016/j.jacc.2004.06.032
Tang, Q., Ren, L., Liu, J., Li, W., Zheng, X., Wang, J., et al. (2020). Withaferin A triggers G2/M arrest and intrinsic apoptosis in glioblastoma cells via ATF4-ATF3-CHOP axis. Cell Prolif. 53 (1), e12706. doi:10.1111/cpr.12706
Tang, V., Fu, S., Rayner, B. S., and Hawkins, C. L. (2019). 8-Chloroadenosine induces apoptosis in human coronary artery endothelial cells through the activation of the unfolded protein response. Redox Biol. 26, 101274. doi:10.1016/j.redox.2019.101274
Vaccarino, V., Almuwaqqat, Z., Kim, J. H., Hammadah, M., Shah, A. J., Ko, Y. A., et al. (2021). Association of mental stress-induced myocardial ischemia with cardiovascular events in patients with coronary heart disease. JAMA 326 (18), 1818–1828. doi:10.1001/jama.2021.17649
Wang, H., Wang, C., Song, X., Liu, H., Zhang, Y., and Jiang, P. (2020). Association of FKBP5 polymorphisms with patient susceptibility to coronary artery disease comorbid with depression. PeerJ 8, e9286. doi:10.7717/peerj.9286
Wang, M., and Kaufman, R. J. (2016). Protein misfolding in the endoplasmic reticulum as a conduit to human disease. Nature 529 (7586), 326–335. doi:10.1038/nature17041
Wu, R. F., Ma, Z., Liu, Z., and Terada, L. S. (2010). Nox4-derived H2O2 mediates endoplasmic reticulum signaling through local Ras activation. Mol. Cell Biol. 30 (14), 3553–3568. doi:10.1128/MCB.01445-09
Xu, M., and Tang, Q. (2016). GW27-e0439 Cdkn1α knockout alleviates pressure overload induced cardiac remodeling. J. Am. Coll. Cardiol., 68(16), C58–C58. doi:10.1016/j.jacc.2016.07.212
Yan, J. J., Ding, X. J., He, T., Chen, A. X., Zhang, W., Yu, Z. X., et al. (2022). A circuit from the ventral subiculum to anterior hypothalamic nucleus GABAergic neurons essential for anxiety-like behavioral avoidance. Nat. Commun. 13 (1), 7464. doi:10.1038/s41467-022-35211-7
Yu, W., Jin, G., Zhang, J., and Wei, W. (2019). Selective activation of cannabinoid receptor 2 attenuates myocardial infarction via suppressing NLRP3 inflammasome. Inflammation 42 (3), 904–914. doi:10.1007/s10753-018-0945-x
Zannas, A. S., Jia, M., Hafner, K., Baumert, J., Wiechmann, T., Pape, J. C., et al. (2019). Epigenetic upregulation of FKBP5 by aging and stress contributes to NF-κB-driven inflammation and cardiovascular risk. Proc. Natl. Acad. Sci. U. S. A. 116 (23), 11370–11379. doi:10.1073/pnas.1816847116
Zanotto, T. M., Quaresma, P. G. F., Guadagnini, D., Weissmann, L., Santos, A. C., Vecina, J. F., et al. (2017). Blocking iNOS and endoplasmic reticulum stress synergistically improves insulin resistance in mice. Mol. Metab. 6 (2), 206–218. doi:10.1016/j.molmet.2016.12.005
Zhang, B. F., Jiang, H., Chen, J., Guo, X., Li, Y., Hu, Q., et al. (2019). Nobiletin ameliorates myocardial ischemia and reperfusion injury by attenuating endoplasmic reticulum stress-associated apoptosis through regulation of the PI3K/AKT signal pathway. Int. Immunopharmacol. 73, 98–107. doi:10.1016/j.intimp.2019.04.060
Zhang, Q., Guan, G., Cheng, P., Cheng, W., Yang, L., and Wu, A. (2021). Characterization of an endoplasmic reticulum stress-related signature to evaluate immune features and predict prognosis in glioma. J. Cell Mol. Med. 25 (8), 3870–3884. doi:10.1111/jcmm.16321
Zhang, T., Hou, D., He, J., Zeng, X., Liu, R., Liu, L., et al. (2022). Oxidative-damaged mitochondria activate GABARAPL1-induced NLRP3 inflammasomes in an autophagic-exosome manner after acute myocardial ischemia. Oxid. Med. Cell Longev. 2022, 7958542. doi:10.1155/2022/7958542
Keywords: chronic psychological stress, endoplasmic reticulum stress, cardiac injuries, CB2R, ceRNA network, functional enrichment analysis
Citation: Qin C, Wang Y, Zhang Y, Zhu Y, Wang Y and Cao F (2023) Transcriptome-wide analysis reveals the molecular mechanisms of cannabinoid type II receptor agonists in cardiac injury induced by chronic psychological stress. Front. Genet. 13:1095428. doi: 10.3389/fgene.2022.1095428
Received: 11 November 2022; Accepted: 27 December 2022;
Published: 10 January 2023.
Edited by:
Bin Liu, Jiangsu Ocean Universiity, ChinaReviewed by:
Faming Zhao, Huazhong University of Science and Technology, ChinaDong Zhang, Qilu Hospital, Shandong University, China
Copyright © 2023 Qin, Wang, Zhang, Zhu, Wang and Cao. 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: Feng Cao, fengcao8828@163.com