Carcinogenesis Advance Access originally published online on August 19, 2005
Carcinogenesis 2006 27(2):319-327; doi:10.1093/carcin/bgi211
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Carcinogenesis vol.27 no.2 © Oxford University Press 2005; all rights reserved.
Low pH induces co-ordinate regulation of gene expression in oesophageal cells
1 Department of Clinical Medicine and 2 Department of Surgery, Institute of Molecular Medicine, Trinity Centre for Health Sciences, St James's Hospital, Dublin 8, Ireland and 3 Department of Pharmacology, Centre for Molecular Medicine, Conway Institute of Biomolecular and Biomedical Research, University College Dublin, Dublin 4, Ireland
* To whom correspondence should be addressed. Tel: +353 1 6082100; Fax: +353 1 6083503; Email: spduggan{at}tcd.ie or kellehdp{at}tcd.ie
| Abstract |
|---|
|
|
|---|
The development of gastro-oesophageal reflux disease (GORD) is known to be a causative risk factor in the evolution of adenocarcinoma of the oesophagus. The major component of this reflux is gastric acid. However, the impact of low pH on gene expression has not been extensively studied in oesophageal cells. This study utilizes a transcriptomic and bioinformatic approach to assess regulation of gene expression in response to low pH. In more detail, oesophageal adenocarcinoma cell lines were exposed to a range of pH environments. Affymetrix microarrays were used for gene-expression analysis and results were validated using cycle limitation and real-time RTPCR analysis, as well as northern and western blotting. Comparative promoter transcription factor binding site (TFBS) analysis (MatInspector) of hierarchically clustered gene-expression data was employed to identify the elements which may co-ordinately regulate individual gene clusters. Initial experiments demonstrated maximal induction of EGR1 gene expression at pH 6.5. Subsequent array experimentation revealed significant induction of gene expression from such functional categories as DNA damage response (EGR1-4, ATF3) and cell-cycle control (GADD34, GADD45, p57). Changes in expression of EGR1, EGR3, ATF3, MKP-1, FOSB, CTGF and CYR61 were verified in separate experiments and in a variety of oesophageal cell lines. TFBS analysis of promoters identified transcription factors that may co-ordinately regulate gene-expression clusters, Cluster 1: Oct-1, AP4R; Cluster 2: NF-kB, EGRF; Cluster 3: IKRS, AP-1F. Low pH has the ability to induce genes and pathways which can provide an environment suitable for the progression of malignancy. Further functional analysis of the genes and clusters identified in this low pH study is likely to lead to new insights into the pathogenesis and therapeutics of GORD and oesophageal cancer.
Abbreviations: ATF3, activating transcription factor 3; BO, Barrett's oesophagus; BEBM, bronchial epithelial basal media; DAF, decay accelerating factor; EGRF (early growth response factor; EGR1, early growth response gene 1; GORD, gastro-oesophageal reflux disease; IL-8, interleukin-8; NF-KB, nuclear factor-kappa-B; NR4A2, nuclear receptor 4A2; PTGS2, prostaglandin synthase-2; TFBS, transcription factor binding site
| Introduction |
|---|
|
|
|---|
Over the past few decades, Western populations have witnessed a 10% annual increase in the incidence of oesophageal adenocarcinoma, whereas the incidence of oesophageal squamous carcinoma remains unchanged (13). Gastro-oesophageal acid reflux has been identified in epidemiological studies as a major risk factor in the pathogenesis of oesophageal adenocarcinoma (4). Oesophageal adenocarcinoma commonly arises from the presence of a premalignant lesion known as Barrett's oesophagus (BO). This is a common condition in which normal squamous epithelium is replaced by a metaplastic intestine-like epithelium containing goblet cells (2,5). Although the exact model for the development of BO has not been fully elucidated, the current paradigm indicates that it most commonly occurs in the setting of gastro-oesophageal reflux disease (GORD), in conjunction with as yet unidentified genetic and environmental factors (2,69).
The GORD refluxate consists of many different agents, such as stomach acid, bile acids (depending on the form of GORD), pepsins and gastrins (10), as well as metabolized oxidative products from the diet. As gastric acid is the major constituent of the reluxate and is postulated to be the cause of the initial inflammatory disease or oesophagitis, we focused our investigation on its modulatory role (1113).
Previous work from our laboratory has demonstrated the ability of acidic or low pH conditions to induce pro-inflammatory signalling pathways, such as nuclear factor-kappaB (NF-KB) transcription factor activation (14,15). This activation is common in oesophageal adenocarcinoma where it inversely correlates with response to neo-adjuvant chemotherapy (14). It was of interest to expand the focus of this work to identify other transcriptional targets of low pH to provide further avenues for treatment. The early growth response gene 1 (EGR1) was an ideal candidate owing to its involvement in other inflammatory conditions and activation by carcinogens.
The EGR gene family of transcription factors encode nuclear proteins with zinc finger domains resembling those of the well characterized Sp1 transcription factor, which has been shown to have increased DNA binding activity under low pH conditions (1618). There is increasing evidence for the involvement of EGR1 in the regulation of immune response both in immune cells and epithelial cells (19,20). EGR1 is rapidly and transiently induced by many extracellular signals, such as growth factors, cytokines and phorbol esters, in a variety of cell types including fibroblasts, epithelial and endothelial cells (21,22). In addition, EGR1 is induced by diverse types of DNA damaging agents and mild oxidative stress. We hypothesized that EGR1 may be induced by low pH in the oesophageal cell lines SKGT4 and OE-33 and may lead to the identification of other related genes involved in low pH signalling. A novel bioinformatic approach for the comparative analysis of transcription factor binding sites (TFBSs) was used to analyse data derived from transcriptomic experiments. This identified distinct subsets of binding elements in the promoter regions that may co-ordinately control the transcription of particular gene clusters. In particular, the NF-KB and early growth response factor (EGRF) binding elements were identified as potentially important mediators in low pH-based gene induction.
| Materials and methods |
|---|
|
|
|---|
Cells
OE-33 cells were obtained from the European Collection of Cell Cultures (Salisbury, UK) and maintained in culture with RPMI 1640 medium supplemented with 10% foetal calf serum, 100 U/ml penicillin, 100 µg/ml streptomycin and 2 mM L-glutamine (Gibco BRL, Grand Island, NY) at 37°C in an atmosphere containing 5% CO2. The SKGT4 cell line was obtained from Schrump et al. (National Cancer Institute, MD) (23) and cultured as above. HET-1A cells were purchased from the American Type Culture Collection (ATCC, Manassas, VA, USA) and cultured in bronchial epithelial basal media (BEBM) media with growth supplements (Cambrex Bioscience, Berkshire, UK). For western blot analysis, cells were seeded in 6-well plates at a density of 5 x 105 cells/ml medium and grown to 80% confluence prior to experimentation. Cell viability was confirmed using acridine orange/ethidium bromide (AO/EB) and propidium iodide (PI) fluorescence staining (Sigma, Poole, UK). For gene-expression studies, 80 cm2 culture flasks were seeded with
5 x 106 cells in order to obtain sufficient RNA for northern blotting, RTPCR and DNA microarray analyses. Acidic media was made by titrating RPMI (without serum at a temperature of 37°C) with hydrochloric acid (0.4 M) to the required pH. Before acid treatment, SKGT4 cells were allowed to rest in serum-free media for 12 h before experimentation. No additive pH buffers were used; however, pH was monitored throughout the experiment and adjusted to within the range required. Resting control cells were mock titrated with equal volumes of water where required for each time point, to control for osmotic effects. RNA samples were taken after pH monitoring from the respective time points up to 4 h and at the same time as the appropriate time point controls. Experiments involving the HET-1A cell line were performed as above using BEBM media without the use of growth supplements.
Western blot analysis
Cells were lysed in 20 mM TrisHCL (pH 7.5), 1% SDS, 150 mM NaCl, 1mM EGTA, 1mM EDTA, 0.5 mM PMSF and 10 µg/ml leupeptin (Sigma), followed by sonication and boiling. Protein extracts were resolved by SDSPAGE and transferred to PVDF membrane (Millipore, Bedford, MA). Blots were blocked with 5% (w/v) dried skim milk in phosphate buffered saline (PBS) for 1 h at room temperature and then incubated for 1 h at room temperature with the primary polyclonal antibody to appropriate protein at a variable dilution in PBS/powdered milk. The immunoblots were then washed and incubated with an anti-rabbit, horseradish peroxidase-conjugated secondary antibody (DAKO, Bucks, UK). Immunodetection was accomplished by enhanced chemiluminescence.
RNA isolation and northern blot analysis
Total RNA was isolated using nucleospin extraction columns (Clontech, Palo Alto, CA) including a DNase step, while maintaining all samples on ice throughout the process. RNA samples (5 µg) were run on 1.5% formaldehyde gels and integrity checked by UV transillumination. The RNA was transferred from formaldehyde gels to nylon membranes (Amersham Biosciences, Bucks, UK) by overnight capillary transfer. The membranes were cross-linked at 80°C for 2 h in an oven. Blots were checked for linear transfer by a methylene blue staining protocol (Sigma). Radiolabelled probes for northern blot analysis were prepared by PCR amplification of an appropriate region of the gene followed by incorporation of the labelled [32P]dCTP nucleotide by random priming (Amersham, Aylesbury, UK). Templates for the probes were produced using the following primers: EGR-1 (forward 5'-CAGCTTGGTCAGTGGCCTAG-3' and reverse 5' AATGTCAGTGTTCGGCGTG-3' from NCBI AJ24325) and GAPDH (forward 5'-CAGAACATCATCCCTGCCTCTAC-3' and reverse 5'-TTGAAGTCAGAGGAGACCACCTG-3' from NCBI J04038
[GenBank]
), synthesized by Invitrogen (Paisley, UK). Membranes were probed in RapidHyb buffer (Amersham Biosciences) at 65°C overnight followed by washing with high stringency salt buffer.
Semi-quantitative RTPCR analysis
Total RNA was isolated as above and was transcribed to cDNA using the RETROscript RTPCR kit (Ambion, Huntingdon, UK). RNA (1 µg) was reverse transcribed in a total volume of 20 µl. Of this cDNA reaction, 5 µl was subjected to PCR amplification with SuperTaq (Ambion) using appropriate primers. Primers used for RTPCR analysis were as detailed in Table I. Assessment of RTPCR-derived products was performed via agarose gel electrophoresis at exponential points within respective reaction cycles.
|
Real-time RTPCR analysis
Extraction of RNA and reverse transcription was performed as above. An aliquot of 1 µl of the reverse transcription reaction was used in the real-time PCR reactions (20 µl final volume) and performed in an ABI Prism 7700 thermocycler (Applied Biosystems, Foster City, CA). Primers for the genes (ATF3, EGR3, CYR61, DUSP5, MAN2C, MNT, C8FW) were obtained from the ABI Assay on Demand website (http://home.appliedbiosystems.com/). Fold inductions were calculated using the cycle threshold

Ct method as previously described (24,25).
DNA microarray analysis
An aliquot of 10 µg of total RNA from treated SKGT4 cells was reverse transcribed into single-stranded cDNA using the SuperScript Choice kit (Invitrogen, Paisley, UK). For this purpose, an oligo-dT primer containing a T7 RNA polymerase promoter (Genset) was utilized. Following double-stranded cDNA synthesis, biotin-labelled cRNA was generated by in vitro transcription using the BioArray RNA labelling kit (Enzo Diagnostics, Farmingdale, N.Y., USA). These complex cRNA targets, which are representative of the transcriptome of a particular sample, were hybridized against U133A arrays (Affymetrix, San Jose, CA). Detection was accomplished via a streptavidin-labelled fluorochrome (phycoerythrin) and laser scanning. The raw microarray data have been submitted to Gene Expression Omnibus for public access and may be obtained using the accession code GSE2144
[NCBI GEO]
.
Inter-array comparison of gene-expression data
Normalization of data and inter-array comparisons of gene-expression profiles was carried out using Microarray Analysis Suite (MAS) version 5.0 (Affymetrix), together with Microsoft Access. In more detail, DNA microarray experiments were analysed using an approach based on the MannWhitney pairwise comparison test (26). To identify differentially expressed genes between any two samples, pairwise comparisons were performed using MAS. Lists of altered transcripts from different pairwise comparisons were sorted via Microsoft Access.
Cluster analysis
Detailed cluster analysis was performed using GeneSpring 4.0. Data were normalized initially by transforming the values <0.0010.001. Each measurement was divided by the 50th percentile of all measurements in that sample. The percentile was calculated with all raw measurements >0 and using all genes not marked as absent. Each gene-expression value was divided by the median of its measurement in all samples. If the median of the raw values was <0.01, then each measurement for that gene was divided by 0.01. Significantly altered genes were identified by viewing the data in the form of an ordered list of significance from which the top 100 were chosen for analysis. Guided analysis in the form of similarity searches were performed, followed by k-means clustering. These lists were then compared with more unguided methods, such as hierarchical clustering, before attempting to identify individual promoter clusters.
Bioinformatic analysis of promoter regions
Promoter regions for genes residing within clusters of interest were identified and retrieved via the use of Genome Browser housed at University of California Santa Cruz (http://genome.ucsc.edu/). Retrieved promoter regions (1500 bp upstream of the predicted start site) from collated clusters of genes were comparatively analysed using the TFBS analysis tool MatInspector (http://www.genomatix.de/) (27). The resultant list of TFBSs identified for each cluster were filtered against each other and against a cluster of randomly selected non-changing genes to give a list of the most significant TFBSs for each cluster.
| Results |
|---|
|
|
|---|
Induction of EGR1 mRNA by low pH in oesophageal cell lines
EGR1 is an early growth response gene whose expression can be induced by many different stimuli including cytokines, growth factors, DNA damaging agents and oxidative stress. To examine EGR1 expression under conditions of low pH, SKGT-4 and OE-33 cells were treated with growth media whose pH had been altered with 0.4 M hydrochloric acid. Figure 1A shows western blotting data for the induction of the 80 kDa EGR1 protein under conditions of varying pH for 1 h in both SKGT4 and OE-33 oesophageal cell lines. The pH ranges between 7 and 6, previously demonstrated to induce NF-KB (14), can clearly be seen to induce EGR1 protein expression. Induction of EGR1 mRNA was observed between pH 6.8 and 6.2 in both cell lines when the pH range was focused (Figure 1B). At pH 6.5 (chosen to reflect gene induction midpoint), maximal induction of EGR1 was seen after 1 h (Figure 1C), which gradually returned to resting levels after 3 h. This is mostly likely reflective of the increase in pH within the media over this time period. A system of pH monitoring was then used to maintain the pH within the range observed for EGR1 gene induction (Figure 1C); in this case, the level of EGR1 mRNA expression was still increased up to 4 h post-initial treatment.
|
Identification of transcriptomic responses to conditions of low pH
To further examine the effect of low pH conditions on a transcriptomic level, we employed DNA microarray-based gene-expression profile analysis. For this purpose, SKGT-4 cells were treated with low pH 6.5, together with pH monitoring, and analysed using Affymetrix U133A arrays, with
20 000 transcripts being examined in this respect. A range of genes were identified as being consistently differentially expressed in conditions of low pH as compared with the resting state, with the following number of transcripts being affected over time: 116 transcripts for 30 min, 480 transcripts for 120 min, 115 transcripts for 180 min and 741 transcripts for 240 min. Figure 2 details those genes whose expression was upregulated across all time periods in comparison with the resting state. Of note, EGR1 was identified from the DNA microarray studies as being one of the significantly altered genes, providing initial validation for the methodology employed. Moreover, additional members of the early growth response gene family, EGR3 and EGR4, were induced by low pH conditions. Amongst the immune response genes, an increased expression of decay accelerating factor (DAF) was observed, which has previously been reported to be upregulated in BO (28). We also noted the expression of interleukin-8 (IL-8) and prostaglandin synthase (PTGS2 or COX2), both of which are upregulated in the Barrett's adenocarcinoma sequence (2931). The transcription factors FOSB and JUN of the AP-1 transcription complex, known to have increased DNA binding potential under conditions of low pH, are upregulated in the differentiation/proliferation group (Figure 2) (32). Many of the identified differentially expressed genes have functional roles in stress response, proliferation, differentiation, inflammation and growth arrest. Other important transcription factors identified using this approach are the Barrett's associated genes activating transcription factor 3 (ATF3) and nuclear receptor 4A2 (NR4A2) (33).
|
Validation of gene-expression data
A cohort of genes were chosen for downstream validation, varying from growth factors (CTGF, CYR61, C8FW) to signalling molecules (MKP-1, DUSP5, MNT, EGR3, ATF3, FOSB). Changes observed in the DNA microarray studies of a subset of genes were confirmed by semi-quantitative RTPCR analysis (Figure 3A). These included MKP-1, which was rapidly upregulated after 30 min, and CTGF, which remained upregulated across the time periods. Additionally, ATF3 and EGR3 were validated by this method, showing rapid induction by 30 min increasing incrementally over time. Real-time RTPCR analysis permitted examination of a much broader range of gene-expression changes (Figure 3B and C). GAPDH was used to normalize all genes in the real-time RTPCR analysis including MAN2C, which was chosen from the array dataset as a non-changing gene to determine baseline fluctuations in the experimental approach. EGR3 and ATF3 have similar induction profiles with maximal induction at 2 h, which may be indicative of similar transcriptional control mechanisms (Figure 3B). The FOSB gene was observed to be strikingly altered, with an
200-fold increase in expression level observed at 120 min post-treatment as compared with that seen in the resting state (Figure 3B). Induction of EGR3 and ATF3 was also demonstrated on a protein level (Figure 3C). The importance of these genes in low pH signalling is highlighted by their responses in other oesophageal cell lines, such as OE-33 and the pre-neoplastic oesophageal cell line HET-1A (Figure 4). The oesophageal adenocarcinoma cell line OE-33 and the oesophageal squamous epithelial cell line HET-1A were treated with low pH 6.5 as before and gene expression was assessed using real-time RTPCR analysis (Figure 4). Overall, the data confirm that the gene inductions observed are not cell line specific and are maintained to different degrees in cell lines of different malignant stages.
|
|
Cluster analysis of DNA microarray data
Similar expression profiles for EGR3 and ATF3 were observed using the GeneSpring software and the subsequent transcript and protein validation (Figure 3), indicating that these genes are under similar transcriptional control. In order to assess this possibility using promoter inspection software, it was necessary to identify distinct groups of profiles for comparison rather than one single profile. These profiles were obtained by clustering the top 100 most significantly altered genes as described by GeneSpring. Figure 5A details the unbiased hierarchical tree cluster of this ordered list, with three groups of genes identified and marked for in silico promoter analysis. We selected the three clusters for further promoter inspection based on the presence of verified and/or Barrett's related genes within the clusters.
|
Common regulatory elements within the promoter regions of individual clusters
A flow chart for the comparative promoter analysis is shown in Figure 5, beginning with the hierarchical tree and selected clusters in Figure 5A. Figure 5B details the expression profile of each of the clusters chosen for promoter analysis. Cluster 1 contains genes displaying a rapid increase in expression in response to low pH, such as EGR3 and ATF3. Cluster 2 contains CYR61 and LDLR genes that reach maximal expression after 2 h, whereas Cluster 3 includes genes of the metallothionein family which have a latent maximal expression at 4 h. Promoter regions of the genes residing in the three clusters were examined for putative transcription factor binding sites which may co-ordinately regulate gene transcription. A fourth cluster of randomly selected non-changing genes was also chosen as an internal control for randomly occurring TFBSs. We identified the increased or significant occurrence of the binding sites within each cluster in comparison with each other and with the control set described. Figure 5C lists the genes residing within each of the cluster described and the most significant set of transcription factor binding elements found in their promoter regions is detailed in Figure 5D. Interestingly, Cluster 1 (EGR genes) may be controlled by the stress response factor Oct-1 and Cluster 2 by factors such as NF-KB and EGRF. Both NF-KB and AP-1 have been shown to have increased DNA binding potential under conditions of low pH (14,32,34). Furthermore, genes upregulated early in response to low pH, such as those of the EGR family (Cluster 1), may be involved in the co-ordinate upregulation of the genes contained in Cluster 2 through the EGRF binding site identified. The co-ordinate relationship between the genes can be verified in correlation searches on the ATF3 gene using Oncomine 2.0 (35). Using the array data from Welsh et al. (36) prostate cancer study (oesophageal correlation not available), a correlation search was performed for the Barrett's related gene ATF3 (Figure 6) (33). A number of the co-regulated genes identified in this study (EGR3, EGR1, FOSB, NR4A1) correlated highly with ATF3 expression, thus indicating that the clusters identified by this study may also occur in vivo.
|
| Discussion |
|---|
|
|
|---|
This study demonstrates that low pH induces a co-ordinate set of transcriptional responses in oesophageal cells, resulting in expression of genes involved in inflammation, stress response, proliferation and differentiation. This series of responses has been mapped to identify the impact of sequential transcriptional events. These findings identify genes induced by low pH which may be of potential importance in the pathogenesis of reflux oesophagitis and indirectly for the Barrett's cancer sequence.
The in silico transcription factor analysis applied in this study represents the first application of comparative binding site analysis to clustered data from a low pH gene induction study. Studies performed on cell lines with only one variable over time, i.e. low pH, are ideal candidates for TFBS analysis. This method identified TFBS sets found in the promoter regions which may co-ordinately regulate the genes within the given clusters. The incorporation of an expression cluster consisting of randomly selected non-changing genes is a method whereby randomly occurring non-functional TFBSs may be controlled and eliminated from the analysis. The discovery of increased presence of both NF-KB and EGRF binding elements in the promoter regions of Cluster 2 is of major interest for a number of reasons. The progress from GORD through to BO and adenocarcinoma passes through a gross inflammatory stage known as oesophagitis, and as with many other forms of inflammatory conditions can be characterized by the activation of the transcription factor NF-KB resulting in the expression of inflammation-related genes, such as IL-8 (13,37). This master regulatory factor is also known to be upregulated in BO and overexpressed in 60% of oesophageal adenocarcinomas, where it inversely correlates with response to neo-adjuvant chemotherapy (14). Chromosome 4 on which NF-KB resides has also been shown to be amplified during oesophageal cancer progression (3840). The above findings and previous reporting of NF-KB activation by low pH and bile acids lends weight to the importance of this transcription factor and its responsive genes in the progression of disease and as a potential therapeutic target. The promoters of Cluster 2 also contain the EGRF element which allows EGR14 to bind and activate transcription. This family of transcription factors are important inflammatory mediators and are often classed as oncogenes owing to their deregulation in many tumour types. EGR1, which is the most characterized member of the family can promote angiogenesis within tumour cores (41,42). The low pH responsive genes EGR1 and EGR3 identified in this study are members of the first expression cluster and this marks the beginning of what is termed a regulatory inference network (43,44). Transcriptional ordering or regulatory network inference is an area of bioinformatics under development as our understanding of gene networks and promoter function grows. We believe the methods employed in this study can add to the tools available to analyse the vast amounts of data generated from transcriptomic experiments and assist in identifying sequential biological events. Further validation and dissection of this regulatory network using such techniques as chromatin immunoprecipitation will improve the current understanding of the relationship between promoter TFBSs and the genes they regulate.
In order to assess the effects of low pH on gene expression, we exposed oesophageal cells to a sustained change in pH over the period of the experiment rather than recapitulating the real-life scenario of intermittent brief exposures over a lifetime. We have previously shown in oesophageal and gastric cells that pH 6.5 can activate NF-KB, a key mediator of inflammation and disease progression (14,15), and the transcription factor EGR1 identified in this study. There have been other reports of AP-1 DNA binding activity at the pH ranges used in this study (32). Therefore, it was of interest to determine the genes downstream of these transcription factors induced by low pH conditions. Sustained exposure to pH levels <6 resulted in loss of induction of NF-KB and it would not have been possible to identify downstream effects of this transcription factor at lower pH levels (15). In addition, preliminary experiments (data not shown) demonstrated that continuous culture in pH 3.5 led to a significant loss in cell viability as determined by PI staining and MTT assay and, as such, was not suitable for detailed molecular analysis of low pH signalling.
The identification of EGR1 as a low pH-responsive gene in oesophageal cells allowed for its use as a positive control gene in the subsequent set of transcriptomic experiments performed. The initial broad pH range used demonstrated EGR1 protein induction from pH 7 to 6; however, when these experiments were performed on a transcriptional level with a more focused pH range, it was discovered that EGR1 transcript induction occurred between pH 6.8 and 6.2. Therefore, our choice of pH 6.5 detailed previously was within the induction pH range for both EGR1 protein and mRNA. The analysis of the data generated from array experimentation was critically assessed using two methods, the MannWhitney protocol developed by Jin et al. (26), followed by complementary analysis using GeneSpring, thus revealing genes whose expression had been altered by low pH conditions.
Subsequent to the array experimentation, the low pH-responsive genes identified were allotted into groups based on similarity in their function as described by literature searches. The data analysis identified a small subset of genes whose expression was downregulated but only by small fold changes that were deemed not statistically significant, hence their omission from Figure 2. The full dataset is available for assessment in the GEO repository (GSE2144 [NCBI GEO] ). Interestingly, in this study the expression of two important phosphatases, DUSP1 and 5, are upregulated. These genes are detailed in the signalling group (Figure 2) and are capable of dephosphorylating members of the SAPK, p38, JNK and MAPK family thus inactivating the pathways. Their induction by low pH at an early time point is indicative of genotoxic stress, environmental stress and MAP kinase pathway activation as this pathway attenuates its own induction in a pulsatile fashion (45,46). Although DUSP1 is found in the cytoplasm where it may inhibit activation of the JNK pathway and thereby inhibit apoptosis, DUSP5 resides in the nucleus and attenuates ERK activation (47,48). The induction and validation of the phosphatases in this study in conjunction with other members of the dual specificity phosphatases found in other low pH signalling models (49) lend weight to the importance of this family of proteins in low pH signalling. The stimulation of IL-8, urokinase-type plasminogen receptor and prostaglandin synthase 2 (PTGS2 or COX2) by low pH allows the activation of the immune system which is often deregulated in many epithelial tumours and BO (13,37). The activation of these genes is a hallmark of NF-KB activity which in turn controls the activation of the inflammatory response to stressful injury (13,40).
DAF, a complement regulatory protein, is another gene in the immune response category (Figure 2) induced by low pH conditions in this study that has previously been used as a biomarker for specialized intestinal metaplasia type BO (28). Other documented transcriptional responses to low pH confirmed in this study to have potential differential expression in BO include CYR61, which is a recently discovered angiogenesis factor, and Amphiregulin, which is an EGF receptor ligand from the differentiation and proliferation category (Figure 2) (33). ATF3 and NR4A2 orphan receptor, from the signalling group are also associated with BO (33). The expression of ATF3, a member of the ATF/CREB subfamily of basic-region leucine zipper (bZIP) proteins, can be induced by a number of stimuli, such as TGF-ß, TNF-
, DNA damage and other stressful conditions (50). To date, this gene has been seen to perform many roles in the form of a transcriptional repressor, for example, the repression of the Id inhibitors of differentiation through such mediators as Smad3 (51). Induction of this gene by low pH conditions in all three oesophageal cell lines used is indicative of oxidative damage to DNA, which can occur over time in acidic environments (50,52). The induction of some genes in both normal and adenocarcinoma cell types indicates that they may not be specific to the metaplasia
adenocarcinoma sequence, but it does not preclude their involvement in oesophageal inflammatory pathogenesis or carcinogenesis.
In conclusion, we have demonstrated that low pH conditions can alter the expression of genes with diverse functions which may have potential association with BO. Novel bioinformatic techniques identified regulators of low pH gene induction, such as NF-KB and EGRF and the clusters of genes that they potentially regulate. Therapeutic intervention targeting these transcription factors and/or their regulatory mechanisms may be used to complement existing therapies to halt the development and progression of BO.
| Acknowledgments |
|---|
We thank Alison Murphy of the Conway Institute for technical assistance and Dr Siobhan Mitchell of the IMM for technical support. This grant was supported by Higher Education Authority (HEA) and the Program for Research in Third Level Institutions (PRTLI).
Conflict of Interest Statement: None declared.
| References |
|---|
|
|
|---|
- McManus,D.T., Olaru,A. and Meltzer,S.J. (2004) Biomarkers of esophageal adenocarcinoma and Barrett's esophagus. Cancer Res., 64, 15611569.
[Abstract/Free Full Text] - Jankowski,J.A., Wright,N.A., Meltzer,S.J., Triadafilopoulos,G., Geboes,K., Casson,A.G., Kerr,D. and Young,L.S. (1999) Molecular evolution of the metaplasia-dysplasia-adenocarcinoma sequence in the esophagus. Am. J. Pathol., 154, 965973.
[Abstract/Free Full Text] - Wijnhoven,B.P., Tilanus,H.W. and Dinjens,W.N. (2001) Molecular biology of Barrett's adenocarcinoma. Ann. Surg., 233, 322337.[CrossRef][ISI][Medline]
- Lagergren,J., Bergstrom,R., Lindgren,A. and Nyren,O. (1999) Symptomatic gastroesophageal reflux as a risk factor for esophageal adenocarcinoma. N. Engl. J. Med., 340, 825831.
[Abstract/Free Full Text] - Katz,P.O. (2004) Management of the patient with Barrett's esophagus: a continuing dilemma for the clinician. Rev. Gastroenterol Disord., 4, 4959.[Medline]
- Eda,A., Osawa,H., Satoh,K., Yanaka,I., Kihira,K., Ishino,Y., Mutoh,H. and Sugano,K. (2003) Aberrant expression of CDX2 in Barrett's epithelium and inflammatory esophageal mucosa. J. Gastroenterol., 38, 1422.[CrossRef][ISI][Medline]
- Tselepis,C., Perry,I., Dawson,C., Hardy,R., Darnton,S.J., McConkey,C., Stuart,R.C., Wright,N., Harrison,R. and Jankowski,J.A. (2002) Tumour necrosis factor-alpha in Barrett's oesophagus: a potential novel mechanism of action. Oncogene, 21, 60716081.[CrossRef][ISI][Medline]
- Riddell,R.H. (2005) The genesis of Barrett esophagus: has a histologic transition from gastroesophageal reflux disease-damaged epithelium to columnar metaplasia ever been seen in humans? Arch. Pathol. Lab. Med., 129, 164169.[ISI][Medline]
- Jankowski,J.A., Harrison,R.F., Perry,I., Balkwill,F. and Tselepis,C. (2000) Barrett's metaplasia. Lancet, 356, 20792085.[CrossRef][ISI][Medline]
- Fitzgerald,R.C. and Abdalla,S. (2003) Gastrin-induced hyperproliferation in Barrett's esophagus. Gastroenterology, 125, 1921; author reply 19211922.[ISI][Medline]
- Fitzgerald,R.C., Omary,M.B. and Triadafilopoulos,G. (1997) Acid modulation of HT29 cell growth and differentiation. An in vitro model for Barrett's esophagus. J. Cell Sci., 110 (Pt 5), 663671.[Abstract]
- Fitzgerald,R.C., Omary,M.B. and Triadafilopoulos,G. (1998) Altered sodium-hydrogen exchange activity is a mechanism for acid-induced hyperproliferation in Barrett's esophagus. Am. J. Physiol., 275, G47G55.[Medline]
- Van Der Woude,C.J., Kleibeuker,J.H., Jansen,P.L. and Moshage,H. (2004) Chronic inflammation, apoptosis and (pre-)malignant lesions in the gastro-intestinal tract. Apoptosis, 9, 123130.[CrossRef][ISI][Medline]
- Abdel-Latif,M.M., O'Riordan,J., Windle,H.J., Carton,E., Ravi,N., Kelleher,D. and Reynolds,J.V. (2004) NF-kappaB activation in esophageal adenocarcinoma: relationship to Barrett's metaplasia, survival, and response to neoadjuvant chemoradiotherapy. Ann. Surg., 239, 491500.[CrossRef][ISI][Medline]
- O'Toole,D., Abdel-Latif,M.M., Long,A., Windle,H.J., Murphy,A.M., Bowie,A., O'Neill,L.A.J., Weir,D.G. and Kelleher,D. (2005) Low pH and Helicobacter pylori increase nuclear factor kappa b binding in gastric epithelial cells: a common pathway for epithelial cell injury? J. Cell Biochem., [Epub ahead of print].
- Torigoe,T., Izumi,H., Yoshida,Y., Ishiguchi,H., Okamoto,T., Itoh,H. and Kohno,K. (2003) Low pH enhances Sp1 DNA binding activity and interaction with TBP. Nucleic Acids Res., 31, 45234530.
[Abstract/Free Full Text] - Skerka,C., Decker,E.L. and Zipfel,P.F. (1995) A regulatory element in the human interleukin 2 gene promoter is a binding site for the zinc finger proteins Sp1 and EGR-1. J. Biol. Chem., 270, 2250022506.
[Abstract/Free Full Text] - Berendji,D., Kolb-Bachofen,V., Zipfel,P.F., Skerka,C., Carlberg,C. and Kroncke,K.D. (1999) Zinc finger transcription factors as molecular targets for nitric oxide-mediated immunosuppression: inhibition of IL-2 gene expression in murine lymphocytes. Mol Med., 5, 721730.[ISI][Medline]
- Abdel-Latif,M.M., Windle,H.J., Fitzgerald,K.A., Ang,Y.S., Eidhin,D.N., Li-Weber,M., Sabra,K. and Kelleher,D. (2004) Helicobacter pylori activates the early growth response 1 protein in gastric epithelial cells. Infect Immun., 72, 35493560.
[Abstract/Free Full Text] - McMahon,S.B. and Monroe,J.G. (1996) The role of early growth response gene 1 (egr-1) in regulation of the immune response. J. Leukoc. Biol., 60, 159166.[Abstract]
- Sukhatme,V.P., Cao,X.M., Chang,L.C. et al. (1988) A zinc finger-encoding gene coregulated with c-fos during growth and differentiation, and after cellular depolarization. Cell, 53, 3743.[CrossRef][ISI][Medline]
- Sukhatme,V.P. (1990) Early transcriptional events in cell growth: the Egr family. J. Am. Soc. Nephrol., 1, 859866.[Abstract]
- Altorki,N., Schwartz,G.K., Blundell,M., Davis,B.M., Kelsen,D.P. and Albino,A.P. (1993) Characterization of cell lines established from human gastric-esophageal adenocarcinomas. Biologic phenotype and invasion potential. Cancer, 72, 649657.[CrossRef][ISI][Medline]
- Heid,C.A., Stevens,J., Livak,K.J. and Williams,P.M. (1996) Real time quantitative PCR. Genome Res., 6, 986994.
[Abstract/Free Full Text] - Gibson,U.E., Heid,C.A. and Williams,P.M. (1996) A novel method for real time quantitative RTPCR. Genome Res., 6, 9951001.
[Abstract/Free Full Text] - Jin,H., Yang,R., Awad,T.A. et al. (2001) Effects of early angiotensin-converting enzyme inhibition on cardiac gene expression after acute myocardial infarction. Circulation, 103, 736742.
[Abstract/Free Full Text] - Quandt,K., Frech,K., Karas,H., Wingender,E. and Werner,T. (1995) MatInd and MatInspector: new fast and versatile tools for detection of consensus matches in nucleotide sequence data. Nucleic Acids Res., 23, 48784884.
[Abstract/Free Full Text] - Hiraoka,S., Mizuno,M., Nasu,J., Okazaki,H., Makidono,C., Okada,H., Terada,R., Yamamoto,K., Fujita,T. and Shiratori,Y. (2004) Enhanced expression of decay-accelerating factor, a complement-regulatory protein, in the specialized intestinal metaplasia of Barrett's esophagus. J. Lab. Clin. Med., 143, 201206.[CrossRef][ISI][Medline]
- Kuramochi,H., Vallbohmer,D., Uchida,K. et al. (2004) Quantitative, tissue-specific analysis of cyclooxygenase gene expression in the pathogenesis of Barrett's adenocarcinoma. J. Gastrointest. Surg., 8, 10071016; discussion 10161017.
- Kawabe,A., Shimada,Y., Soma,T., Maeda,M., Itami,A., Kaganoi,J., Kiyono,T. and Imamura,M. (2004) Production of prostaglandinE2 via bile acid is enhanced by trypsin and acid in normal human esophageal epithelial cells. Life Sci., 75, 2134.[CrossRef][ISI][Medline]
- Shirvani,V.N., Ouatu-Lascar,R., Kaur,B.S., Omary,M.B. and Triadafilopoulos,G. (2000) Cyclooxygenase 2 expression in Barrett's esophagus and adenocarcinoma: Ex vivo induction by bile salts and acid exposure. Gastroenterology, 118, 487496.[CrossRef][ISI][Medline]
- Shimokawa,N., Kumaki,I., Qiu,C.H., Ohmiya,Y., Takayama,K. and Koibuchi,N. (2005) Extracellular acidification enhances DNA binding activity of MafG-FosB heterodimer. J. Cell Physiol., 205, 7785.
- Barrett,M.T., Yeung,K.Y., Ruzzo,W.L., Hsu,L., Blount,P.L., Sullivan,R., Zarbl,H., Delrow,J., Rabinovitch,P.S. and Reid,B.J. (2002) Transcriptional analyses of Barrett's metaplasia and normal upper GI mucosae. Neoplasia, 4, 121128.[CrossRef][ISI][Medline]
- Xu,L., Fukumura,D. and Jain,R.K. (2002) Acidic extracellular pH induces vascular endothelial growth factor (VEGF) in human glioblastoma cells via ERK1/2 MAPK signaling pathway: mechanism of low pH-induced VEGF. J. Biol. Chem., 277, 1136811374.
[Abstract/Free Full Text] - Rhodes,D.R., Yu,J., Shanker,K., Deshpande,N., Varambally,R., Ghosh,D., Barrette,T., Pandey,A. and Chinnaiyan,A.M. (2004) ONCOMINE: a cancer microarray database and integrated data-mining platform. Neoplasia, 6, 16.[ISI][Medline]
- Welsh,J.B., Sapinoso,L.M., Su,A.I., Kern,S.G., Wang-Rodriguez,J., Moskaluk,C.A., Frierson,H.F.Jr and Hampton,G.M. (2001) Analysis of gene expression identifies candidate markers and pharmacological targets in prostate cancer. Cancer Res., 61, 59745978.
[Abstract/Free Full Text] - Konturek,P.C., Nikiforuk,A., Kania,J., Raithel,M., Hahn,E.G. and Muhldorfer,S. (2004) Activation of NFkappaB represents the central event in the neoplastic progression associated with Barrett's esophagus: a possible link to the inflammation and overexpression of COX-2, PPARgamma and growth factors. Dig. Dis. Sci., 49, 10751083.[CrossRef][ISI][Medline]
- Croft,J., Parry,E.M., Jenkins,G.J., Doak,S.H., Baxter,J.N., Griffiths,A.P., Brown,T.H. and Parry,J.M. (2002) Analysis of the premalignant stages of Barrett's oesophagus through to adenocarcinoma by comparative genomic hybridization. Eur. J. Gastroenterol Hepatol., 14, 11791186.[CrossRef][ISI][Medline]
- Doak,S.H., Jenkins,G.J., Parry,E.M., D'Souza,F.R., Griffiths,A.P., Toffazal,N., Shah,V., Baxter,J.N. and Parry,J.M. (2003) Chromosome 4 hyperploidy represents an early genetic aberration in premalignant Barrett's oesophagus. Gut, 52, 623628.
[Abstract/Free Full Text] - Jenkins,G.J., Harries,K., Doak,S.H., Wilmes,A., Griffiths,A.P., Baxter,J.N. and Parry,J.M. (2004) The bile acid deoxycholic acid (DCA) at neutral pH activates NF-kappaB and induces IL-8 expression in oesophageal cells in vitro. Carcinogenesis, 25, 317323.
[Abstract/Free Full Text] - Adamson,E.D. and Mercola,D. (2002) Egr1 transcription factor: multiple roles in prostate tumor cell growth and survival. Tumour Biol., 23, 93102.[CrossRef][Medline]
- Krones-Herzig,A., Adamson,E. and Mercola,D. (2003) Early growth response 1 protein, an upstream gatekeeper of the p53 tumor suppressor, controls replicative senescence. Proc. Natl Acad. Sci. USA, 100, 32333238.
[Abstract/Free Full Text] - Haverty,P.M., Hansen,U. and Weng,Z. (2004) Computational inference of transcriptional regulatory networks from expression profiling and transcription factor binding site identification. Nucleic Acids Res., 32, 179188.
[Abstract/Free Full Text] - Haverty,P.M., Frith,M.C. and Weng,Z. (2004) CARRIE web service: automated transcriptional regulatory network inference and interactive analysis. Nucleic Acids Res., 32, W213W216.
[Abstract/Free Full Text] - Ryser,S., Massiha,A., Piuz,I. and Schlegel,W. (2004) Stimulated initiation of mitogen-activated protein kinase phosphatase-1 (MKP-1) gene transcription involves the synergistic action of multiple cis-acting elements in the proximal promoter. Biochem. J., 378, 473484.[CrossRef][ISI][Medline]
- Barry,O.P., Mullan,B., Sheehan,D., Kazanietz,M.G., Shanahan,F., Collins,J.K. and O'Sullivan,G.C. (2001) Constitutive ERK1/2 activation in esophagogastric rib bone marrow micrometastatic cells is MEK-independent. J. Biol. Chem., 276, 1553715546.
[Abstract/Free Full Text] - Mizuno,R., Oya,M., Shiomi,T., Marumo,K., Okada,Y. and Murai,M. (2004) Inhibition of MKP-1 expression potentiates JNK related apoptosis in renal cancer cells. J. Urol., 172, 723727.[CrossRef][ISI][Medline]
- Small,G.W., Shi,Y.Y., Edmund,N.A., Somasundaram,S., Moore,D.T. and Orlowski,R.Z. (2004) Evidence that mitogen-activated protein kinase phosphatase-1 induction by proteasome inhibitors plays an anti-apoptotic role. Mol. Pharmacol., 66,14781490.
[Abstract/Free Full Text] - Morgan,C., Alazawi,W., Sirieix,P., Freeman,T., Coleman,N. and Fitzgerald,R. (2004) In vitro acid exposure has a differential effect on apoptotic and proliferative pathways in a Barrett's adenocarcinoma cell line. Am. J. Gastroenterol., 99, 218224.[CrossRef][ISI][Medline]
- Fan,F., Jin,S., Amundson,S.A. et al. (2002) ATF3 induction following DNA damage is regulated by distinct signaling pathways and over-expression of ATF3 protein suppresses cells growth. Oncogene, 21, 74887496.[CrossRef][ISI][Medline]
- Kang,Y., Chen,C.R. and Massague,J. (2003) A self-enabling TGFbeta response coupled to stress signaling: Smad engages stress response factor ATF3 for Id1 repression in epithelial cells. Mol. Cell., 11, 915926.[CrossRef][ISI][Medline]
- Jolly,A.J., Wild,C.P. and Hardie,L.J. (2004) Acid and bile salts induce DNA damage in human oesophageal cell lines. Mutagenesis, 19, 319324.
[Abstract/Free Full Text]
![]()
CiteULike
Connotea
Del.icio.us What's this?
This article has been cited by other articles:
![]() |
K. Dvorak, C. M Payne, M. Chavarria, L. Ramsey, B. Dvorakova, H. Bernstein, H. Holubec, R. E Sampliner, N. Guy, A. Condon, et al. Bile acids in combination with low pH induce oxidative stress and oxidative DNA damage: relevance to the pathogenesis of Barrett's oesophagus Gut, June 1, 2007; 56(6): 763 - 771. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||






