- Open Access
Differential expression of apoptotic genes PDIA3 and MAP3K5 distinguishes between low- and high-risk prostate cancer
Molecular Cancer volume 8, Article number: 130 (2009)
Despite recent progress in the identification of genetic and molecular alterations in prostate cancer, markers associated with tumor progression are scarce. Therefore precise diagnosis of patients and prognosis of the disease remain difficult. This study investigated novel molecular markers discriminating between low and highly aggressive types of prostate cancer.
Using 52 microdissected cell populations of low- and high-risk prostate tumors, we identified via global cDNA microarrays analysis almost 1200 genes being differentially expressed among these groups. These genes were analyzed by statistical, pathway and gene enrichment methods. Twenty selected candidate genes were verified by quantitative real time PCR and immunohistochemistry. In concordance with the mRNA levels, two genes MAP3K5 and PDIA3 exposed differential protein expression. Functional characterization of PDIA3 revealed a pro-apoptotic role of this gene in PC3 prostate cancer cells.
Our analyses provide deeper insights into the molecular changes occurring during prostate cancer progression. The genes MAP3K5 and PDIA3 are associated with malignant stages of prostate cancer and therefore provide novel potential biomarkers.
Prostate cancer is the most frequent cancer diagnosed in men (20.3% of the total), followed by lung (17.2%) and colorectal cancer (12.8%) . Measuring prostate specific antigen (PSA) has been a matter of routine to detect prostate cancer, but is insufficient to distinguish between different tumor grades. The Gleason Grading System is commonly used for histology-based grading of prostate cancer tissue . Since prostate tumors are often multifocal, the Gleason Score (GS) is the sum of the two most prevalent tumor patterns, which are graded 1 (CA1) as the most differentiated and 5 (CA5) as the least differentiated pattern of cancerous glands. Other methods for sub-classification have been described in recent reports . These indicate that translocations fusing the strong androgen-responsive gene TMPRSS2 with ERG or other oncogenic ETS factors may facilitate prostate cancer development. It has been proposed that the presence or absence of this genetic rearrangement may be used, much like the Gleason grading system, as a diagnostic tool to extract prognostically relevant sub-classifications of this cancer .
The discrimination between different tumor grades is important with respect to treatment decisions: Currently, many men who are diagnosed with GS 6 prostate cancer are often "over"-treated and risk suffering from urinary and sexual dysfunction . Therefore, it is important to develop a sensitive and specific diagnostic tool to distinguish between different tumor grades. To address this problem, many groups have recently started to profile gene expression levels in prostate tumor tissues to identify deregulated genes during disease progression. However, although many of these have addressed the question of molecular differences between normal, tumor, benign prostatic hyperplasia (BPH), and the putative precursor lesion prostatic intraepithelial neoplasia (PIN), little is still known about molecular changes between low- and high-risk tumors [6–9].
In the present study, we performed microarray-based gene expression profile analysis of 65 microdissected tissues comprising 25 samples of GS 6, 27 of GS 8-10 and 13 non cancerous samples. We sought to identify biological markers of distinct functional groups for the discrimination between low- and high-risk tumors. Overall, we found 20 genes with a significant alteration in expression between high-risk compared to low-risk tumors. Two of these genes exhibited Gleason grade associated protein expression in tumor tissues, which could serve as a valuable diagnostic tool in the future.
mRNA expression analysis revealed large expression differences between GS 6 and GS 8-10 tumors
To selectively isolate pure populations of prostate epithelial cancer cells with different Gleason Scores, we first applied laser-capture microdissection. We monitored the gene expression levels by hybridization of twice-amplified RNA to cDNA microarrays representing ~37500 mapped genes. In total, we hybridized 65 RNA samples derived from 13 benign and 52 prostate cancer tissue comprising 25 samples with Gleason Score (GS) 6 and 27 samples with GS 8-10 (Table 1). After quality assessment of microarray hybridizations, we subjected gene expression profiles to SAM . Numbers of deregulated genes identified by SAM analyses are summarized in Table 2, and complete gene lists are provided (see additional file 1 and additional file 2).
For the identification of grade-discriminating genes, we compared the expression levels of GS 6 with GS 8-10 tumors. SAM analysis revealed 1141 up-regulated and 54 down-regulated non-redundant genes in advanced tumors (FDR 5%; see additional file 1). For validation, we compared our data with an independent study from True and coworkers, who reported 86 genes as deregulated during tumor progression from low to high GS . Of these, we identified 24 genes (28%) which all displayed the same tendency as in the original report (see additional file 3). Another comparison to the study of Lapointe and coworkers , who described 41 genes to be associated to a higher Gleason score revealed an overlap of six genes (BGN, COL1A2, COL3A1, PLA2G2A, SPARC, VCAN).
To identify biological processes associated with tumor progression, we performed gene ontology analysis with genes differentially regulated between low- and high-risk tumors. In order to extract highly significant canonical pathways, each gene symbol was mapped to its corresponding gene object in the IPA Knowledge Base, and networks were generated. Significant canonical pathways were related to actin-mediated processes, e.g. regulation of actin-based motility mediated by Rho-family GTPases (11 of 92 annotated genes, p = 5.23E-03) and actin cytoskeleton signaling (19 of 221 annotated genes, p = 1.4E-02). In addition, various metabolic processes, including oxidative phosphorylation (22 of 158 annotated genes; p = 8.95E-06) and protein ubiquitination (19 of 205 annotated genes; p = 5.6E-03) were deregulated in high-risk tumors. In order to extract as many genes as possible involved in apoptotic processes we used two further GO analysis tools (FatiGO  and GOstat ) and identified a set of 46 genes associated with apoptosis (GOstat p = 0.00033; see additional file 4).
Additionally, we compared gene signatures between normal and tumor tissue, which lead to the identification of > 2500 deregulated genes (FDR 5%) of which 2390 genes were down and 243 up-regulated (see additional file 2). We performed separate SAM analyses between normal and GS 6 or GS 8-10 and revealed 2016 and 2001 down-regulated as well as 463 and 454 up-regulated non redundant genes (Table 2). Of these, 1197 genes were deregulated with the same tendency in both tumor groups. Interestingly, three genes (VCAN, CLK1 and TMEM16G) revealed opposite expression levels in these comparisons. VCAN and CLK1 were found to be significantly down-regulated in GS 6, but up-regulated in high-risk tumors in comparison to normal tissue. In agreement with these results, over-expression of VCAN and CLK1 between primary prostate cancer and metastatic cancer has recently been described . In contrast TMEM16G, for which the microarray findings were supported by qRT-PCR data, showed the opposite trend. TMEM16G was found up-regulated between normal and GS 6 tissue, but down-regulated between normal and GS 8 as well as between GS 6 and GS 8 tissues. In a recent study, TMEM16G was described as a prostate-specific plasma membrane protein promoting cell-cell contact in the prostate cancer cell line LNCaP .
Validation of selected genes by qRT-PCR
Based on the microarray gene expression differences between GS 6 and GS 8-10 tumors (see additional file 1), 68 genes were validated in the same cohort via quantitative real time PCR (see additional file 5). B2M (β2-microglobulin) was used as a housekeeping gene due to its even expression in all analyzed patient groups (Figure 1A). We focused on a selection of genes that are linked to cancer-relevant gene ontology categories like apoptosis, cell morphology, metabolism and ubiquitylation. For example, 24/68 selected genes are functionally associated to apoptosis (see additional file 4).
In the qRT-PCR analysis of the 68 genes, 23 were significantly deregulated between low- and high-risk tumors (Wilcoxon p < 0.1; see additional file 5). Twenty of these (87%; Table 3) were in concordance with the microarray results. Thirty-three genes were determined as non-significant and 12 genes exhibited Ct values below detection level. Deregulation of the apoptotic process plays a major role in tumorigenesis and influences therapeutic outcome . In total, 11 of 20 significantly verified genes are involved in apoptotic processes (PDIA3, MAP3K5, ANXA5, VDAC1, NGFRAP1, TEGT, NPM1, VCP, TRAF4, HMGB1, and ROCK1). The two genes MAP3K5 and PDIA3 were selected for in-depth analysis because of their association to apoptosis and previous findings in cancer studies. QRT-PCR expression patterns of PDIA3 and MAP3K5 are given in Figure 1B and 1C.
Immunohistochemistry demonstrates Gleason-grade associated protein expression of MAP3K5 and PDIA3
To confirm our data at the protein level, we performed immunohistochemical analysis of the proteins MAP3K5 and PDIA3 representing the largest functional group (apoptosis) of validated genes. Specificity of antibodies was controlled using western blotting (data not shown). Protein expression levels in tumor tissue samples were scored according to a 4 point scoring system. Lowest expression levels of MAP3K5 and PDIA3 proteins were seen in benign epithelial cells (Figure 2). In agreement with the transcript analyses, MAP3K5 exhibited a significant Gleason grade-associated protein expression (p < 0.01, Wilcoxon signed rank test). Highest expression levels were observed in Gleason pattern 4 regions (mean 1.8) whereas Gleason pattern 3 (mean 1.2) and Gleason pattern 5 tumor regions (mean 1.19) displayed lower immunoreactivity. Of note, we also observed significant protein overexpression in prostate intraepithelial neoplasia (PIN) and in regions of inflammation (data not shown), which is in agreement with the described involvement of MAP3K5 with inflammation processes. Immunostaining was observed in the cytoplasm.
PDIA3 also showed a grade-associated protein expression (p < 0.065, Wilcoxon signed rank test; Figure 2). Expression was significantly increased in tumor cells compared to benign epithelium. Among different Gleason patterns in the tumors, Gleason pattern 4 (GP4) showed the highest expression level. Interestingly, expression levels of PDIA3 in Gleason pattern 3 (CA3) seem to depend on the accompanying Gleason pattern in the tumor. PDIA3 expression was higher in presence of higher Gleason patterns (CA3 in GS 8 tumors, e.g. CA3 associated with CA5; mean score 1.97) than in tumors with uniform Gleason pattern 3 (CA3 in GS6 tumors; mean score 1.74). Immunostainings showed cytoplasmic and perinuclear localization of PDIA3. In advanced tumors, PDIA3 and MAP3K5 displayed a heterogeneous staining pattern and the variances of intensity distributions were higher. Representative pictures of PDIA3 and MAP3K5 in different tissue regions are shown in Figure 3, where AMACR serves as a positive control for tumor cells .
Decreased apoptotic activity upon knockdown of PDIA3 in prostate cancer cell lines
MAP3K5 and PDIA3 were found to be proteins associated with apoptotic processes via pathway analysis. Unlike MAP3K5, whose pro-apoptotic and inflammatory role in context of tumorigenesis is well established [16, 17], the function for PDIA3 in apoptosis has been largely unexplored. To investigate an apoptosis-related function of PDIA3 we performed siRNA-based knockdown in the human prostate cancer cell lines PC3 and LNCaP. 48 hours after siRNA treatment (20 nM or 40 nM) the knockdown efficiency was determined by qRT-PCR (Figure 4A). Induction of apoptosis was mediated by three different stimuli. Staurosporine (STS), Fenretinide (FenR) and Tapsigargin (TG) are known to activate apoptosis via distinct mechanisms [18–20]. Each stimulus activated the apoptotic pathway reflected by activation of caspase 3 and/or caspase 7 (Figure 4B). PDIA3 siRNA treatment revealed a significant decrease of caspase activation in PC3 cells with all stimuli in comparison to control siRNA treated cells. In LNCaP cells similar results were obtained, but were only significant after STS induction (data not shown). These results indicate a novel, pro-apoptotic role for PDIA3 in prostate cancer cells.
Recent studies showed that it is important to include different tumor stages of prostate cancer in gene expression analyses to be able to find new diagnostic and prognostic markers [6, 8, 9]. Here, we generated gene expression profiles of tissues from low-risk (GS 6) and high-risk prostate tumors (GS 8-10) tissues. In contrast to most other published studies, all tissue samples were carefully microdissected before RNA isolation. The comparison of these profiles revealed that both tumor subgroups differ by a large number of genes, most of which are up-regulated in high GS tumors. A comparison with another published data set  suggested that these results reflect a general trend in transcriptional activation in advanced prostate tumors. However, it cannot be fully ruled out that systematic changes introduced by different extents of stromal cells  in the microdissected material as well as two rounds of RNA amplification contribute to this bias.
Among the detected genes many are involved in pathways, which are known to be altered in tumor progression such as apoptosis, morphologic changes, metabolism, and ubiquitin-mediated protein degradation. Twenty representatives of these processes were verified via qRT-PCR, and a set of genes discriminating between less and more aggressive tumor forms was identified. One of the hallmarks of aggressive cancer is the imbalance between cell survival and apoptosis. Our gene ontology analysis revealed that a pronounced number of apoptosis-related genes exhibited expression changes between low- and high-risk prostate tumors. Thus, for validation, we focused our analysis on up-regulated anti-apoptotic genes and key players of apoptotic signaling and verified the expression level of 11 selected genes. These data were supported by the analysis of protein expression levels for MAP3K5 and PDIA3 by IHC. MAP3K5 and PDIA3 were also identified in other prostate cancer profiling studies [6, 7, 9, 21]. MAP3K5 (also known as apoptosis signal-regulating kinase 1; ASK1) has been widely accepted as one of the key components regulating reactive oxygen species (ROS) - induced JNK and p38 activation leading to differentially regulated apoptosis . ROS - dependent activation of MAP3K5 also plays a critical role in innate immune responses through production of proinflammatory cytokines . There is considerable evidence suggesting that oxidative stress contributes to the pathogenesis of prostate cancer [24, 25]. Given that mitochondria are a major source of reactive oxygen species (ROS), altered mitochondrial bioenergetics might induce MAP3K5 over-expression and contribute to the malignant progression of prostate tumors. In concordance with this hypothesis, we also found a significant number of deregulated genes involved in oxidative phosphorylation and mitochondrial dysfunction.
Like MAP3K5, PDIA3 (protein disulfide isomerase A3) is a member of the endoplasmatic reticulum stress signaling pathway also known as unfolded protein response (UPR), and its expression level increases in response to cellular stress due to its function as a chaperone . Recently published data connected PDIA3 to the apoptotic process and demonstrated an anti-apoptotic effect of PDIA3 in the melanoma cell line A375 after induction of ER stress . In contrast, our study suggested a decrease of caspase activity due to down regulation of PDIA3 in prostate cancer cell lines. This result suggests that the observed increase of PDIA3 in this study is most likely due to elevated cellular stress. But besides the role as a chaperone, PDIA3 might function as a pro-apoptotic protein in the prostate. Taking our IHC data of PDIA3 into account, PDIA3 protein concentration decreases significantly in CA5 compared to CA4 tissues and expression data comparing localized with metastatic prostate cancer showed a down-regulation of PDIA3 . These findings support the idea that down regulation of PDIA3 might play a role in late onset of prostate cancer progression. A lack of PDIA3 expression also correlates with increased tumor invasion and advanced stage of gastric cancer and has therefore been proposed to be a negative prognostic marker .
In addition to its role in the ER stress pathway, PDIA3 has recently gained attention due to its function as a component of the peptide-loading complex of the major histocompatibility complex (MHC) class I pathway [29, 30]. In PDIA3 deficient mice this complex is impaired and negatively influences presentation of antigenic peptides. This may help tumors to escape from immune surveillance by cytotoxic T cells .
The results of our IHC analysis point to a potential use of PDIA3 as a diagnostic marker: PDIA3 expression of Gleason pattern 3 tumors is higher in the presence of a Gleason pattern 5 tumor than in presence of another Gleason pattern 3 tumor. Additionally, PDIA3 and MAP3K5 have been found to be significantly (FDR 5%) up-regulated in tumors harboring a TMPRSS2 - fusion protein (data not shown). This underlines the opportunity to use PDIA3 and MAP3K5 as discriminating biomarkers in respect to histological grading system and gene arrangement classification.
In summary, this study validated a set of 20 genes, which discriminated between low and high Gleason grade prostate tumors. These genes comprise important functional processes well known to be involved in tumor progression such as apoptosis, morphological changes, metabolism and others. In addition, we show a grade-associated protein expression of MAP3K5 and PDIA3. In particular, high PDIA3 protein levels in Gleason pattern 3 cancers may indicate the presence of more aggressive tumor foci in the same tissue and could be of diagnostic value, possibly as part of a larger molecular signature.
Materials and methods
Frozen and paraffin-embedded prostate tissue samples were obtained from previously untreated patients who had undergone radical prostatectomy after tumor diagnosis in a PSA based screening program performed in Tyrol by the Department of Urology, Medical University of Innsbruck . The study was approved by the ethics committee at the Medical University of Innsbruck. Immediately after surgery, the prostate specimens were cooled in ice/water and brought to the pathologist who performed a rapid section and isolated a prostate slice that was embedded in Tissue-Tek OCT Compound (Sakura, Tokyo, Japan), snap frozen in liquid nitrogen and stored at -80°C until use. The rest of the prostate was fixed and paraffin-embedded according to standard procedures.
For isolation of total RNA, frozen sections were stained with hematoxylin and eosin for pathological analysis and exact localization of the tumors. Parallel unstained slides were used for microdissection. These were pre-treated for 1 min in each of the following pre-cooled solutions: 75% ethanol, RNase-free water, 100% ethanol (twice) and xylene (twice), and air dried. Laser-capture microdissection was performed on a Pix Cell II microdissection microscope (Arcturus, Sunnyvale, CA, USA) using 2,000-5,000 laser impulses corresponding to approximately 15 000 - 30 000 cells for each sample. Tumor samples were isolated from a cohort of Gleason score 6 tumors (Gleason pattern 3) and a group of Gleason score 8 - 10 tumors (Gleason patterns 4 and 5). Benign epithelial cell samples were microdissected apart from tumor foci from histopathologically normal regions of the same specimens. After microdissection, total RNA was isolated using the PicoPure isolation kit (MDC, Sunnyvale, USA) according to the protocol of the supplier. Quality control was done employing the Agilent Bioanalyzer 2100 system (Agilent Technologies, Waldbronn, Germany).
20 ng of RNA isolated from laser-capture microdissected epithelial prostate cells of tissues from patients who had undergone radical prostatectomy were subjected to a two-round amplification using the MessageAmpTM II aRNA Amplification Kit (Applied Biosystems/Ambion, Austin, USA). The quality of amplified RNA (aRNA) fragments was assessed by microcapillary electrophoresis using the Agilent Bioanalyzer 2100 system. Two micrograms of aRNA were subjected to microarray hybridization as described in . Briefly, aRNA was reverse transcribed using SuperScript II reverse transcriptase (Invitrogen, San Diego, CA, USA) and labeled with Cy5-dUTP. Each sample was compared to a common reference (Universal Human Reference RNA; Stratagene, La Jolla, CA, USA) labeled with Cy3-dUTP. Hybridizations were done using the platform Human Unigene3.1 cDNA Array 37.5K v1.0 (NCBI, GEO, GPL3050) representing estimated 22,000 transcripts. Data were analyzed using the GenePix 4.0 software (Axon Instruments, Foster City, USA). Low quality measurements were excluded from further analysis. Raw expression values were pre-processed using ArrayMagic  and thereby normalized using the VSN method .
Data analysis and data mining
Significance Analysis of Microarrays (SAM) was applied to identify genes differentially regulated between normal tissue, tumor tissue GS 6 and GS ≥ 8 . A two class unpaired SAM test with 1000 permutations was used. The False Discovery Rate (FDR) was set below 5%. Results from the SAM analysis were imported into FatiGO  and Ingenuity Pathways Analysis (IPA) software (Ingenuity Systems, Redwood City, CA, USA) to identify gene ontologies that were significantly over- or under-represented. MatchMiner software  was used to match gene entries between different microarray studies.
Quantitative real-time RT-PCR validation
Verification of expression of selected genes was performed via quantitative real-time RT-PCR using the ABI Prism 7900 HT Sequence Detection System (Applied Biosystems, Foster City, CA, USA) and the Universal Probe Library System (Roche, Basel, Switzerland). Ct values were extracted by using the SDS-software (Applied Biosystems). The expression level of the housekeeping gene β-2-microglobulin was used for normalization, calculated with the 2-ΔΔCt method . Gene expression differences between GS 6 and GS 8-10 tumors were analyzed using t-test. A list of examined genes including mean Ct values of each analyzed group and corresponding primers is given in additional file 5.
For validation of expression at the protein level, we used immunohistochemistry (IHC) on corresponding paraffin-embedded tissue specimens from the same patient cohort. Immunohistochemistry was performed with 5 μm paraffin tissue sections employing the Ventana Discovery - XT staining automat (Roche). Standard CC1 pre-treatment and antigen retrieval was followed by incubation with antibody solution for 1 hr, choice of amplification kit, universal antibody solution for 60 min, staining with DAP map kit and counter stain for 4 min with haematoxylin II bluing reagent (all from Roche). For expression analysis of MAP3K5 (ASK1) protein the mouse monoclonal antibody EP553Y (Abcam Limited, Cambridge, MA, USA) was used at a dilution of 1:40 for the PDIA3 protein the mouse monoclonal antibody MaPERp571 (Abcam) at a dilution of 1:10. Specificity of staining was controlled by including a control antibody (DAKO Cytomation, Glostrup, Denmark). Immunoreactivity was then scored by a uropathologist and stratified according to the histology and the Gleason pattern of the specimens using a 4 point scaling system: 0, no staining, 1, weak staining, 2, intermediate staining, 3, strong staining.
RNA interference and apoptosis assay
The human prostate cancer cell line PC3 and LNCaP were purchased from ATCC (Manassas, VA, USA) and cultured in RPMI 1640 or HAMs F12 medium, respectively. The medium was supplemented with 50 units/ml penicillin, 50 μg/ml streptomycin sulphate, 1% nonessential amino acids, and 10% FBS (all from GIBCO/BRL, Gaithersburg, MD, USA).
PDIA3 siRNA (target sequences, see additional file 6) was purchased from Dharmacon (Lafayette, CO, USA) and evaluated against a scrambled siRNA control (Qiagen, Hilden, Germany). The siRNA knockdown experiments were performed by plating 0.8 × 104 cells PC3 cells in a 96-well plate (NUNC, Roskilde, Denmark) overnight. For transfection, siRNA and Lipofectamine 2000 (Invitrogen, Carlsbad, CA, USA) were diluted separately and incubated for 5 min at room temperature. The two solutions were mixed and incubated for 20 min at room temperature. siRNA-Lipofectamine 2000 mixture was then added to the cells, and the plate was mixed by gentle rocking. Transfected cells were incubated at 37°C and 5% CO2 for 48 h. Knockdown efficiency was verified by qRT PCR.
Induction of apoptosis was performed by adding the indicated amounts of Staurosporine (Roche, Mannheim, Germany), Fenretinide (Sigma, Munich, Germany) or Tapsigargin (Sigma) for 6 and 24 hours, respectively. Control cells were left untreated. Activation of apoptosis was determined by measuring caspase 3 and 7 activities using Caspase-Glo 3/7® Assay (Promega, Madison, WI, USA) following the manufacturer's protocol.
The microarray data sets reported in this paper have been deposited MIAME compliant to NCBI Gene Expression Omnibus (GEO) database (accession no. GSE15484).
Ferlay J, Autier P, Boniol M, Heanue M, Colombet M, Boyle P: Estimates of the cancer incidence and mortality in Europe in 2006. Ann Oncol. 2007, 18: 581-592. 10.1093/annonc/mdl498
Gleason DF, Mellinger GT: Prediction of prognosis for prostatic adenocarcinoma by combined histological grading and clinical staging. J Urol. 1974, 111: 58-64.
Tomlins SA, Laxman B, Varambally S, Cao X, Yu J, Helgeson BE, Cao Q, Prensner JR, Rubin MA, Shah RB: Role of the TMPRSS2-ERG gene fusion in prostate cancer. Neoplasia. 2008, 10: 177-188. 10.1593/neo.07822
Rajput AB, Miller MA, De Luca A, Boyd N, Leung S, Hurtado-Coll A, Fazli L, Jones EC, Palmer JB, Gleave ME: Frequency of the TMPRSS2:ERG gene fusion is increased in moderate to poorly differentiated prostate cancers. 2007, 60: 1238-1243.
Klotz L: Active surveillance versus radical treatment for favorable-risk localized prostate cancer. Curr Treat Options Oncol. 2006, 7: 355-362. 10.1007/s11864-006-0003-z
True L, Coleman I, Hawley S, Huang CY, Gifford D, Coleman R, Beer TM, Gelmann E, Datta M, Mostaghel E: A molecular correlate to the Gleason grading system for prostate adenocarcinoma. Proc Natl Acad Sci USA. 2006, 103: 10991-10996. 10.1073/pnas.0603678103
Lapointe J, Li C, Higgins JP, Rijn van de M, Bair E, Montgomery K, Ferrari M, Egevad L, Rayford W, Bergerheim U: Gene expression profiling identifies clinically relevant subtypes of prostate cancer. Proc Natl Acad Sci USA. 2004, 101: 811-816. 10.1073/pnas.0304146101
Tomlins SA, Mehra R, Rhodes DR, Cao X, Wang L, Dhanasekaran SM, Kalyana-Sundaram S, Wei JT, Rubin MA, Pienta KJ: Integrative molecular concept modeling of prostate cancer progression. Nat Genet. 2007, 39: 41-51. 10.1038/ng1935
Dhanasekaran SM, Barrette TR, Ghosh D, Shah R, Varambally S, Kurachi K, Pienta KJ, Rubin MA, Chinnaiyan AM: Delineation of prognostic biomarkers in prostate cancer. Nature. 2001, 412: 822-826. 10.1038/35090585
Tusher VG, Tibshirani R, Chu G: Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci USA. 2001, 98: 5116-5121. 10.1073/pnas.091062498
Al-Shahrour F, Diaz-Uriarte R, Dopazo J: FatiGO: a web tool for finding significant associations of Gene Ontology terms with groups of genes. 2004, 20: 578-580.
Beissbarth T, Speed TP: GOstat: find statistically overrepresented Gene Ontologies within a group of genes. Bioinformatics. 2004, 20: 1464-1465. 10.1093/bioinformatics/bth088
Das S, Hahn Y, Nagata S, Willingham MC, Bera TK, Lee B, Pastan I: NGEP, a prostate-specific plasma membrane protein that promotes the association of LNCaP cells. Cancer Res. 2007, 67: 1594-1601. 10.1158/0008-5472.CAN-06-2673
Igney FH, Krammer PH: Death and anti-death: tumour resistance to apoptosis. Nat Rev Cancer. 2002, 2: 277-288. 10.1038/nrc776
Cooper CS, Campbell C, Jhavar S: Mechanisms of Disease: biomarkers and molecular targets from microarray gene expression studies in prostate cancer. Nat Clin Pract Urol. 2007, 4: 677-687. 10.1038/ncpuro0946
Takeda K, Noguchi T, Naguro I, Ichijo H: Apoptosis signal-regulating kinase 1 in stress and immune response. Annu Rev Pharmacol Toxicol. 2008, 48: 199-225. 10.1146/annurev.pharmtox.48.113006.094606
Iriyama T, Takeda K, Nakamura H, Morimoto Y, Kuroiwa T, Mizukami J, Umeda T, Noguchi T, Naguro I, Nishitoh H: ASK1 and ASK2 differentially regulate the counteracting roles of apoptosis and inflammation in tumorigenesis. Embo J. 2009, 28: 843-853. 10.1038/emboj.2009.32
Hail N, Kim HJ, Lotan R: Mechanisms of fenretinide-induced apoptosis. Apoptosis. 2006, 11: 1677-1694. 10.1007/s10495-006-9289-3
Joshi B, Li L, Taffe BG, Zhu Z, Wahl S, Tian H, Ben-Josef E, Taylor JD, Porter AT, Tang DG: Apoptosis induction by a novel anti-prostate cancer compound, BMD188 (a fatty acid-containing hydroxamic acid), requires the mitochondrial respiratory chain. Cancer Res. 1999, 59: 4343-4355.
Furuya Y, Lundmo P, Short AD, Gill DL, Isaacs JT: The role of calcium, pH, and cell proliferation in the programmed (apoptotic) death of androgen-independent prostatic cancer cells induced by thapsigargin. Cancer Res. 1994, 54: 6167-6175.
Welsh JB, Sapinoso LM, Su AI, Kern SG, Wang-Rodriguez J, Moskaluk CA, Frierson HF, Hampton GM: Analysis of gene expression identifies candidate markers and pharmacological targets in prostate cancer. Cancer Res. 2001, 61: 5974-5978.
Ichijo H, Nishida E, Irie K, ten Dijke P, Saitoh M, Moriguchi T, Takagi M, Matsumoto K, Miyazono K, Gotoh Y: Induction of apoptosis by ASK1, a mammalian MAPKKK that activates SAPK/JNK and p38 signaling pathways. Science. 1997, 275: 90-94. 10.1126/science.275.5296.90
Matsuzawa A, Saegusa K, Noguchi T, Sadamitsu C, Nishitoh H, Nagai S, Koyasu S, Matsumoto K, Takeda K, Ichijo H: ROS-dependent activation of the TRAF6-ASK1-p38 pathway is selectively required for TLR4-mediated innate immunity. Nat Immunol. 2005, 6: 587-592. 10.1038/ni1200
Dakubo GD, Parr RL, Costello LC, Franklin RB, Thayer RE: Altered metabolism and mitochondrial genome in prostate cancer. J Clin Pathol. 2006, 59: 10-16. 10.1136/jcp.2005.027664
Kumar B, Koul S, Khandrika L, Meacham RB, Koul HK: Oxidative stress is inherent in prostate cancer cells and is required for aggressive phenotype. Cancer Res. 2008, 68: 1777-1785. 10.1158/0008-5472.CAN-07-5259
Ni M, Lee AS: ER chaperones in mammalian development and human diseases. FEBS Lett. 2007, 581: 3641-3651. 10.1016/j.febslet.2007.04.045
Corazzari M, Lovat PE, Armstrong JL, Fimia GM, Hill DS, Birch-Machin M, Redfern CP, Piacentini M: Targeting homeostatic mechanisms of endoplasmic reticulum stress to increase susceptibility of cancer cells to fenretinide-induced apoptosis: the role of stress proteins ERdj5 and ERp57. Br J Cancer. 2007, 96: 1062-1071. 10.1038/sj.bjc.6603672
Leys CM, Nomura S, LaFleur BJ, Ferrone S, Kaminishi M, Montgomery E, Goldenring JR: Expression and prognostic significance of prothymosin-alpha and ERp57 in human gastric cancer. Surgery. 2007, 141: 41-50. 10.1016/j.surg.2006.05.009
Garbi N, Hammerling G, Tanaka S: Interaction of ERp57 and tapasin in the generation of MHC class I-peptide complexes. Curr Opin Immunol. 2007, 19: 99-105. 10.1016/j.coi.2006.11.013
Panaretakis T, Joza N, Modjtahedi N, Tesniere A, Vitale I, Durchschlag M, Fimia GM, Kepp O, Piacentini M, Froehlich KU: The co-translocation of ERp57 and calreticulin determines the immunogenicity of cell death. Cell Death Differ. 2008, 15: 1499-1509. 10.1038/cdd.2008.67
Garbi N, Tanaka S, Momburg F, Hammerling GJ: Impaired assembly of the major histocompatibility complex class I peptide-loading complex in mice deficient in the oxidoreductase ERp57. Nat Immunol. 2006, 7: 93-102. 10.1038/ni1288
Bartsch G, Horninger W, Klocker H, Pelzer A, Bektic J, Oberaigner W, Schennach H, Schafer G, Frauscher F, Boniol M: Tyrol Prostate Cancer Demonstration Project: early detection, treatment, outcome, incidence and mortality. BJU Int. 2008, 101: 809-816. 10.1111/j.1464-410X.2008.07502.x
Boer JM, Huber WK, Sultmann H, Wilmer F, von Heydebreck A, Haas S, Korn B, Gunawan B, Vente A, Fuzesi L: Identification and classification of differentially expressed genes in renal cell carcinoma by expression profiling on a global human 31, 500-element cDNA array. Genome Res. 2001, 11: 1861-1870.
Buness A, Huber W, Steiner K, Sultmann H, Poustka A: arrayMagic: two-colour cDNA microarray quality control and preprocessing. Bioinformatics. 2005, 21: 554-556. 10.1093/bioinformatics/bti052
Huber W, von Heydebreck A, Sultmann H, Poustka A, Vingron M: Variance stabilization applied to microarray data calibration and to the quantification of differential expression. Bioinformatics. 2002, 18 (Suppl 1): S96-104.
Bussey KJ, Kane D, Sunshine M, Narasimhan S, Nishizuka S, Reinhold WC, Zeeberg B, Ajay W, Weinstein JN: Matchminer: a tool for batch navigation among gene and gene product identities. Genome Biol. 2003, 4: R27- 10.1186/gb-2003-4-4-r27
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: 402-408. 10.1006/meth.2001.1262
This article is dedicated to the memory of Professor Annemarie Poustka, who was the founder and head of the Division Molecular Genome Analysis at the DKFZ. She was an inspiring scientist and a wonderful person.
This study was supported by a grant of the Austrian Nationalstiftung and the Austria Wirtschaftsservice GmbH in the framework of the IMGuS research program (Institute for Medical Genome Research and Systems Biology, Wien). We thank Sabrina Balaguer Marcello Schifani, Christian Pohl, Christof Seifarth and Heidi Huebl for excellent technical assistance, Holger Fröhlich for statistical analysis and Sabine Elowe for critical reading of the manuscript.
The authors declare that they have no competing interests.
NCP participated in the design and coordination of the study, contributed to assay designs and methods development, interpretation of data, and drafted the manuscript. HS and AP directed the design and coordination of the study and contributed to drafting the manuscript. HK, GS and ES identified the cancer cases, provided microdissected samples, performed immunohistochemistry and drafted the manuscript. GB participated in the design and coordination of the study. VDL carried out the microarray study, participated in the study design and methods development. RK participated in data analysis and manuscript drafting. MR performed statistical data analysis.
All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: Differentially expressed genes between GS 6 and GS 8-10 prostate cancer. A two class unpaired SAM test with 1000 permutations was performed to identify genes differentially regulated between GS 6 and GS ≥ 8 tumors. The False Discovery Rate (FDR) was set below 5%. Genes were assigned using RZPD ID, gene symbol and gene name. (PDF 58 KB)
Additional file 2: Differentially expressed genes between prostate normal and cancer tissue. A two class unpaired SAM test with 1000 permutations was performed to identify genes differentially regulated between normal and cancer tissue. The False Discovery Rate (FDR) was set below 5%. Genes were assigned using RZPD ID, gene symbol and gene name. (PDF 109 KB)
Additional file 3: Differentially expressed genes between tumors of high versus low Gleason score across independent datasets. Intersection of differentially expressed genes between high- and low risk prostate cancer with microarray data of True et al. 2006  and Lapointe et al. 2004 . (PDF 7 KB)
Additional file 4: Selected genes linked to apoptosis and differentially expressed genes between tumors of high versus low Gleason score. Gene ontology analyses (FatiGO, Ingenuity, GOstat) were performed to link differentially expressed genes between high versus low Gleason score to the gene ontology apoptosis (GO:0006915, GO:0008219). Statistical data from the microarray experiment were included. (PDF 11 KB)
Additional file 5: Genes analyzed by quantitative real-time PCR. In total, 68 genes were analyzed by qRT-PCR in order to validate expression differences between tumors of high (GS8-10) and low (GS6) Gleason scores. p-values from the statistical analysis are shown. In addition, fold change, mean Ct values of both tumor groups and the primer sequences of the gene-specific assays are indicated. (PDF 12 KB)
About this article
Cite this article
Pressinotti, N.C., Klocker, H., Schäfer, G. et al. Differential expression of apoptotic genes PDIA3 and MAP3K5 distinguishes between low- and high-risk prostate cancer. Mol Cancer 8, 130 (2009) doi:10.1186/1476-4598-8-130
- Prostate Cancer
- Gleason Score
- Prostate Intraepithelial Neoplasia
- Human Prostate Cancer Cell Line