A systems approach identifies co-signaling molecules of early growth response 1 transcription factor in immobilization stress
BMC Systems Biology volume 8, Article number: 100 (2014)
Adaptation to stress is critical for survival. The adrenal medulla, the major source of epinephrine, plays an important role in the development of the hyperadenergic state and increased risk for stress associated disorders, such as hypertension and myocardial infarction. The transcription factor Egr1 plays a central role in acute and repeated stress, however the complexity of the response suggests that other transcription factor pathways might be playing equally important roles during acute and repeated stress. Therefore, we sought to discover such factors by applying a systems approach.
Using microarrays and network analysis we show here for the first time that the transcription factor signal transducer and activator of transcription 3 (Stat3) gene is activated in acute stress whereas the prolactin releasing hormone (Prlh11) and chromogranin B (Chgb) genes are induced in repeated immobilization stress and that along with Egr1 may be critical mediators of the stress response.
Our results suggest possible involvement of Stat3 and Prlh1/Chgb up-regulation in the transition from short to repeated stress activation.
The adrenal medulla plays a key role in the response to acute and chronic stress. It is the major site of biosynthesis of epinephrine (Epi) in the periphery. Upon exposure to stress, the release of adrenomedullary Epi and norepinephrine (NE) are among the most rapid response to handle the emergency situation. This is crucial for activation of the “fight or flight” response to deal with a threat to homeostasis.
When stress is prolonged or repeated the adrenal medulla exhibits crucial adaptive and subsequently maladaptive responses. These include important changes in gene expression. The best characterized are the up-regulation of expression of catecholamine biosynthetic enzymes (reviewed in ). Exposure to single immobilization stress triggers a manifold elevation in transcription and expression of mRNAs of the catecholamine biosynthetic enzymes, tyrosine hydroxylase (TH), the first and major rate limiting enzyme in catecholamine biosynthesis as well as of phenylethanolamine N-methyltransferase (PNMT), the enzyme which catalyses the conversion of NE to Epi −. This rise is transient and returns to normal within one day and is not sufficient for substantial increase in their activity. Following repeated IMO stress the increase in gene expression of catecholamine biosynthetic enzymes is now more sustained with prolonged elevation in catecholamine biosynthetic enzyme activity (reviewed in ).
To determine the repertoire of changes and the mechanism of alterations in gene expression mediating the response of adrenal medulla to single and repeated stress, microarray profiling was performed . Following single exposure to 2 hr IMO stress there was altered expression, of greater than 2 fold, in nearly 4% of the total transcripts. Transcription factors and cell signaling genes displayed the most prevalent changes. Approximately 20% of the transcripts up-regulated by single IMO were transcription factors. Not only was Egr1 mRNA markedly induced in the adrenal medulla by single as well as repeated exposure to IMO, but pathway analysis indicated that Egr1 likely plays a central role .
Egr1 (Zif268, NGFI-A, TIS8 or Krox24) is a transcription factor with three zinc fingers of the Cys2His2 class (reviewed by ,). Egr1 binds to a GC-rich motif (5’-GCG (T/G) GGGCG-3’) through its three zinc finger DNA binding domains  and modulates transcription of a number of genes that participate in various cellular functions (reviewed by ,). Egr1 plays critical roles in divergent cellular processes. For example, Egr1 and Stat3 have been implicated in neuronal differentiation, specifically during neurite outgrowth (reviewed in ,, in tumor development -, oxidant stress , immune responses  and in insulin signaling and in nutrition ).
Egr1 target genes include catecholamine biosynthetic enzymes. Transcription of both TH and PNMT is up-regulated by Egr1 −. We have previously shown that Egr1 is markedly induced in the adrenal medulla by IMO stress . While barely expressed under basal conditions, immunofluorescence demonstrated widespread expression in the nucleus of TH expressing chromaffin cells in the adrenal medulla after IMO stress . However the molecules that form the core of the signaling cascade inducing these responses are not well understood. Because complex biological behaviors arise from the coordinated behavior of sets of genes acting in concert (gene modules), we hypothesized that genes that are co-expressed with Egr1 during single or repeated IMO stress might provide insights into to significant signaling pathways that participate in stress signaling. Here we employed Gene Set Enrichment Analysis to identify Egr1 co-expressed genes from IMO microarrays, extracted their interactors and all their interrelationships and reconstructed Egr1 networks. From their network properties, we have identified the transcription factor Stat3 and the peptide Prlh1 in short and prolonged stress respectively as Egr1 neighbors in the adrenal medulla implicating them for the first time in stress signaling.
Gene sets that enrich with Egr1 expression in acute and repeated stress
Acute and repeated stress responses are accompanied by different patterns of gene expression, particularly of transcription factor genes, suggesting an interplay of transcription factors and the gene programs they control. In order to identify novel genes and their products that might be instrumental in networks leading from acute to repeated stress, we applied a strategy (Figure 1) that allowed us the re-construction of Egr1-centered networks and the extraction of network neighbors from expression profiles identical to that of Egr1. We executed this strategy in two steps: First, we used Kolmogorov-Smirnov analysis of acute (1 ×) and repeated (6×) IMO microarray expression data in order to rank expression levels of all genes in the microarrays. Second we used gene set enrichment analysis (GSEA) and computed gene module enrichment scores (ES) for each module of genes that are either co-expressed with Egr1 (positive ES) or anti-coexpressed (negative ES) (Additional file 1: Figure S1). Specifically, using Egr1 as an index gene in GSEA we extracted the top fifty genes (Egr1_POS module, top positive ES score, control vs. 1 ×) that are co-expression neighbors of Egr1 (See Additional file 2, Computational and Bioinformatic Methods, for a full account of methods). Second, by categorical class analysis (control vs. 1×) we extracted the top fifty down-regulated genes and top fifty up-regulated genes in control and in 1× samples. We repeated the same analyses for 6× IMO data. Third, we extracted the top negatively correlated genes from the top gene sets (Egr1_NEG module, top negative ES scores). Finally, we mined a list of 300 rat genes from the COXPRESdb database reported to be co-expressed with Egr1 (Data are available online as Additional files 3 and 4).
Extraction of all PPIs/interrelationships and Reconstruction of Egr1-centered networks
In the next step of our strategy we combined the three previously discussed gene/protein lists and extracted all genetic and physical interactions and then we text-mined all published interrelationships from public databases with statistical tools that are incorporated within the websites, using an expectation value (E value) greater than 0.7. We verified that data reflected real interrelationships and rejected data that were mere co-incidences of textual referral. Thus, interaction data with an E value less than 0.7 were rejected (see Additional file 5).
Re-construction of 1× and 6× IMO stress networks
In the final step of our strategy we re-constructed 1× and 6× networks as described (Additional file 2: Computational and Bioinformatic Methods). The 1× network contains 1717 genes/proteins (nodes) and 6554 interactions (edges) and is whereas the 6× contains 1313 nodes and 5203 edges. It can be seen that the network parameters of 1× and 6× networks are similar to the parameters of other biological networks such as the yeast proteome or the human HTFN however, they are quite dissimilar to random ER networks (Table 1). We then re-displayed networks 1× and 6× around Egr1 and extracted the Egr1 network neighborhoods (Figure 2, left upper and right panels respectively, and Figure 3, left panel) within the HUBBA website by calculating the intersection between Egr1 and top hub and bottleneck nodes (see Additional file 6, “top hubs” Excel sheet). Notably, the 1× network is organized around Egr1, a result consistent with the central role of this transcription factor in acute stress, and also around Stat3, which is a novel observation (Figure 2, lower panel). In contrast, while Egr1 remains, as expected, a top hub in the 6× network, the other top hubs in this network, Apoa1 and Hspd1 are not neighbors of Egr1, suggesting that the Egr1 neighborhood is organized differently (Figure 3, right panel). In order to further analyze Egr1’s network neighbors in 1× and 6× networks, and infer possible links between them and Egr1 we extracted the top 10 motifs of interacting proteins with the MCODE algorithm within Cytoscape and identified top GO functional classifications and KEGG pathways with the GATHER algorithm for nodes in the neighborhood of Egr1. The Egr1-centered, 1× IMO network neighborhood motifs are enriched for cell cycle, MAPK kinase and cytokine pathways, consistent with their role in acute stress (Additional file 6 and Table 2), whereas the 6× Egr1-centered network motifs is enriched for purine metabolism, insulin signaling, MAPK and the neuroactive signaling pathways (Additional file 7 and Table 3).
We further narrowed our focus on the shortest path neighbors of Egr1 within the top motifs in both acute and repeated stress. Egr1 and Stat3 also are top bottlenecks as well as network neighbors in the 1× network and they are in motif 4 (Additional file 6). Calculation of shortest paths in the neighborhood of Egr1within the HUBBA website revealed that Stat3 is a network neighbor of Egr1 (Figure 2, left upper panel) and Additional file 6, see “top hubs” sheet). Also, Egr1 is a top hub in 6× IMO samples, and is a member of motif 4 along with Prlh and Chgb (Figure 3, right panel). Few of the other nodes in other motifs were of sufficient interest for further functional analysis. We therefore extracted the intersection between top hubs and/bottlenecks and motif 4 members (Additional file 7) within HUBBA by calculating the shortest path between top hubs/bottlenecks and members of motif 4 which also contained Egr1 as well (Additional file 7, see sheet “intersection”). Notably, among Egr1’s neighbors, expression of Prlh and Chgb were also up-regulated but not that of CREB (Additional file 7, ‘intersection’ Excel sheet).
Gene ontology analysis of 1× network motifs and the Egr1/Stat3 neighborhood showed that it is enriched for cell cycle, cell proliferation and cytokine-cytokine receptor (Sub-network motifs 4, 6, 7, 8) and Jak/Stat pathway genes (Additional file 6) whereas 6× network motifs are enriched for genes belonging to the insulin signaling, neuroactive signaling and the purine metabolism pathways. Prlh11 is a network neighbor of Egr1 (Figure 3, right panel). Strikingly, the MAPK signaling pathway is enriched for genes in all motifs except 5 and 7 (1× network) an observation that is consistent with previous results. This finding is consistent with a role of the ERK-STAT3-Egr1 pathway in neurite outgrowth. Specifically, MEK-ERK1/2 is required for FGF1-induced neurite outgrowth, pSTAT3 (S727) and Egr1 expression in PC12 cells .
In order to confirm our findings, we tested independent samples of 1× and 6× Egr1 network neighbors for expression changes with qRT-PCR or immunoblots (Tables 4 and 5). Several of these changes were further verified in independent IMO stress samples at the mRNA level by qRT-PCR or at the protein level by western blot analysis or immunocytochemistry (Table 4). These findings confirm for the first time that Stat 3 expression (Figure 4, panel A) is significantly up-regulated with 1× IMO. Expression of Prlh1 (and of Chgb) is up-regulated in both 1× and in 6× (Figure 4, panel B). Strikingly, we observed for the first time down regulation of gene PIAS3 (Inhibitor of Stat 3, Table 4, bottom) implicating Stat3 signaling in acute IMO stress activation. PIAS3 is a small E3-type small ubiquitin-like modifier (SUMO) ligase that plays a critical role in regulating the Stat3 signaling pathway by inhibiting Stat3-DNA binding ,. The Egr1 neighbors Prlh11 (and Chgb) are significantly represented in changes in the 6× IMO samples, suggesting that they might have important functions in long-term stress signaling (Figure 4, panel B).
We have sought to identify novel members in pathways, particularly transcription factor pathways, and gene products that might be regulated in IMO stress. We approached this by performing microarray analyses of samples derived from the adrenal medullae of rats subjected to acute (1×) or repeated (6×) IMO stress in order to identify genes/gene products that are co-expressed or co-regulated with the transcription factor Egr1, which is critical for stress responses in the adrenal medulla of rats. Following identification of gene modules that are co-expressed with Egr1 with KS ranking and GSEA analysis, we mined all physiologically relevant interactions and interrelationships with text-mining, and we established that Stat3 and Prlh are significantly up-regulated in 1× and in 6× IMO stress responses suggesting for the first time that they are likely participants. Also, we confirmed earlier observations that Chgb is activated in 6× IMO stress . Stat3 is co-expressed and a network neighbor of Egr1 in the re-constructed 1× network, implicating Stat3 involvement in IMO stress. Notably, the 1× network is organized around Egr1 and Stat3 since they are top hubs (highest number of interactions/interrelationships, Figure 2, lower panel). This contention is supported by the concomitant down-regulation of the Stat3 inhibitor PIAS3 which inhibits the DNA activity of Stat3 therefore shutting down part of Stat3-mediated genomic expression changes. The observed enrichment for MAPK and cytokine genes in 1× IMO stress (see Figure 2) is consistent with the involvement of cytokines in regulating the stress response. Stat3 belongs to the STAT (signal transducers and activators of transcription) family of transcription factors that feed into the Jak/Stat signaling cascade. Several cytokine receptors regulate this cascade and this is consistent with the observed GO enrichment for cytokine receptors. Stats are phosphorylated by Jaks and activated, allowing them to translocate to the nucleus. The Jak/Stat pathway is regulated by phosphorylation/dephosphorylation by kinase/phosphatase enzymes, by Stat gene activation antagonists such as SOCS (suppressors of cytokine signaling) and by PIAS (Protein Inhibitors of Activated Stats) ,.
Prlh1 is a peptide widely distributed in the CNS and involved in mediating stress responses and activating the HPA axis . Prlh1 was found to be co-expressed with TH and PNMT in Epi synthesizing cells of the adrenal medulla , however its function in the adrenal is unclear.
In repeated stress (6× IMO), the Egr1-centered network is organized differently with Stat3 being absent and with Prlh and Chgb being close network neighbors of Egr1 (Figure 3, right panel). Moreover, the enriched motifs are strikingly different compared to the 1× Egr1 neighborhoods. The predominant motifs include members of the insulin signaling pathway and members of pathways for purine metabolism, neuroactive signaling and MAPK pathway. The latter observation is in agreement with the central role of the MAPK pathway in both acute and repeated stress responses. Chgb is up-regulated confirming previous experiments (see ref ). Induction of chromogranin B (Chgb) gene expression selectively with 6× IMO stress but not with 1×, is especially intriguing since Chgb functions in the biogenesis of secretory granules and in the sorting of proteins to the regulated secretory pathway . Chgb deficient mice display reduced levels of catecholamines released per quanta . The induction of Chgb with repeated IMO suggests formation of additional neurosecretory vesicles or production of larger quantal release which may help provide additional neurosecretory strength to adapt to further demands of chronically repeated stress . The induction of Chgb with repeated IMO suggests that formation of additional neurosecretory vesicles or production of larger quantal release may help provide additional neurosecretory strength to adapt to further demands of chronically repeated stress.
In addition to the transcription factor Egr1 which is critical for IMO stress induction in the rat adrenal medulla, we here provide evidence for the first time that the gene encoding the transcription factor Stat3 and the gene encoding the peptide Prlh1 are activated in acute (1×) and in repeated stress (6×) respectively. The data suggest that the transcription factor Egr1 has different roles in acute and in repeated stress, indicated by different networks and network neighbors and furthermore that Stat3 signaling in 1× and Prlh (and Chgb) activation cascades in 6× might also be important in transitioning from acute to repeated stress responses.
The stress procedures, isolation of RNA and Afflymetrix analysis were as previously described . Briefly, male, murine pathogen-free, Sprague−Dawley rats (280−2320 g), obtained from Taconic Farms (Germantown, NY, USA), were maintained under controlled conditions of a 12 h light−dark cycle (lights on from 6×am to 6×pm) at 23α2°C with food and water ad libitum. All animal experiments were performed in accordance with the NIH Guide for the Care and Use of Laboratory Animals and approved by the Institutional Animal Care and Use Committee.
Immobilization stress (IMO) was performed as previously described ,,. For acute stress rats were subjected to immobilization stress for 2 hrs once (1×IMO). For repeated stress, the animals were immobilized for 2 hrs daily for 6 consecutive days (6×IMO). Following the last IMO, rats were euthanized by decapitation the adrenal medulla dissected. Control groups were not exposed to stress (absolute controls). All animal manipulations were performed between 8 AM and 1×PM to control for circadian variations.
For the microarray, RNA was isolated from two separate immobilization experiments. To minimize sample variability caused by individual differences among animals, each sample was pooled from left adrenal medulla from 4 individual rats. There were 3 pooled samples for each group. RNA was extracted using Absolutely RNA Miniprep Kit (Stratagene, La Jolla, CA). The integrity of the RNA was assessed by the A260/A280 ratio which was close to 2.0 and by electrophoresis (Agilent Bioanalyzer 2100).
Gene expression analyses were performed by the NIH Neuroscience Microarray Consortium at UCLA Medical Center (Los Angeles, CA). Total RNA (≥4μ g) from each group was converted to cDNA by using superscript reverse transcriptase and the T7-Oligo (dT) promoter primer kit (Affymetrix, Inc). Following RNase H-mediated second-strand cDNA synthesis, the double-stranded cDNA were purified and served as a template in the subsequent in vitro transcription reaction (Affymetrix, Inc.). The in vitro transcription reaction was carried by T7 RNA polymerase and a biotinylated nucleotide analog/ribonucleotide mix (Affymetrix, Inc.). The biotinylated cRNA targets were purified and fragmented. Each cRNA was hybridized to an individual Affymetrix GeneChip Rat Array Expression 230 2.0 (RAE 230 2.0 array) which was subsequently processed for washing and staining with the antibody stain solution with streptavidin phycoerythrin and the arrays were scanned on the GeneChip Scanner 3000. The raw pixel data have been deposited to the Gene Expression Omnibus (GEO) database Series # GSE8184.
Computational and bioinformatic methods
Top genes in gene sets that might be co-regulated with Egr1 were identified with Kolmogorov-Smirnov (KS) statistics and gene set enrichment analysis (GSEA) , (For full description, see Additional file 2) on the genome-wide gene expression microarray data sets obtained from samples of the adrenal medulla of rats subjected to 1× or 6× IMO stress, each data set containing three samples. The strategy employed extraction of genes that are co-expressed with Egr1 in single (1×) or repeated (6×) exposure to stress (Figure 1).
Animal experiments with n = 4 per group were from at least two separate experiments. Data were analyzed with ANOVA followed by Bonferroni post-hoc analysis using GraphPad Prism 4 software (GraphPad Software, Inc., La Jolla, CA). A value of p = ≤ 0.05 was considered significant.
Analysis of variance
Epidermal growth factor
Epidermal growth factor receptor
Early growth factor 1
Egr1 positive gene index phenotype
Fibroblast growth factor 2
Gene annotation tool to help explain relationships
Gene set enrichment analysis
- HTFN :
Human transcription factor network
Hub objects analyzer
Single immobilization stress/acute IMO stress
Immobilization stress repeated daily for 6 consecutive days/prolonged IMO stress
Matrix Gia Protein
Suppressors of cytokine signaling
Prolactin releasing peptide
Protein inhibitor of activated Stat3
Real time polymerase Chain reaction
Signal transducer and activator of transcription 3
Kvetnansky R, Sabban EL, Palkovits M: Catecholaminergic systems in stress: structural and molecular genetic approaches. Physiol Rev. 2009, 89: 535-606. 10.1152/physrev.00042.2006.
McMahon A, Kvetnanský R, Fukuhara K, Weise VK, Kopin IJ, Sabban EL: Regulation of tyrosine hydroxylase and dopamine beta-hydroxylase mRNA levels in rat adrenals by a single and repeated immobilization stress. J Neurochem. 1992, 58 (6): 2124-2130. 10.1111/j.1471-4159.1992.tb10954.x.
Viskupic E, Kvetnansky R, Sabban EL, Fukuhara K, Weise VK, Kopin IJ, Schwartz JP: Increase in rat adrenal phenylethanolamine N-methyltransferase mRNA level caused by immobilization stress depends on intact pituitary-adrenocortical axis. J Neurochem. 1994, 63 (3): 808-814. 10.1046/j.1471-4159.1994.63030808.x.
Wong DL, Her S, Tai TC, Bell RA, Rusnák M, Farkas R, Kvetňanský R, Shih JC: Stress-Induced Expression of Phenylethanolamine N-Methytransferase: Normal and Knock out Animals. Stress: Neural, Endocrine and Molecular Studies. Edited by: McCarty R, Aguilera G, Sabban EL, Kvetnansky R. 2002, Taylor and Francis, London, 129-135.
Sabban EL, Kvetnansky R: Stress-triggered activation of gene expression in catecholaminergic systems: dynamics of transcriptional events. Trends Neurosci. 2001, 24 (2): 91-98. 10.1016/S0166-2236(00)01687-8.
Liu X, Serova L, Kvetnanský R, Sabban EL: Identifying the stress transcriptome in the adrenal medulla following acute and repeated immobilization. Ann N Y Acad Sci. 2008, 1148: 1-28. 10.1196/annals.1410.082.
O’Donovan KJ, Tourtellotte WG, Millbrandt J, Baraban JM: The EGR family of transcription-regulatory factors: progress at the interface of molecular and systems neuroscience. Trends Neurosci. 1999, 22 (4): 167-173. 10.1016/S0166-2236(98)01343-5.
Gashler A, Sukhatme VP: Early growth response protein 1 (Egr-1): prototype of a zinc-finger family of transcription factors. Prog Nucleic Acid Res Mol Biol. 1995, 50: 191-224. 10.1016/S0079-6603(08)60815-6.
Christy B, Nathans D: DNA binding site of the growth factor-inducible protein Zif268. Proc Natl Acad Sci U S A. 1989, 86 (22): 8737-8741. 10.1073/pnas.86.22.8737.
Thiel G, Cibelli G: Regulation of life and death by the zinc finger transcription factor Egr-1. J Cell Physiol. 2002, 193 (3): 287-292. 10.1002/jcp.10178.
Silverman ES, Collins T: Pathways of Egr-1-mediated gene transcription in vascular biology. Am J Pathol. 1999, 154 (3): 665-670. 10.1016/S0002-9440(10)65312-6.
Knapska E, Kaczmarek L: A gene for neuronal plasticity in the mammalian brain: Zif268/Egr-1/NGFI-A/Krox-24/TIS8/ZENK?. Prog Neurobiol. 2004, 74 (4): 183-211. 10.1016/j.pneurobio.2004.05.007.
Bozon B, Davis S, Laroche S: Regulated transcription of the immediate-early gene Zif268: mechanisms and gene dosage-dependent function in synaptic plasticity and memory formation. Hippocampus. 2002, 12 (5): 570-577. 10.1002/hipo.10100.
Tarcic G, Avraham R, Pines G, Amit I, Shay T, Lu Y, Zwang Y, Katz M, Ben-Chetrit N, Jacob-Hirsch J, Virgilio L, Rechavi G, Mavrothalassitis G, Mills GB, Domany E, Yarden Y: EGR1 and the ERK-ERF axis drive mammary cell migration in response to EGF. FASEB J. 2012, 26 (4): 1582-1592. 10.1096/fj.11-194654.
Zwang Y, Oren M, Yarden Y: Consistency test of the cell cycle: roles for p53 and EGR1. Cancer Res. 2012, 72 (5): 1051-1054. 10.1158/0008-5472.CAN-11-3382.
Baron V, Adamson ED, Calogero A, Ragona G, Mercola D: The transcription factor Egr1 is a direct regulator of multiple tumor suppressors including TGFbeta1, PTEN, p53, and fibronectin. Cancer Gene Ther. 2006, 13 (2): 115-124. 10.1038/sj.cgt.7700896.
Cubero FJ, Nieto N: Arachidonic acid stimulates TNFalpha production in Kupffer cells via a reactive oxygen species-pERK1/2-Egr1-dependent mechanism. Am J Physiol Gastrointest Liver Physiol. 2012, 303 (2): G228-G239. 10.1152/ajpgi.00465.2011.
Waki N, Yamane M, Yamamoto S, Okazaki M, Sugimoto S, Matsukawa A, Oto T, Miyoshi S: Egr1: a novel target for ameliorating acute allograft rejection in an experimental lung transplant model. Eur J Cardiothorac Surg. 2012, 41 (3): 669-675. 10.1093/ejcts/ezr030.
Yu X, Shen N, Zhang ML, Pan FY, Wang C, Jia WP, Liu C, Gao Q, Gao X, Xue B, Li CJ: Egr-1 decreases adipocyte insulin sensitivity by tilting PI3K/Akt and MAPK signal balance in mice. EMBO J. 2011, 30 (18): 3754-3765. 10.1038/emboj.2011.277.
Papanikolaou NA, Sabban EL: Ability of Egr1 to activate tyrosine hydroxylase transcription in PC12 cells. Cross-talk with AP-1 factors. J Biol Chem. 2000, 275 (35): 26683-26689.
Nakashima A, Ota A, Sabban EL: Interactions between Egr1 and AP1 factors in regulation of tyrosine hydroxylase transcription. Brain Res Mol Brain Res. 2003, 112 (1×2): 61-69. 10.1016/S0169-328X(03)00047-0.
Ebert SN, Balt SL, Hunter JP, Gashler A, Sukhatme V, Wong DL: Egr-1 activation of rat adrenal phenylethanolamine N-methyltransferase gene. J Biol Chem. 1994, 269 (33): 20885-20898.
Morita K, Ebert SN, Wong DL: Role of transcription factor Egr-1 in phorbol ester-induced phenylethanolamine N-methyltransferase gene expression. J Biol Chem. 1995, 270 (19): 11161-11167. 10.1074/jbc.270.19.11161.
Tai TC, Morita K, Wong DL: Role of Egr-1 in cAMP-dependent protein kinase regulation of the phenylethanolamine N-methyltransferase gene. J Neurochem. 2001, 76 (6): 1851-1859. 10.1046/j.1471-4159.2001.00189.x.
Papanikolaou NA, Sabban EL: Sp1/Egr1 motif: a new candidate in the regulation of rat tyrosine hydroxylase gene transcription by immobilization stress. J Neurochem. 1999, 73 (1): 433-436. 10.1046/j.1471-4159.1999.0730433.x.
Liu X, Kvetnansky R, Serova L, Sollas A, Sabban EL: Increased susceptibility to transcriptional changes with novel stressor in adrenal medulla of rats exposed to prolonged cold stress. Brain Res Mol Brain Res. 2005, 141 (1): 19-29. 10.1016/j.molbrainres.2005.07.019.
Lin WF, Chen CJ, Chang YJ, Chen SL, Chiu IM, Chen L: SH2B1beta enhances fibroblast growth factor 1 (FGF1)-induced neurite outgrowth through MEK-ERK1/2-STAT3-Egr1 pathway. Cell Signal. 2009, 21 (7): 1060-1072. 10.1016/j.cellsig.2009.02.009.
Chung CD, Liao J, Liu B, Rao X, Jay P, Berta P, Shuai K: Specific inhibition of Stat3 signal transduction by PIAS3. Science. 1997, 278 (5344): 1803-1805. 10.1126/science.278.5344.1803.
Husby J, Todd AK, Haider SM, Zinzalla G, Thurston DE, Neidle S: Molecular dynamics studies of the STAT3 homodimer: DNA complex: relationships between STAT3 mutations and protein-DNA recognition. J Chem Inf Model. 2012, 52 (5): 1179-1192. 10.1021/ci200625q.
Sabban EL, Nankova BB, Serova LI, Kvetnansky R, Liu X: Molecular regulation of gene expression of catecholamine biosynthetic enzymes by stress: sympathetic ganglia versus adrenal medulla. Ann N Y Acad Sci. 2004, 1018: 370-377. 10.1196/annals.1296.046.
Nankova B, Devlin D, Kvetnanský R, Kopin IJ, Sabban EL: Repeated immobilization stress increases the binding of c-Fos-like proteins to a rat dopamine beta-hydroxylase promoter enhancer sequence. J Neurochem. 1993, 61 (2): 776-779. 10.1111/j.1471-4159.1993.tb02188.x.
Nankova BB, Rivkin M, Kelz M, Nestler EJ, Sabban EL: Fos-related antigen 2: potential mediator of the transcriptional activation in rat adrenal medulla evoked by repeated immobilization stress. J Neurosci. 2000, 20 (15): 5647-5653.
Sabban EL, Liu X, Serova L, Gueorguiev V, Kvetnansky R: Stress triggered changes in gene expression in adrenal medulla: transcriptional responses to acute and chronic stress. Cell Mol Neurobiol. 2006, 26 (4-6): 845-854.
Kvetnansky R, Kubovcakova L, Tillinger A, Micutkova L, Krizanova O, Sabban EL: Gene expression of phenylethanolamine N-methyltransferase in corticotropin-releasing hormone knockout mice during stress exposure. Cell Mol Neurobiol. 2006, 26 (4-6): 733-752. 10.1007/s10571-006-9063-7.
Sabban EL, Tillinger A, Nostramo R, Serova L: Stress triggered changes in expression of genes for neurosecretory granules in adrenal medulla. Cell Mol Neurobiol. 2012, 32 (5): 795-800. 10.1007/s10571-011-9785-z.
Ye L, Wang X, Metzger DS, Riedel E, Montaner LJ, Ho W: Upregulation of SOCS-3 and PIAS-3 impairs IL-12-mediated interferon-gamma response in CD56 T cells in HCV-infected heroin users. PLoS ONE. 2010, 5 (3): e9602-10.1371/journal.pone.0009602.
Greenhalgh CJ, Hilton DJ: Negative regulation of cytokine signaling. J Leukoc Biol. 2001, 70 (3): 348-356.
Sun B, Fujiwara K, Adachi S, Inoue K: Physiological roles of prolactin-releasing peptide. Regul Pept. 2005, 126 (1×2): 27-33. 10.1016/j.regpep.2004.08.008.
Fujiwara K, Matsumoto H, Yada T, Inoue K: Identification of the prolactin-releasing peptide-producing cell in the rat adrenal gland. Regul Pept. 2005, 126 (1×2): 97-102. 10.1016/j.regpep.2004.08.039.
Natori S, Huttner WB: Chromogranin B (secretogranin I) promotes sorting to the regulated secretory pathway of processing intermediates derived from a peptide hormone precursor. Proc Natl Acad Sci U S A. 1996, 93 (9): 4431-4436. 10.1073/pnas.93.9.4431.
Borges R, Díaz-Vera J, Domínguez N, Arnau MR, Machado JD: Chromogranins as regulators of exocytosis. J Neurochem. 2010, 114 (2): 335-343. 10.1111/j.1471-4159.2010.06786.x.
Kvetnansky R, Mikulaj L: Adrenal and urinary catecholamines in rats during adaptation to repeated immobilization stress. Endocrinology. 1970, 87 (4): 738-743. 10.1210/endo-87-4-738.
Mootha VK, Lindgren CM, Eriksson KF, Subramanian A, Sihag S, Lehar J, Puigserver P, Carlsson E, Ridderstråle M, Laurila E, Houstis N, Daly MJ, Patterson N, Mesirov JP, Golub TR, Tamayo P, Spiegelman B, Lander ES, Hirschhorn JN, Altshuler D, Groop LC: PGC-1alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat Genet. 2003, 34 (3): 267-273. 10.1038/ng1180.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, Mesirov JP: Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005, 102 (43): 15545-15550. 10.1073/pnas.0506580102.
We are grateful to Dr. Lidia Serova, New York Medical College for her assistance with the stress experiments and very useful suggestions. The microarray study was supported by NIH grant NS28869 (to ES).
The authors declare that they have no competing interests.
NAP conceived the idea, executed all computational and bioinformatic methods (KS statistics, module analysis with GSEA, GO analysis and network re-construction) and oversaw all aspects of the work and the manuscript. AT and XL performed all experiments except for microarray experiments which were performed at the NIH Neuroscience Microarray Consortium, UCLA Medical Center, CA. AGP has assisted with interpretation of the results and with the manuscript and ELS has managed all aspects of the work. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Figure S1.: Schematic representation of the KS ranking and GSEA procedures. The normalized enrichment score (ES) can be positive (for gene sets that are enriched and therefore correlate with the expression profile) or negative for anti-correlating profiles. (JPEG 60 KB)
Additional file 2:Computational and bioinformatic methods, specifically for Kolomogorov-Smirniv statistics and for gene set expression analysis.(DOCX 37 KB)
Additional file 3:GSEA 1× IMO data. This file contains the GSEA parameters used in analyzing 1× microarray samples, the ranked gene list and their KS scores, the top Egr1 gene index and categorical gene lists and their corresponding top gene sets with their ES scores. (XLS 3 MB)
Additional file 4:GSEA 6× IMO data. Same as for Additional file 3, except the data are for 6× samples. (XLS 3 MB)
Additional file 5:Gene Interaction-Interrelationship Modules. This file contains the two gene module lists generated with GSEA for 1× and 6× samples respectively, and the list of Egr1 co-expressed genes extracted from the CoxpressDb. (XLS 50 KB)
Additional file 6:1× IMO Network Motif GO Categories and Network Data. Additional file 6 contains the top motifs for 1× IMO samples, the Gene Ontology enrichment data for the motifs, motif figures, top hubs and bottlenecks. (XLS 602 KB)
About this article
Cite this article
Papanikolaou, N.A., Tillinger, A., Liu, X. et al. A systems approach identifies co-signaling molecules of early growth response 1 transcription factor in immobilization stress. BMC Syst Biol 8, 100 (2014). https://doi.org/10.1186/s12918-014-0100-8
- Adrenal medulla