- Research article
- Open Access
The duration of gastrin treatment affects global gene expression and molecular responses involved in ER stress and anti-apoptosis
BMC Genomics volume 14, Article number: 429 (2013)
How cells decipher the duration of an external signal into different transcriptional outcomes is poorly understood. The hormone gastrin can promote a variety of cellular responses including proliferation, differentiation, migration and anti-apoptosis. While gastrin in normal concentrations has important physiological functions in the gastrointestine, prolonged high levels of gastrin (hypergastrinemia) is related to pathophysiological processes.
We have used genome-wide microarray time series analysis and molecular studies to identify genes that are affected by the duration of gastrin treatment in adenocarcinoma cells. Among 403 genes differentially regulated in transiently (gastrin removed after 1 h) versus sustained (gastrin present for 14 h) treated cells, 259 genes upregulated by sustained gastrin treatment compared to untreated controls were expressed at lower levels in the transient mode. The difference was subtle for early genes like Junb and c-Fos, but substantial for delayed and late genes. Inhibition of protein synthesis by cycloheximide was used to distinguish between primary and secondary gastrin regulated genes. The majority of gastrin upregulated genes lower expressed in transiently treated cells were primary genes induced independently of de novo protein synthesis. This indicates that the duration effect of gastrin treatment is mainly mediated via post-translational signalling events, while a smaller fraction of the differentially expressed genes are regulated downstream of primary transcriptional events. Indeed, sustained gastrin treatment specifically induced prolonged ERK1/2 activation and elevated levels of the AP-1 subunit protein JUNB. Enrichment analyses of the differentially expressed genes suggested that endoplasmic reticulum (ER) stress and survival is affected by the duration of gastrin treatment. Sustained treatment exerted an anti-apoptotic effect on serum starvation-induced apoptosis via a PKC-dependent mechanism. In accordance with this, only sustained treatment induced anti-apoptotic genes like Clu, Selm and Mcl1, while the pro-apoptotic gene Casp2 was more highly expressed in transiently treated cells. Knockdown studies showed that JUNB is involved in sustained gastrin induced expression of the UPR/ER stress related genes Atf4, Herpud1 and Chac1.
The duration of gastrin treatment affects both intracellular signalling mechanisms and gene expression, and ERK1/2 and AP-1 seem to play a role in converting different durations of gastrin treatment into distinct cellular responses.
Although the cellular response to one particular ligand is usually specific to a given cell type , some studies have shown that the duration [2–5] or concentration [4, 6] of the external stimuli may also influence the transcriptional program and cell fate decision. Gene expression profiling in a pancreatic beta cell model treated with glucose and cAMP has shown that the majority of genes regulated in sustained treated cells were not regulated by transient (1 h) treatment, indicating that beta cells can produce drastically different transcriptional outputs in response to different durations of metabolic stimuli . Glauser et al.  later showed that sustained elevated glucose promotes long-term phosphorylation of extracellular signal-regulated kinases (ERK1/2) followed by sustained changes in the activator protein 1 (AP-1 ) subunits composition as well as in AP-1 controlled gene expression. Hence, the authors suggested that in transiently (1 h) treated cells the ERK1/2 activation is too short to stabilize downstream gene expression. In human keratinocytes (HaCaT), sustained transforming growth factor beta (TGFβ) treatment induces a more persistent phosphorylation of the primary intracellular mediator SMAD2 than short pulses, and this may be critical for cell fate determination like cell growth arrest in this cell system . Others have observed that a group of tumour necrosis factor (TNF) upregulated genes which remain at high levels in the sustained mode for 10 h, quickly returned to baseline if TNF was removed after 6 h . However, in which way the duration of external signals affects molecular and biological responses has only been explored in a few cell systems and the mechanisms are still not well characterized.
The peptide hormone gastrin is the central regulator of gastric acid secretion and plays a prominent role in regulation of growth and differentiation of gastric and colonic mucosa [7, 8]. Gastrin signals via the gastrin/cholecystokinin-2 (CCK2) receptor , and promotes a variety of cell or tissue specific outcomes including proliferation, survival, anti-apoptosis, differentiation and migration [7, 10]. In normal physiological conditions, gastrin levels are transiently upregulated in response to a meal. Prolonged elevated blood levels of gastrin (hypergastrinemia) can occur as a consequence of e.g. atrophic gastritis or pharmacologic inhibition of gastric acid secretion, which interrupts negative feedback mechanisms on gastrin producing G-cells residing in the gastric mucosa . Recently it was shown that the gastrin promotor can be activated by disease associated Helicobacter pylori strains via the EGFR/Raf/MEK/ERK cascade . Both hypergastrinemia and mutational activation of the CCK2 receptor have been linked to development of neuroendocrine gastrointestinal tumours (carcinoids) ; and gastrin and CCK2R are reported to be upregulated in human pancreatic adenocarcinoma [12, 13].
Since transiently increased gastrin levels have important physiological functions in the gastrointestine, while sustained high gastrin levels (hypergastrinemia) are related to pathophysiological processes [7, 8, 10, 14], it is of interest to examine how the duration of gastrin treatment affects gene expression and molecular responses. We have therefore conducted the present study to examine how adenocarcinoma cells respond to transient versus sustained gastrin signalling and to identify characteristic differences between downstream biological responses. To do this, we analysed genome-wide time series data of gastrin-regulated gene expression which compared treatment in the transient (1 h) versus sustained (14 h) mode. We used data from additional gastrin response time series experiments to identify dependence on de novo protein synthesis related to transcriptional timing (early, delayed, late) for genes that are upregulated in the sustained mode (compared to untreated controls) and lower expressed in transiently treated cells. Enrichment analysis of genes differentially regulated in transient versus sustained gastrin signalling suggested that pathways related to ER stress, cell survival and anti-apoptosis were affected by duration of gastrin signalling. Indeed, several genes known to be involved in these pathways were expressed at lower levels in cells treated in the transient mode as compared to cells treated in the sustained mode. Furthermore, the anti-apoptotic effect of gastrin on serum starvation-induced apoptosis was dependent on sustained treatment, and our results indicate that the anti-apoptotic effect of gastrin involves PKC/ERK1/2 signalling. Sustained gastrin treatment induced prolonged ERK1/2 activation and elevated levels of the AP-1 subunit protein JUNB. The important role of ERK1/2 and AP-1 in converting different durations of gastrin treatment into distinct transcriptional responses was strengthened by our findings that knock down of JUNB reduced gastrin-mediated transcriptional activation of genes related to the ER stress and survival pathways. Overall, our work may contribute to a better understanding on how a cell deciphers the durations of gastrin induced signalling into specific cellular responses.
Results and discussion
Identification of genes that are differentially regulated by transient versussustained gastrin signalling
To investigate the mRNA transcriptome in response to varying duration of gastrin treatment, we performed genome-wide microarray time series experiments in the adenocarcinoma cell line AR42J. We treated the cells in a sustained mode (14 h of continuous presence of gastrin) and in a transient mode (gastrin was removed after 1 h of treatment), as illustrated in Figure 1A. Statistical analysis of the time series gene expression responses identified 259 genes with lower and 144 genes with higher expression in cells treated with gastrin in transient versus sustained mode (Additional file 1: Table S1). The 403 differentially expressed genes were visualized in a cluster analysis heat map (Figure 1B). The upper part of the heat map shows that the 144 genes expressed at higher levels in the transient mode are downregulated by sustained gastrin treatment compared to untreated controls, while the 259 genes with lower expression in the transient mode (lower part of heat map) are upregulated by sustained gastrin treatment. The heat map further illustrates that differences in gene expression between cells treated in the transient and sustained mode were subtle at 1.5 h and 2 h, while substantial differences were apparent at time points after 2 h. The most striking overall trend for these 403 genes is that their transient mode expression levels after 4 h are similar to their levels in untreated cells, while they are clearly elevated or reduced in the sustained mode. These observations comply with our results from an earlier microarray 24 h time series study which showed that the regulation of several hundred genes was affected by gastrin signalling duration and that globally, gene expression levels returned earlier to baseline in cells subjected to transient (2 h) treatment compared to sustained treatment (Additional file 2: Accession number: E-MTAB-123). Thus, our results indicate that sustained gastrin signalling induces a transcriptional programme that differs significantly from the transcriptional response to transient gastrin signalling.
To search for molecular responses affected by the duration of gastrin treatment, the 403 genes exhibiting differential regulation in transiently versus sustained treated cells (Additional file 1) were subjected to pathway enrichment analysis using the MetaCore tool of GeneGo package . We found that high scoring canonical pathway and networks were related to apoptosis-survival-development as well as to stress responses with a focus on the unfolded protein response (UPR), endoplasmic reticulum (ER)-stress, DNA damage and cell cycle (Table 1). Gene Ontology (GO) biological process enrichment analysis of the 259 genes with lower expression in transiently versus sustained treated cells also identified cellular responses to unfolded proteins as significantly enriched (Additional file 3: Table S2). Enriched GO processes associated with the 144 genes with higher expression in the transient mode were related to cell cycle phases and mitosis (Additional file 3: Table S3). Thus, it seems that sustained gastrin treatment is required for activation of genes involved in ER stress/survival (Additional file 3: Table S2), while many genes that are less downregulated by transient treatment encode proteins that promote cell proliferation (Additional file 3: Table S3).
ER stress/UPR is increasingly recognized to play a role in tumourigenesis as well as in cell homeostasis [16–18]. Thus, our novel findings that gastrin i) regulate ER stress/UPR genes and ii) does so in a signal duration-specific manner are highly interesting. Depending upon the duration and the degree of ER stress, the UPR can provide either survival signals by activating adaptive and anti-apoptotic pathways, or death signals by inducing cell death programs . As illustrated by the selected genes shown in Figure 2, sustained gastrin treatment was required for increased expression of the main regulator of the UPR, Hspa5 (BiP/Grp78), the transcription factor Atf4 and the downstream mediators Herpud1, Ddit3 (Chop) and Chac1[18–21]. Since HERPUD1 is reported to mediate cell protection via the stabilization of ER Ca2+ homeostasis , its activation indicates a possible mechanism for enhanced cell survival by sustained gastrin treatment. In contrast, the ATF4 target gene Ddit3 (Chop), and DDIT3 target gene Chac1 encode proteins associated with the pro-apoptotic effect of the UPR [21, 23]. Our results demonstrate that sustained gastrin treatment is required for induction of several ER stress/UPR genes and that the proteins encoded by these genes are associated with both pro-survival and pro-apoptotic effects. Thus, it cannot be directly deduced from these results in which direction the subsequent cellular response will be affected.
The anti-apoptotic effect of gastrin is well documented [24–27], and it is therefore of interest to take a closer look at the relation to apoptosis among the genes differentially regulated in the transient versus sustained mode. The pro-apoptotic Casp2 exhibits higher expression levels in the transient mode than in the sustained mode, while the anti-apoptotic genes Mcl1, Itpr1, Selm and Clu were more strongly increased by sustained gastrin treatment (Figure 3). We have previously shown that CLU is required for the anti-apoptotic effect of gastrin in our AR42J adenocarcinoma cell line model . As a versatile stress-induced chaperone, CLU is suggested to be involved in protein homeostasis via unfolded protein and ER stress responses , and was recently shown to function cooperatively with GRP78 (Bip) to mediate anti-apoptotic effect in the mitochondria pathway . MCL1, a member of the Bcl-2 protein family which inhibits release of cytochrome c from mitochondria , has previously been shown to be involved in the anti-apoptotic effect of gastrin in human adenocarcinoma cells (AGS-GR) stably transfected with the gastrin/CCK-2 receptor . The ITPR1 receptor mediates Ca2+ release from ER and may be involved in mediating anti-apoptosis by interacting with anti-apoptotic members of the Bcl-2 family, including MCL1 , while SELM is shown to decrease Ca2+ release from ER in response to oxidative stress and reduce apoptotic cell death . Taken together, our gene expression microarray analysis indicates that sustained but not transient gastrin treatment affects the balance between apoptosis and survival via several cross linked mechanisms including communication between ER and mitochondria which is critical for cellular decision making .
Sustained gastrin treatment inhibits serum starvation-induced apoptosis via a PKC-dependent mechanism
To examine the effect of gastrin signalling duration on cell survival, we induced apoptosis in AR42J cells by serum starvation for 72 h and measured caspase activity in cells subjected to sustained or transient gastrin treatment. We found that sustained gastrin treatment significantly reduced effector caspase activity (Figure 4A), which is in accordance with our previous results using both TUNEL and caspase assay for detection of apoptotic cells . This anti-apoptotic effect was not detected in cells treated with gastrin in a transient mode (Figure 4B). Thus, manifestation of the gastrin-induced anti-apoptotic effect which is well documented in literature [24–27], requires sustained gastrin treatment and is not induced when cells are exposed to gastrin for only 1 hour.
Since both the PI3K/AKT (PKB) pathway and ERK1/2 activation have been shown to be important in mediating the biological effects of gastrin, including anti-apoptosis [25, 26], we next set out to determine the involvement of these two signalling pathways in our anti-apoptosis model system. Involvement of ERK1/2-signalling was assayed by the use of a PKC-inhibitor, since ERK1/2 activation has been shown to be mediated by PKC-Src/Ras/Raf/MEK/ERK1/2 cascade and/or by direct activation of Raf by PKC in the gastrin response . The anti-apoptotic effect of sustained gastrin treatment was not influenced by the PI3K inhibitor (Figure 4C). In the presence of the PKC inhibitor gastrin still resulted in a small decrease in caspase 3/7 activity. However, this decrease was smaller than in the absence of the PKC-inhibitor, and the anti-apoptotic effect was no longer statistically significant (Figure 4D). We therefore conclude that PKC and its downstream effector ERK1/2 may play a decisive role in the sustained gastrin-induced anti-apoptotic effect. This observation complies with results reported by others showing that a specific inhibitor of the ERK1/2 activator MEK1 blocked gastrin-induced anti-apoptosis in AR42J cells .
Genes that require sustained gastrin treatment for late upregulation include both primary and secondary responders
Cellular decision making in response to growth factors like gastrin is driven by transcriptional cascades involving primary response genes dominated by regulators of transcription and signal transduction and secondary response genes which depend on de novo protein synthesis and which include a high fraction of biological response effector genes [34, 35]. To investigate primary and secondary transcriptional response mechanisms underlying the distinct biological responses to transient versus sustained gastrin signalling, we focused on genes that were markedly upregulated by gastrin and expressed at lower levels in the transient mode. This 181 gene subset can be grouped according to their time profiles as follows: 15 early genes peaking before 2 h (Figure 5 a-b); 50 delayed genes with peak expression at 2–4 h (Figure 5, c-d); and 116 late genes upregulated 4–14 h (Figure 5, e-f)(see Additional file 4: Table S4 for details). Genes dependent on de novo protein synthesis were identified by investigating the effect of the protein synthesis inhibitor CHX on gastrin-induced gene expression responses. This was done using data from a genome wide time series experiment where the cells were treated by gastrin in the absence and presence of CHX. To minimize misinterpretation due to the confounding effect of CHX [34, 36, 37], the temporal profiles (7 time points; 0–10 h) of treated and control cells were manually evaluated: For those genes where the gastrin response in the presence of CHX was higher than the response with CHX alone, the gene was classified as a primary gene not dependent on de novo protein synthesis for gastrin-induced expression. Gastrin induced genes that were not upregulated in the presence of CHX were classified as secondary genes, i.e., genes whose transcriptional activation depends on de novo protein synthesis of one or several factors that must be induced by gastrin previous to these secondary genes. In cases where the effect of gastrin could not be separated from the confounding effect of CHX, the gene was classified as uncertain. Among the subset of 181 genes upregulated in sustained treated cells and lower expressed in transiently treated cells, 115 were classified as primary genes, 48 as secondary and 18 as uncertain. The fraction of primary genes was successively reduced from 100% among the early genes to 40% among the late genes, while 43 of the 48 secondary genes were found among the late gastrin induced genes (Table 2 and Additional file 4: Table S4). Temporal expression data for genes in the subset are shown in Additional file 5.
Even though we found that the fractions of primary genes are reduced throughout the time course, most (~70%) of the genes in the subset were classified as primary genes that do not depend on de novo protein synthesis of upstream regulators (Table 2 and Additional file 4: Table S4). This finding is somewhat surprising, since previous reports have indicated that signal-responsive genes are induced in waves where the group of genes that requires de novo protein synthesis for expression and that is normally expressed in later waves, is far more numerous than the group of early primary response genes . One explanation for our apparently conflicting observations implying that a large fraction of delayed and late genes are indeed primary genes, is that we recorded the genome-wide transcriptome at a high number of time points up to 14 h and in this way are able to display the gastrin-induced transcriptional cascades at high temporal resolution. This approach revealed that in addition to early genes transcriptionally activated within 2 h, also a high proportion of the delayed and late gastrin induced genes were upregulated independently of de novo protein synthesis. The observation that many delayed induced genes were primary genes are in accordance by results from a global time series gene expression analysis in human glioblastoma cells where Tullai et al.  showed that of the 133 genes induced within 4 h of human platelet-derived growth factor (PGDF) treatment, 49 were primary immediate early genes and 58 were primary delayed genes. Detailed analyses of the relation between temporal gene expression profiles beyond 4 h and dependence on de novo protein synthesis have to our knowledge not been reported. Thus, our present study contributes to shed light on these aspects of growth factor response transcriptional dynamics.
To further characterize the subset of primary and secondary gastrin induced genes that were lower expressed in the transient mode, we annotated their molecular function based on Gene Ontology (GO) and literature and grouped them into five main categories: regulation of gene expression, signal transduction, transporters, protein binding and enzymes. Regulators of gene expression are mainly found among the early and delayed genes, while genes encoding signal transduction proteins are spread throughout all temporal profile groups. Late genes are dominated by effectors involved in transport, enzyme and chaperone activity (see Table 2 and Additional file 4 for details).
The group of late genes upregulated only by sustained gastrin treatment comprises 65 primary and 43 secondary mRNA transcripts. Both primary and secondary late gene products are involved in transport and protein binding including chaperone activity. In addition, 35 late genes are classified as enzymes (Table 2 and Additional file 4: Table S4). Among these, enzymes are possible effector genes specifically associated with processes in the exocrine pancreas: The protease coding genes Ctrc, Prss1, Prss3, Sec11c and Cela3b are all highly expressed in the pancreatic acinar derived AR42J cell-line. Sustained but not transient gastrin treatment significantly enhanced the mRNA levels of these enzymes (Additional file 5). The genes were not induced in the presence of CHX and were therefore classified as secondary gastrin-induced genes. To our knowledge, gastrin has not previously been shown to activate gene expression of these proteases.
Proteases can contribute to cancer development by several mechanisms . For example, in human pancreatic cancer, PRSS3 upregulates the ligand VEGFA via the proteinase-activated receptor 1 (PAR1)-mediated ERK1/2 pathway; and blockade of PRSS3-triggered ERK signalling is shown to delay the progression of metastasis and prolong the survival of mice bearing PRSS3-positive human pancreatic tumours . Others have shown that gastrin enhances VEGFA gene expression in human colon cancer cells stably transfected with a wild-type CCK2 receptor . Moreover, somatic mutations that increase gastrin/CCK2 receptor activity increase the secretion of VEGF and promote cell migration and angiogenesis in colorectal and gastric cancer cells . Our findings that upregulation of both Prss3 and Vegfa require sustained gastrin signalling indicate that the proteins encoded by these genes may be involved in hypergastrinemia associated pathophysiological processes (see Figure 6 and discussion below).
Overall, scrutiny of the 181 selected temporal gene expression profiles revealed that the differences in transiently versus sustained treated cells were subtle for early response genes and more substantial for delayed or late gastrin induced genes; and genes that are only upregulated in the sustained mode comprise both primary and secondary genes (see Additional files 4 and 5 for extensive analysis). Exemplified time profiles for primary and secondary genes are shown in Figure 6. Primary early genes like c-Fos and Junb were upregulated in both sustained and transiently gastrin treated cells. However, their mRNA levels returned to baseline earlier in cells treated in a transient versus a sustained mode (Figure 6A). Delayed gastrin-responsive genes like primary Hdac5 and secondary Vegfa were also temporarily upregulated in the transient mode compared to unstimulated control at early time points, but only sustained gastrin treatment induced high and prolonged expression (Figure 6B). Interestingly, both late primary genes like Maged2 and late secondary genes like Prss1 and Prss3 were only detected in cells treated in a sustained mode (Figure 6C). This suggests that gastrin directs duration-dependent gene expression via two different routes: i) via signal transduction mechanisms that directly trigger transcriptional activation of primary genes throughout both early and late stages of the time course and ii) via mechanisms that depend on increased expression of transcriptional regulators responsible for upregulation of the late secondary genes.
Sustained gastrin treatment is required for extended ERK1/2 activation
Since our data indicated that a large fraction of gastrin induced genes differentially expressed in transiently versus sustained treated cells are regulated by post-translational direct signalling mechanisms, it was of interest to examine the effect of the duration of gastrin treatment on signalling pathways known to be activated by the gastrin/CCK2 receptor . We therefore analysed AKT and ERK1/2 phosphorylation in cells treated with gastrin in a transient or a sustained mode; starting at 75 minutes, i.e. 15 min after gastrin was removed in the transient protocol. Our results show that whereas ERK1/2 activity remains elevated for several hours in cells treated in a sustained mode, the level of phosphorylated ERK1/2 returned to base-line within 60 min after gastrin was removed in transiently treated cells (Figure 7A). The level of phosphorylated AKT was low in both untreated and gastrin treated cells in this time period (Figure 7B), which is in accordance with finding by others showing that gastrin induced AKT phosphorylation peaks 5–30 min after treatment in AR42J cells  (i.e., in a time interval before the 75 minutes start of our observations). Thus, while AKT phosphorylation is not prolonged by sustained gastrin treatment, ERK1/2 phosphorylation, which is known to involve upstream PKC , is strongly elevated at later time points in sustained but not transiently treated cells.
It is well known that the duration of ERK1/2 activation is involved in directing cellular outcomes. For example, nerve growth factor (NGF) and heregulin (HRG) induce prolonged phosphorylation (i.e., activation) of ERK1/2 and subsequent differentiation in PC-12 and MCF-7 cells, respectively. Epidermal growth factor (EGF), on the other hand, induces transient phosphorylation of ERK1/2 and proliferation in the same cells [6, 48–52]. Similarly, our findings that sustained gastrin treatment is required for the sustained activation of the PKC /ERK1/2-pathway as well as for PKC-dependent activation of anti-apoptosis (Figure 4), indicate that the duration of ERK1/2 activation translates duration of gastrin signalling into pro-survival cell fate in the AR42J adenocarcinoma cell model system.
The magnitude of ERK activity can be regulated by scaffold proteins . We have recently shown that the gastrin induced MAPK scaffold protein MEK partner 1 (MP1) is important for gastrin induced phosphorylation of ERK1 and ERK2 in the human gastric adenocarcinoma cell line AGS-GR. Qualitative and quantitative differences in ERK activity can also be regulated at the receptor level . For example, the EGF receptor has been reported to undergo rapid internalization and degradation upon EGF stimulation with subsequent transient ERK activation . Less is known about how prolonged gastrin stimulation affect the CCK2 (gastrin) receptor. However in our genome-wide time series experiments we find that gastrin induces gene expression of CCK2R and that the induction of Cck2r mRNA is somewhat higher in sustained versus transiently treated cells. The Cck2r was highly expressed in the presence of CHX which indicates that Cck2r is a primary gastrin-induced gene and negatively feedback-regulated at the mRNA level (data not shown). This is in accordance with findings by others showing that gastrin increases CCK2R expression at the mRNA and protein levels, both in cell cultures [54, 55] and in vivo. More investigation is needed to address the status of the gastrin receptor upon different duration of gastrin treatment. However, our finding that sustained gastrin treatment induces prolonged ERK1/2 activation suggests that the gastrin receptor desensitizes relatively slowly compared to the EGF receptor.
Molecular responses dependent on sustained gastrin treatment involve AP-1
Sustained MEK/ERK activation is known to enhance the activity of several members of the JUN- and FOS-family of proteins [3, 48] which are involved in transcription regulation through homo- or heterodimeric AP-1 complexes that can also involve proteins from the closely related CREB/ATF- and Maf-families . Thus, it was of interest to examine the potential involvement of AP-1 in gastrin-mediated transcriptional activation of the primary and secondary genes described in Figure 5 and Table 2. We found that binding sites for AP1, ATF and CREB were overrepresented (Bonferroni corrected p value≤0.05) among primary upregulated genes (Additional file 6: Table S5). None of these transcription factor binding sites were overrepresented among the secondary genes. Interestingly, binding sites for AP-1 were overrepresented also among primary genes with late expression profiles (d-f, Figure 5). These findings indicate that AP-1 transcription factors may play a central role in gastrin induced gene expression that requires sustained treatment and is mediated via post translational signalling events independent of de novo protein synthesis.
Detailed analysis of JUNB mRNA and protein showed that although the difference at mRNA level was subtle (Figure 7C), protein levels were substantially higher and more prolonged in cells treated in the sustained mode (Figure 7D). A potential mechanism underlying this observation may be stabilization of JUNB by phosphorylation due to prevention of its degradation via proteasome pathways. Such inhibition of proteasome degradation by sustained ERK1/2 activation has been reported in other model systems [3, 48]. Our results suggest that enhanced JUNB protein levels downstream of ERK1/2 signalling may play a decisive role in differential gene expression responses to sustained versus transient gastrin treatment.
To further examine the role of JUNB, we studied gastrin-induced gene expression in cells where JUNB had been knocked down by retrovirus-based shRNA (JUNB KD). The efficiency of the RNAi was demonstrated by the fact that Junb mRNA was suppressed 70% in JUNB KD compared to control cells (Firefly KD) (Figure 8). By focusing our analysis on UPR related genes differentially expressed in transiently versus sustained gastrin treated cells (Figure 2A), we found that JUNB knock down suppressed gastrin induced transcriptional activation of the transcription factor Atf4, involved in the PERK-eIF2α pathway of the UPR, and the downstream mediator genes Herpud1 and Chac1. Gastrin-induced Clu mRNA levels were similar in JUNB KD and control cells indicating that JUNB is not critically involved in the regulation of this gene (Figure 8).
These results suggest that JUNB is involved in sustained gastrin induced expression of Atf4 which may subsequently affect downstream ATF4 mediated gene expression, thus contributing to the differences seen between sustained and transient gastrin treatment. JUNB can regulate gene transcription both positively and negatively as partners in AP-1 dimers [57, 58]. The stress induced ATF4 transcription factor has numerous dimerization partners including members of AP-1 and C/EBP family of proteins . In addition to Herpud1 and Chac1, ATF4 is known to mediate upregulation of several of the sustained gastrin regulated genes identified in this work, e.g. Atf3, Mcl1 and Vegfa[59, 60]. Recently, JUNB was found to inhibit stress-induced apoptosis , and ATF3 has been identified as a downstream target of JUNB in the survival mechanism . Taken together, our results indicate that the distinct mRNA expression patterns in cells treated in transient and sustained modes are, at least in part, linked to the duration of ERK1/2 activation and expression of the AP-1 components like JUNB as well as the ATF subfamily proteins ATF4 and ATF3 for which the mRNA expression levels are prolonged only in the sustained mode.
Genome-wide time series analyses are an important approach to capture the dynamics of stimuli-induced gene expression responses over time . The major contribution of this study is that we by extensive time series gene expression analysis and molecular studies characterize the distinct effects of duration of gastrin treatment on intracellular signalling events and gene expression in adenocarcinoma cells. We further show that the differences in global gene expression were reflected in distinct cellular responses in transiently and sustained gastrin treated cells since our results point to the possibility that a sustained high level of gastrin may affect the survival/apoptosis balance. Our findings may contribute to a better understanding of how gastrin-mediated regulation of gene expression relates to pathophysiological processes in response to high and sustained levels of gastrin. Although we have limited our investigation to one cell line in the present study, our findings provide interesting new hypotheses for further studies in vivo as well as other model systems.
Cells and reagents
The AR42J cells (rat pancreatic acinar cell derived, ATCC, Rockville, MD) were cultured at 37°C in humidified 5% CO2 Dulbecco’s modified Eagle’s medium (DMEM) containing 4.5 g/l glucose supplemented with 1% (v/v) penicillin-streptomycin, 0.1 mg/ml L-glutamine, 1 μg/ml fungizone, 1mM sodium pyruvate (all Gibco, Invitrogen), and 15% fetal bovine serum (FBS, Lonza BioWhittaker, Basel, Switzerland). Gastrin-17 (G-17) was purchased from Sigma-Aldrich (St. Louis, MO). LY 294002 and GF 109203× were obtained from Calbiochem (La Jolla, CA).
Time-series experiments in AR42J cells
The rat pancreatic acinar cell derived cell line AR42J has been used as a model system to study exocrine pancreas . Furthermore, AR42J cells express gastrin receptors endogenously and can therefore be used as a model system to study gastrin responses like proliferation , differentiation  and apoptosis [24–27]. In the time-series experiment we treated cells with 10 nM gastrin in accordance with several other studies investigating gastrin-induced responses in the AR42J cell line [25, 54, 67]. The cells were grown in 6-well plates (3 ×105 cells/well) for 72 h. Then the growth medium was replenished with 2 ml serum-free DMEM, and the cells serum starved for 20–24 h before adding gastrin. Treated and untreated cells were grown in parallel and harvested (pool of 2–3 technical replicates) at several time points, as indicated in the figures. In experiments with the protein synthesis inhibitor cycloheximide (CHX), pre-treatments with CHX (10 μg/ml) were initiated 30 min before gastrin was added. In experiments with transient versus sustained gastrin treatment, the growth medium of untreated and gastrin treated cells was removed 1 h after gastrin treatment; the cells were then washed with serum-free medium before fresh serum-free medium with gastrin (sustained gastrin treated cells) or without gastrin (transiently gastrin treated or untreated cells) was added. The cells were harvested for RNA isolation or protein lysate as described below.
Total RNA from AR42J cells and frozen oxyntic samples was isolated using RNeasy Mini kit (Qiagen, Germantown, MD) according to the manufacturer’s instructions. RNA quantity and purity was assessed using a Nanodrop spectrophotometer (Nanodrop Technologies, Rockland, DE), and RNA integrity controlled by measuring RIN values using a Bioanalyser capillary gel electrophoresis assay (Agilent Technologies, Palo Alto, CA). The RNA samples were kept at −80°C until further processing.
Genome-wide gene expression analysis on Illumina ExpressionBeadChips
RNA amplifications and hybridization were performed at the NTNU Genomics Core Facilely (GCF). Briefly, RNA was amplified with Ambion's Illumina® TotalPrep RNA Amplification kit (cat no AMIL1791) using 400 ng of total RNA as input material. The in vitro transcription (IVT) amplification that incorporated biotin-labeled nucleotides was performed overnight (14 hours) at 37°C. After the amplifications the cRNA concentrations where checked with NanoDrop ND-1000 and cRNA quality was controlled by BioRad’s Experion electrophoresis station. A total of 750 ng of each biotin-labeled cRNA sample was hybridized to Illumina’s RaRef-12-v1 Expression BeadChips at 58°C overnight (17 h) according to the manufacturer’s instructions . The hybridized biotinylated cRNA was detected with 1ug/ml Cyanine3-streptavidine, (GE Healthcare Biosciences) and the Beadchips were scanned with Illumina BeadArray Reader (Factor=1, PMT=521, Filter=100%). Numerical results were extracted with Bead Studio v18.104.22.168 without any normalization or background subtraction.
Data of gastrin-regulated gene expression which compared treatment in the transient (1 h) versus sustained (14 h) mode is a time course reference design with a non-treated reference measured at each sampled time point. The data was normalised by loess adjustment within time points and average quantile normalised between time points. The data was analysed using the Limma (ver. 3.12.1) Bioconductor package . Replicated treated and untreated (reference) samples were compared using two separate linear models. The first model covers each time point from 0–4 h and the second and main model encompasses the remaining samples (6–14 h). Both models were tested for overall significance across time using a moderated F-test. Genes with a FDR adjusted p-value < 0.05 were taken as significant. The ratio between treated and untreated (reference) sample at each time point was averaged between biological replicates, and 403 significant genes were visualised in a heat map (Figure 1B) where genes (rows) are ordered by hierarchical clustering using Euclidean distance and Ward agglomeration. Mapping between Illumina probe identifiers and gene symbol was fetched from official Illumina annotation files (RatRef-12_V1_0_R5_11222119).
Database submission of microarray data
In addition to the time series experiment comparing treatment in transient (1 h) versus sustained (14 h) mode, we used data from two other time series experiments to analyse temporal time profiles of early and late gastrin responsive genes and assess gene expression in the presence and absence of the protein synthesis inhibitor cyclohexamide (CHX) to distinguish primary and secondary genes. The microarray data were prepared according to minimum information about a microarray experiment (MIAME) recommendations  and deposited in the Array Express . Detailed information about the microarray designs and raw data files from the experiments are accessible by use of these accession numbers: GSE32869, and E-MTAB-1268 (Illumina platform).
Enrichment analysis of pathways, GO-annotations and transcription factor binding site (TFBS)
Canonical pathways and networks most significantly associated by the genes differentially expressed in transienly versus sustained gastrin stimulated cells were determined by GeneGo MetaCoreTM enrichment analysis with the general p-value threshold p<0.05 for the data inputs. Molecular functions were defined based on data from MetaCore tool of GeneGo package i.e., Network objects (object types), Gene Ontology Annotation (UniProt-GOA) Database GOA  and literature. Transcription factor binding site (TFBS) enrichment was assessed for human orthologues of the rat genes by using DAVID, Database for Annotation, Visualization and Integrated Discovery, version 6.7 [72, 73].
Real time qRT-PCR
cDNA synthesis was performed using Transcriptor First Strand cDNA Synthesis Kit from Roche Diagnostics GmbH (Mannheim, Germany), as described earlier . After the cDNA synthesis reaction, the samples were diluted 1:3 in sterile water. Quantitative RT-PCR (qRT-PCR) was performed using StepOnePlus, Real-Time PCR System from Applied biosystems and B-R SYBR Green SuperMix for iQ (Quanta Bioscience, Gaithersburg, MD). Each reaction contained SYBR Green SuperMix (2×) (12.5 μl), sense primer (2.5 μl of 3 μM), antisense primer (2.5 μl of 3 μM), cDNA template (2.5 μl), and sterile water (5.0 μl). Quantitative PCR thermal cycling program: 40 thermal cycles of 15 s at 95°C, 30 s at 60°C, and 40 s at 72°C. A dissociation curve was made to confirm primer specificity. Primers: (Rattus Norvegicus): Clu, sense: 5′-GCTCCATAGCCCAGCTTTAC-3′, antisense: 5′-ACTTCTCACACTGGCCCTTC-3′ ; B2m: sense: 5′-CGAGACCGATGTATATGCTTGC-3′, antisense: 5′-GTCCAGATGATTCAGAGCTCCA-3′ ; Atf4, sense: 5′-GTTGGTCAGTGCCTCAGACA-3′, antisense: 5′-CATTCGAAACAGAGCATCGA-3′ ; Junb, sense: 5′-AGCTAGCCTCCACGGAACT-3′, antisense: 5′-CTCCTGCTCCTCGGTGAC-3′; Herpud1, sense: 5′-TTGCACCTCGTGTGCAATGTGAGG-3′, antisense: 5′-ACTAGTGTTGTCCGGCTGCTCTGT-3′; Chac1, sense: 5′-CCTTCCACAGGGGCAGCGATAAGAT-3′, antisense: 5′-AACCTGGTATGCCACACCCCAAGTG-3′. Primers for Junb, Herpud1 and Chac1 were designed using Clone Manager 9 (Scientific & Educational Software, Cary, NC) or Primer-Blast . All samples were run in triplicates, and relative mRNA expression levels were quantified using the ΔΔCt-method  with B2m as reference gene.
Western blot analysis
Cells were washed twice in PBS, and harvested in 400 μl cold RIPA buffer (Thermo Scientific, Rockford, IL) containing 8 μl Halt™ Protease Inhibitor Single-Use Cocktail (100×) (Thermo Scientific) and 4 μl Halt™ Phosphatase Inhibitor Single-Use Cocktail (100×) (Thermo Scientific) per well. The lysate was homogenized using a syringe and needle (21 G), and cell debris removed by centrifugation (15 min, 14000 ×g, 4°C). The supernatant was stored at −80°C until further processing. The protein concentrations were measured using the Pierce® BCA Protein Assay Kit (Thermo Scientific) according to the manufacturer’s instructions. Gel electrophoresis and blotting were performed using the NuPAGE system from Invitrogen according to the manufacturer’s instructions; as previously described . The membranes were blocked in TBST (50 mM Tris·HCl pH 7.5, 150 mM NaCl, 0.1% Tween 20) with 5% bovine serum albumin (BSA) (Sigma-Aldrich) for 1 h at room temperature, incubated with the primary antibody diluted in TBST with 1% BSA overnight at 4°C, and incubated with the secondary HRP-conjugated antibody diluted in TBST with 1% BSA for 2 h at room temperature. For detection of α-tubulin (loading control), blocking, washing, and antibody incubation were performed with SNAP i.d. Protein Detection System (Millipore, Billerica, MA) according to the manufacturer’s instructions, using TBST with 1% BSA for blocking and antibody dilution. Binding of secondary antibodies was visualized by use of the SuperSignal West Femto Maximum Sensitive Substrate (Thermo Scientific) and Kodak Image Station 2000R (Kodak, Pittsburgh, PA).
The following antibodies were used: monoclonal (mouse) anti-human α-Tubulin (reacts with mouse, rat and human) (1:300; sc-5286, Santa Cruz Biotechnology), polyclonal horseradish peroxidase (HRP) conjugated (goat) anti-mouse (1:3000; PO447, DAKO), polyclonal (rabbit) phospho-p44/42 MAPK (Erk1/2) (Thr202/Tyr204) (1:1000; #9101, Cell Signaling), polyclonal (rabbit) p44/42 MAPK (Erk1/2) (1:1000; #9102, Cell Signaling), monoclonal (rabbit) JunB (1:1000; #3753 Cell Signaling), polyclonal (rabbit) phospho-AKT (Ser473) (1:1000; #9271, Cell Signaling) and polyclonal (rabbit) AKT (1:1000;#9272, Cell Signaling), HRP-conjugated goat anti-rabbit (1:1000; #7074, Cell Signaling).
Stable knockdown (KD) of JUNB by retrovirus-based RNAi
Eight 97-mer shRNA oligoes were designed to target Junb using the algorithm from . The oligoes were cloned into MSCV P2Gm FF retroviral transfer vector (Addgene) using XhoI and EcoRI restriction enzymes (New England Biolabs). Correct clones were verified by sequencing. Packaged virus was obtained by cotransfection of 293FT cells (ATCC) with 1μg MSCV P2Gm FF transfer vector (Addgene), 0.5 μg MLV gag-pol (Addgene), and 0.5 μg CMV-VSVg (Addgene) expression vectors using Fugene6 (Roche) according to manufacturer’s protocol and as previously described . Medium was replaced 24 h and 48 h after transfection, 293FT cell conditioned medium was collected, filtered through a 0.45μm filter, and applied to AR42J cells with 4 μg/ml polybrene (Sigma) . Infected AR42J cells were GFP positive and selected with 2 μg/ml puromycin (Sigma). The cell lines with most efficient knock down were chosen for further analysis.
Hairpin sequence: shJunb
AR42J cells (1.5 × 104) were seeded in white-walled 96-well plates (Perkin Elmer). The cells were incubated for one day before the growth medium was replenished with 150 μl serum-free DMEM with or without different doses of gastrin (1–10 nM). After 1 h, the medium was again replenished and gastrin added to cells treated in the sustained mode. Cells were further incubated for 72 h to induce apoptosis. Chemical inhibitors were added 30 min before the gastrin treatment. Caspase activity was measured using the Caspase-Glo 3/7 assay from Promega (Madison, WI) according to the manufacturer’s descriptions. Luminescence was measured using Wallac 1420 Victor3TM plate reader (Perkin Elmer). All the data were log-transformed before statistical analysis. For presentation, the data were transformed back to the original scale and plotted as means with 95% confidence interval (CI) as error bars. Student's t-test was used to evaluate statistical significance, and we performed a Bonferroni correction for multiple testing at 0.05 level in each analysis.
Linn-Karina M Selvik and Christina S. Fjeldbo should be regarded as joint First Authors.
Miller-Jensen K, Janes KA, Brugge JS, Lauffenburger DA: Common effector processing mediates cell-specific responses to stimuli. Nature. 2007, 448 (7153): 604-608. 10.1038/nature06001.
Glauser DA, Schlegel W: Mechanisms of transcriptional regulation underlying temporal integration of signals. Nucleic Acids Res. 2006, 34 (18): 5175-5183. 10.1093/nar/gkl654.
Glauser DA, Schlegel W: Sequential actions of ERK1/2 on the AP-1 transcription factor allow temporal integration of metabolic signals in pancreatic beta cells. FASEB J. 2007, 21 (12): 3240-3249. 10.1096/fj.06-7798com.
Zi Z, Feng Z, Chapnick DA, Dahl M, Deng D, Klipp E, Moustakas A, Liu X: Quantitative analysis of transient and sustained transforming growth factor-beta signaling dynamics. Mol Syst Biol. 2011, 7: 492-
Hao S, Baltimore D: The stability of mRNA influences the temporal order of the induction of genes encoding inflammatory molecules. Nature immunology. 2009, 10 (3): 281-288. 10.1038/ni.1699.
Nagashima T, Shimodaira H, Ide K, Nakakuki T, Tani Y, Takahashi K, Yumoto N, Hatakeyama M: Quantitative transcriptional control of ErbB receptor signaling undergoes graded to biphasic response for cell differentiation. J Biol Chem. 2007, 282 (6): 4045-4056.
Watson SA, Grabowska AM, El-Zaatari M, Takhar A: Gastrin - active participant or bystander in gastric carcinogenesis?. Nat Rev Cancer. 2006, 6 (12): 936-946. 10.1038/nrc2014.
Dockray GJ: Clinical endocrinology and metabolism. Gastrin. Best Pract Res Clin Endocrinol Metab. 2004, 18 (4): 555-568. 10.1016/j.beem.2004.07.003.
Dufresne M, Seva C, Fourmy D: Cholecystokinin and gastrin receptors. Physiol Rev. 2006, 86 (3): 805-847. 10.1152/physrev.00014.2005.
Burkitt MD, Varro A, Pritchard DM: Importance of gastrin in the pathogenesis and treatment of gastric tumors. World J Gastroenterol. 2009, 15 (1): 1-16. 10.3748/wjg.15.1.
Wiedemann T, Hofbaur S, Tegtmeyer N, Huber S, Sewald N, Wessler S, Backert S, Rieder G: Helicobacter pylori CagL dependent induction of gastrin expression via a novel alphavbeta5-integrin-integrin linked kinase signalling complex. Gut. 2012, 61 (7): 986-996. 10.1136/gutjnl-2011-300525.
Caplin M, Savage K, Khan K, Brett B, Rode J, Varro A, Dhillon A: Expression and processing of gastrin in pancreatic adenocarcinoma. Br J Surg. 2000, 87 (8): 1035-1040. 10.1046/j.1365-2168.2000.01488.x.
Goetze JP, Nielsen FC, Burcharth F, Rehfeld JF: Closing the gastrin loop in pancreatic carcinoma: coexpression of gastrin and its receptor in solid human pancreatic adenocarcinoma. Cancer. 2000, 88 (11): 2487-2494. 10.1002/1097-0142(20000601)88:11<2487::AID-CNCR9>3.0.CO;2-E.
Chao C, Hellmich MR: Gastrin, inflammation, and carcinogenesis. Curr Opin Endocrinol Diabetes Obes. 2010, 17 (1): 33-39.
Hetz C: The unfolded protein response: controlling cell fate decisions under ER stress and beyond. Nat Rev Mol Cell Biol. 2012, 13 (2): 89-102.
Walter P, Ron D: The unfolded protein response: from stress pathway to homeostatic regulation. Science. 2011, 334 (6059): 1081-1086. 10.1126/science.1209038.
Li X, Zhang K, Li Z: Unfolded protein response in cancer: the physician's perspective. J Hematol Oncol. 2011, 4: 8-10.1186/1756-8722-4-8.
Gorman AM, Healy SJ, Jager R, Samali A: Stress management at the ER: Regulators of ER stress-induced apoptosis. Pharmacol Ther. 2012, 134 (3): 306-316. 10.1016/j.pharmthera.2012.02.003.
Ma Y, Hendershot LM: Herp is dually regulated by both the endoplasmic reticulum stress-specific branch of the unfolded protein response and a branch that is shared with other cellular stress pathways. J Biol Chem. 2004, 279 (14): 13792-13799. 10.1074/jbc.M313724200.
Mungrue IN, Pagnon J, Kohannim O, Gargalovic PS, Lusis AJ: CHAC1/MGC4504 is a novel proapoptotic component of the unfolded protein response, downstream of the ATF4-ATF3-CHOP cascade. J Immunol. 2009, 182 (1): 466-476.
Belal C, Ameli NJ, El Kommos A, Bezalel S, Al'Khafaji AM, Mughal MR, Mattson MP, Kyriazis GA, Tyrberg B, Chan SL: The homocysteine-inducible endoplasmic reticulum (ER) stress protein Herp counteracts mutant alpha-synuclein-induced ER stress via the homeostatic regulation of ER-resident calcium release channel proteins. Hum Mol Genet. 2012, 21 (5): 963-977. 10.1093/hmg/ddr502.
Allagnat F, Fukaya M, Nogueira TC, Delaroche D, Welsh N, Marselli L, Marchetti P, Haefliger JA, Eizirik DL, Cardozo AK: C/EBP homologous protein contributes to cytokine-induced pro-inflammatory responses and apoptosis in beta-cells. Cell Death Differ. 2012, 19: 1836-1846. 10.1038/cdd.2012.67.
Fjeldbo CS, Bakke I, Erlandsen SE, Holmseth J, Laegreid A, Sandvik AK, Thommesen L, Bruland T: Gastrin upregulates the prosurvival factor secretory clusterin in adenocarcinoma cells and in oxyntic mucosa of hypergastrinemic rats. Am J Physiol Gastrointest Liver Physiol. 2012, 302 (1): G21-33. 10.1152/ajpgi.00197.2011.
Todisco A, Ramamoorthy S, Witham T, Pausawasdi N, Srinivasan S, Dickinson CJ, Askari FK, Krametter D: Molecular mechanisms for the antiapoptotic action of gastrin. Am J Physiol Gastrointest Liver Physiol. 2001, 280 (2): G298-307.
Stepan V, Ramamoorthy S, Pausawasdi N, Logsdon CD, Askari FK, Todisco A: Role of small GTP binding proteins in the growth-promoting and antiapoptotic actions of gastrin. Am J Physiol Gastrointest Liver Physiol. 2004, 287 (3): G715-725. 10.1152/ajpgi.00169.2003.
Ramamoorthy S, Stepan V, Todisco A: Intracellular mechanisms mediating the anti-apoptotic action of gastrin. Biochem Biophys Res Commun. 2004, 323 (1): 44-48. 10.1016/j.bbrc.2004.08.059.
Zoubeidi A, Gleave M: Small heat shock proteins in cancer therapy and prognosis. Int J Biochem Cell Biol. 2012, 44 (10): 1646-56. 10.1016/j.biocel.2012.04.010.
Li N, Zoubeidi A, Beraldi E, Gleave ME: GRP78 regulates clusterin stability, retrotranslocation and mitochondrial localization under ER stress in prostate cancer. Oncogene. 2012, 32: 1933-1942.
Michels J, Johnson PW, Packham G: Mcl-1. Int J Biochem Cell Biol. 2005, 37 (2): 267-271. 10.1016/j.biocel.2004.04.007.
Pritchard DM, Berry D, Przemeck SM, Campbell F, Edwards SW, Varro A: Gastrin increases mcl-1 expression in type I gastric carcinoid tumors and a gastric epithelial cell line that expresses the CCK-2 receptor. Am J Physiol Gastrointest Liver Physiol. 2008, 295 (4): G798-805. 10.1152/ajpgi.00015.2008.
Eckenrode EF, Yang J, Velmurugan GV, Foskett JK, White C: Apoptosis protection by Mcl-1 and Bcl-2 modulation of inositol 1,4,5-trisphosphate receptor-dependent Ca2+ signaling. J Biol Chem. 2010, 285 (18): 13678-13684. 10.1074/jbc.M109.096040.
Reeves MA, Bellinger FP, Berry MJ: The neuroprotective functions of selenoprotein M and its role in cytosolic calcium regulation. Antioxid Redox Signal. 2010, 12 (7): 809-818. 10.1089/ars.2009.2883.
Choi HJ, Hwang S, Lee SH, Lee YR, Shin J, Park KS, Cho YM: Genome-wide identification of palmitate-regulated immediate early genes and target genes in pancreatic beta-cells reveals a central role of NF-kappaB. Mol Biol Rep. 2012, 39 (6): 6781-6789. 10.1007/s11033-012-1503-5.
Amit I, Citri A, Shay T, Lu Y, Katz M, Zhang F, Tarcic G, Siwak D, Lahad J, Jacob-Hirsch J, et al: A module of negative feedback regulators defines growth factor signaling. Nat Genet. 2007, 39 (4): 503-512. 10.1038/ng1987.
Glauser DA, Brun T, Gauthier BR, Schlegel W: Transcriptional response of pancreatic beta cells to metabolic stimulation: large scale identification of immediate-early and secondary response genes. BMC Mol Biol. 2007, 8: 54-10.1186/1471-2199-8-54.
Cullingford TE, Markou T, Fuller SJ, Giraldo A, Pikkarainen S, Zoumpoulidou G, Alsafi A, Ekere C, Kemp TJ, Dennis JL, et al: Temporal regulation of expression of immediate early and second phase transcripts by endothelin-1 in cardiomyocytes. Genome Biol. 2008, 9 (2): R32-10.1186/gb-2008-9-2-r32.
Gene Ontology Annotation (UniProt-GOA) Database. http://www.ebi.ac.uk/GOA/,
Fowler T, Sen R, Roy AL: Regulation of primary response genes. Mol Cell. 2011, 44 (3): 348-360. 10.1016/j.molcel.2011.09.014.
Tullai JW, Schaffer ME, Mullenbrock S, Sholder G, Kasif S, Cooper GM: Immediate-early and delayed primary response genes are distinct in function and genomic architecture. J Biol Chem. 2007, 282 (33): 23981-23995. 10.1074/jbc.M702044200.
Koblinski JE, Ahram M, Sloane BF: Unraveling the role of proteases in cancer. Clin Chim Acta. 2000, 291 (2): 113-135. 10.1016/S0009-8981(99)00224-7.
Jiang G, Cao F, Ren G, Gao D, Bhakta V, Zhang Y, Cao H, Dong Z, Zang W, Zhang S, et al: PRSS3 promotes tumour growth and metastasis of human pancreatic cancer. Gut. 2010, 59 (11): 1535-1544. 10.1136/gut.2009.200105.
Ellrichmann M, Ritter PR, Schrader H, Schmidt WE, Meier JJ, Schmitz F: Gastrin stimulates the VEGF-A promotor in a human colon cancer cell line. Regul Pept. 2010, 165 (2–3): 146-150.
Willard MD, Lajiness ME, Wulur IH, Feng B, Swearingen ML, Uhlik MT, Kinzler KW, Velculescu VE, Sjoblom T, Markowitz SD, et al: Somatic Mutations in CCK2R Alter Receptor Activity that Promote Oncogenic Phenotypes. Mol Cancer Res. 2012, 10 (6): 739-749. 10.1158/1541-7786.MCR-11-0483.
Peixoto P, Castronovo V, Matheus N, Polese C, Peulen O, Gonzalez A, Boxus M, Verdin E, Thiry M, Dequiedt F, et al: HDAC5 is required for maintenance of pericentric heterochromatin, and controls cell-cycle progression and survival of human cancer cells. Cell Death Differ. 2012, 19 (7): 1239-1252. 10.1038/cdd.2012.3.
Ellis LM, Hicklin DJ: VEGF-targeted therapy: mechanisms of anti-tumour activity. Nat Rev Cancer. 2008, 8 (8): 579-591. 10.1038/nrc2403.
Modlin IM, Kidd M, Latich I, Zikusoka MN, Eick GN, Mane SM, Camp RL: Genetic differentiation of appendiceal tumor malignancy: a guide for the perplexed. Ann Surg. 2006, 244 (1): 52-60. 10.1097/01.sla.0000217617.06782.d5.
Mullenbrock S, Shah J, Cooper GM: Global expression analysis identified a preferentially nerve growth factor-induced transcriptional program regulated by sustained mitogen-activated protein kinase/extracellular signal-regulated kinase (ERK) and AP-1 protein activation during PC12 cell differentiation. Journal of Biological Chemistry. 2011, 286 (52): 45131-45145. 10.1074/jbc.M111.274076.
Ebisuya M, Kondoh K, Nishida E: The duration, magnitude and compartmentalization of ERK MAP kinase activity: mechanisms for providing signaling specificity. Journal of Cell Science. 2005, 118 (Pt 14): 2997-3002.
Sasagawa S, Ozaki Y, Fujita K, Kuroda S: Prediction and validation of the distinct dynamics of transient and sustained ERK activation. Nat Cell Biol. 2005, 7 (4): 365-373. 10.1038/ncb1233.
Santos SD, Verveer PJ, Bastiaens PI: Growth factor-induced MAPK network topology shapes Erk response determining PC-12 cell fate. Nat Cell Biol. 2007, 9 (3): 324-330. 10.1038/ncb1543.
Katz M, Amit I, Yarden Y: Regulation of MAPKs by growth factors and receptor tyrosine kinases. Biochim Biophys Acta. 2007, 1773 (8): 1161-1176. 10.1016/j.bbamcr.2007.01.002.
Steigedal TS, Prestvik WS, Selvik LK, Fjeldbo CS, Bruland T, Laegreid A, Thommesen L: Gastrin-induced proliferation involves MEK partner 1 (MP1). In Vitro Cell Dev Biol Anim. 2013, 49 (3): 162-169. 10.1007/s11626-013-9588-2.
Gunther R, Carstens OC, Schmidt WE, Folsch UR: Transient agonist-induced regulation of the cholecystokinin-A and cholecystokinin-B receptor mRNA levels in rat pancreatic acinar AR42J cells. Pancreatology. 2003, 3 (1): 47-54. 10.1159/000069142.
Ashurst HL, Varro A, Dimaline R: Regulation of mammalian gastrin/CCK receptor (CCK2R) expression in vitro and in vivo. Exp Physiol. 2008, 93 (2): 223-236.
Lopez-Bergami P, Lau E, Ronai Z: Emerging roles of ATF2 and the dynamic AP1 network in cancer. Nat Rev Cancer. 2010, 10 (1): 65-76. 10.1038/nrc2681.
Vesely PW, Staber PB, Hoefler G, Kenner L: Translational regulation mechanisms of AP-1 proteins. Mutat Res. 2009, 682 (1): 7-12. 10.1016/j.mrrev.2009.01.001.
Shaulian E: AP-1–The Jun proteins: Oncogenes or tumor suppressors in disguise?. Cell Signal. 2010, 22 (6): 894-899. 10.1016/j.cellsig.2009.12.008.
Ameri K, Harris AL: Activating transcription factor 4. Int J Biochem Cell Biol. 2008, 40 (1): 14-21. 10.1016/j.biocel.2007.01.020.
Hu J, Dang N, Menu E, De Bryune E, Xu D, Van Camp B, Van Valckenborgh E, Vanderkerken K: Activation of ATF4 mediates unwanted Mcl-1 accumulation by proteasome inhibition. Blood. 2012, 119 (3): 826-837. 10.1182/blood-2011-07-366492.
Gurzov EN, Ortis F, Bakiri L, Wagner EF, Eizirik DL: JunB inhibits ER stress and apoptosis in pancreatic beta cells. PLoS One. 2008, 3 (8): e3030-10.1371/journal.pone.0003030.
Gurzov EN, Barthson J, Marhfour I, Ortis F, Naamane N, Igoillo-Esteve M, Gysemans C, Mathieu C, Kitajima S, Marchetti P, et al: Pancreatic beta-cells activate a JunB/ATF3-dependent survival pathway during inflammation. Oncogene. 2012, 31 (13): 1723-1732. 10.1038/onc.2011.353.
Androulakis IP, Yang E, Almon RR: Analysis of time-series gene expression data: methods, challenges, and opportunities. Annu Rev Biomed Eng. 2007, 9: 205-228. 10.1146/annurev.bioeng.9.060906.151904.
Gonzalez A, Santofimia-Castaño P, Salido GM: Culture of pancreatic AR42J cell for use as a model for acinar cell function. The Pancreapedia: Exocrine Pancreas Knowledge Base. http://www.pancreapedia.org/?q=node/769,
Stepan VM, Dickinson CJ, del Valle J, Matsushima M, Todisco A: Cell type-specific requirement of the MAPK pathway for the growth factor action of gastrin. Am J Physiol. 1999, 276 (6 Pt 1): G1363-1372.
Hofsli E, Thommesen L, Norsett K, Falkmer S, Syversen U, Sandvik A, Laegreid A: Expression of chromogranin A and somatostatin receptors in pancreatic AR42J cells. Mol Cell Endocrinol. 2002, 194 (1–2): 165-173.
Daulhac L, Kowalski-Chauvel A, Pradayrol L, Vaysse N, Seva C: Gastrin stimulates the formation of a p60Src/p125FAK complex upstream of the phosphatidylinositol 3-kinase signaling pathway. FEBS Lett. 1999, 445 (2–3): 251-255.
httpIllumina RatRef-12 Expression BeadChip. http://support.illumina.com/array/array_kits/ratref-12_expression_beadchip_kit.ilmn,
Smyth GK: Limma: linear models for microarray data. Bioinformatics and Computational Biology Solutions using R and Bioconductor. Edited by: Gentleman R, Carey V, Dudoit S, Irizarry R, Huber W. 2005, New York: Springer, 397-420.
Brazma A, Hingamp P, Quackenbush J, Sherlock G, Spellman P, Stoeckert C, Aach J, Ansorge W, Ball CA, Causton HC, et al: Minimum information about a microarray experiment (MIAME)-toward standards for microarray data. Nat Genet. 2001, 29 (4): 365-371. 10.1038/ng1201-365.
Parkinson H, Sarkans U, Shojatalab M, Abeygunawardena N, Contrino S, Coulson R, Farne A, Lara GG, Holloway E, Kapushesky M, et al: ArrayExpress--a public repository for microarray gene expression data at the EBI. Nucleic Acids Res. 2005, 33 (Database issue): D553-555.
da Huang W, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009, 4 (1): 44-57.
da Huang W, Sherman BT, Lempicki RA: Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009, 37 (1): 1-13. 10.1093/nar/gkn923.
Cardozo AK, Ortis F, Storling J, Feng YM, Rasschaert J, Tonnesen M, van Eylen F, Mandrup-Poulsen T, Herchuelz A, Eizirik DL: Cytokines downregulate the sarcoendoplasmic reticulum pump Ca2+ ATPase 2b and deplete endoplasmic reticulum Ca2+, leading to induction of endoplasmic reticulum stress in pancreatic beta-cells. Diabetes. 2005, 54 (2): 452-461.
Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) Method. Methods. 2001, 25 (4): 402-408. 10.1006/meth.2001.1262.
RNAi Central shRNA retriever. http://cancan.cshl.edu/RNAi_central/RNAi.cgi?type=shRNA,
Stern P, Astrof S, Erkeland SJ, Schustak J, Sharp PA, Hynes RO: A system for Cre-regulated RNA interference in vivo. Proc Natl Acad Sci U S A. 2008, 105 (37): 13895-13900. 10.1073/pnas.0806907105.
This work was supported by The Liaison Committee between the Central Norway Regional Health Authority (RHA) and the Norwegian University of Science and Technology (NTNU); Sør-Trøndelag University College (HiST); The Norwegian Research Council; The microarray and bioinformatics services were provided by the Genomics Core Facility, Norwegian University of Science and Technology, and NMC - a national technology platform supported by the functional genomics program (FUGE) of the Research Council of Norway. We thank Anne Kristensen and Dr. Kirsti Berg for skilful technical assistance.
The authors declare that they have no competing interests.
LKMS carried out cell and molecular studies and helped to draft the manuscript. CSF participated in design and interpretation of microarray experiments, carried out cell and molecular studies and helped to draft the manuscript. AF participated in design of microarray experiments and performed statistical analysis and bioinformatics. TSS carried out cell and molecular studies and participated in interpretation of microarray experiments. KM carried out cell and molecular studies and participated in interpretation of microarray experiments. EA participated in design of microarray experiments and performed statistical analysis and bioinformatics. BD participated in design and analysis of microarray experiments and carried out initial microarray experiments. ML participated in design and statistical analysis of initial microarray experiments. ST carried out cell and molecular studies and participated in enrichment analysis. VB participated in design and analysis of microarray experiments and participated in enrichment analysis. AL raised funding, participated in study design and coordination, participated in interpretation of results, and helped to draft the manuscript. LT raised funding, participated in study design and coordination, participated in interpretation of results, and helped to draft the manuscript. TB participated in study design and coordination, participated in design, experimental work, analysis and interpretation of microarray experiments, carried out cell and molecular studies and drafted the manuscript. All authors read and approved the final manuscript.
Linn-Karina M Selvik, Christina S Fjeldbo contributed equally to this work.
Electronic supplementary material
Additional file 2: The file includes results from an initial cDNA microarray time series experiment (Accession number: E-MTAB-123). Figure S1A: Schematic representation of stimulation protocol. Figure S1B: Graphical representation of experimental design of the two-colour cDNA microarray hybridizations. Figure S2: Number of differentially expressed genes. Figure S3: Analysis of microarray time series by dimension reduction methods. (PDF 496 KB)
Additional file 3: Table S2: Enrichment analysis of 259 genes significantly lower expressed in transiently versus sustained gastrin treated cells. Table S3: Enrichment analysis of 144 genes significantly higher expressed in transiently versus sustained gastrin treated cells. All 403 genes are shown in the heat map in Figure 1B and Additional file 1. (PDF 115 KB)
Additional file 4: Table S4: The file lists a subset of 181 gastrin-induced probes (177 unique IDs) with differing expression patterns in transient versus sustained mode which were used to further characterize the temporal profiles including information about gene expression in the presence and absence of a protein synthesis inhibitor (CHX) as well as molecular functions. The data are extracted from independent time series experiments in Accession numbers GSE32869 and E-MTAB-1268. (XLSX 26 KB)
Additional file 5: Temporal profiles of gastrin-induced genes with differing expression patterns in transient versus sustained mode described in Figure 5 and Additional file 4 : Table S4. Time profiles as log2 expression data from three independent time series experiments are shown for each individual gene. Upper panels: transiently or sustained treated cells; Middle panels: sustained treated or untreated control cells. Lower panels: gene expression in the presence and absence of the protein synthesis inhibitor CHX. (PDF 172 KB)
Additional file 6: Table S5: Identification of overrepresented transcription factor binding sites (TFBS) among the subset of upregulated genes lower expressed in transiently gastrin treated cells (described in Figure 5 and Additional file 4: Table S4). TFBS enrichment was assessed for human orthologues of the rat genes by using DAVID version 6.7 (Bonferroni corrected p-values). (XLSX 13 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Selvik, L.M., Fjeldbo, C.S., Flatberg, A. et al. The duration of gastrin treatment affects global gene expression and molecular responses involved in ER stress and anti-apoptosis. BMC Genomics 14, 429 (2013) doi:10.1186/1471-2164-14-429
- Temporal Profiles
- Primary-and Secondary Genes
- UPR/ER Stress