| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
Priority Reports |
1 Medical Biotechnology, VTT Technical Research Centre of Finland and University of Turku, Turku, Finland; 2 Department of Cancer Prevention, Institute for Cancer Research, Rikshospitalet-Radiumhospitalet Medical Center, Oslo, Norway; and 3 Department of Urology, Radboud University Nijmegen Medical Center, the Netherlands
Requests for reprints: Olli Kallioniemi, Medical Biotechnology Centre, VTT Technical Research Centre of Finland and University of Turku, Turku, Finland. Phone: 358-2478-8602; Fax: 358-2478-8601; E-mail: Olli.Kallioniemi{at}vtt.fi.
| Abstract |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
| Materials and Methods |
|---|
|
|
|---|
Sample preparation, gene copy number, and expression data. DNA was extracted after overnight proteinase K treatment using standard protocols, whereas total RNA was extracted using either TRIzol (Invitrogen, Carlsbad, CA) or LiCl/urea isolationbased methods. Genome-wide data for DNA copy number and matching gene expression data from the same samples will be described in detail elsewhere.4 The data were merged to analyze the corresponding copy number and expression levels for the ERG, ETV4, ETV1, and TMPRSS2 loci. Briefly, array-based CGH was done using human genome CGH 44A and 44B oligo microarrays according to the version 2 protocol provided by Agilent Technologies (Palo Alto, CA), with minor modifications. Male genomic DNA (Promega, Madison, WI) was used as a reference in all hybridizations. Three micrograms of digested tumor DNA and reference DNA were labeled using Cy5-dUTP and Cy3-dUTP incorporation (Perkin-Elmer, Wellesley, MA) and the Bioprime Array CGH genomic labeling module (Invitrogen). A laser confocal scanner (Agilent Technologies) was used to obtain signal intensities from targets, and Agilent Feature Extraction software (version 8.1.1.1) was applied using the manufacturer's recommended settings. To analyze the aCGH data, we used the CGH Analytics software (version 3.2.32; Agilent Technologies).
Gene expression levels were measured using Affymetrix GeneChip Human Genome U133 Plus 2.0 arrays (Affymetrix, Santa Clara, CA). Sample processing and labeling were done according to the protocol provided by Affymetrix. Three micrograms of total RNA from each sample was used for the initial one-cycle cDNA synthesis. Arrays were scanned immediately after staining using a GeneChip scanner (Affymetrix).
Reverse transcription-PCR to identify TMPRSS2-ERG, TMPRSS2-ETV1, and TMPRSS2-ETV4 fusion transcripts. Reverse transcription used 50 ng of total RNA, Sensiscript reverse transcription kit (Qiagen, Valencia, CA) and oligo dT-primers. The fusion transcripts were PCR-amplified, using gene-specific primers for TMPRSS2 (5'-CAGAGCTGCTAACAGGAGGCGGAGGCGGA-3'), ERG (5'-CATAGTAGTAACGGAGGGCGC-3'), ETV1 (5'-TTGTGGTGGGAAGGGGATGTTT-3'), and ETV4 (5'-CGAAGTCCGTCTGTTCCTGT-3') cDNAs. The PCR was done with Phusion High-Fidelity DNA polymerase (Finnzymes, Espoo, Finland). All PCR experiments included reverse transcription (RT)-negative controls and a blank with no template. PCR products were isolated from agarose gels, treated with Taq polymerase to generate polyadenylic acid overhangs, and cloned into pCRII-TOPO cloning vector (Invitrogen). Sequencing reactions using the same primers as for amplification were prepared by using the ABI BigDye Terminator V3.1 cycle sequencing kit, according to the manufacturer's instructions and analyzed on the ABI 3100 genetic Analyzer (Applied Biosystems, Foster City, CA).
Data filtering and normalization, clustering, and statistical analysis. Affymetrix U133 Plus 2.0 arrays were normalized using R (8) and the RMA (9) implementation in Bioconductor package affy. Multiple probe sets mapping to the same genes were combined using mean values. Both genes and samples were clustered hierarchically using Euclidean distance and complete linkage analyses.
In silico analysis of potential ERG target genes. We assembled expression data from 410 human prostate tissue samples, consisting of 178 normal samples and 232 tumors and metastases (Supplementary Table S2). We analyzed the patterns of ERG coexpressed genes by four independent methods. First, significance analysis of microarrays (SAM; ref. 10) was done to identify genes that correlate with ERG in prostate tissues. Prostate samples from the HG-U95A platform were divided into ERG-positive samples (group 2, n = 16), expressing ERG at higher levels than any of the normal samples, and into ERG-negative samples including both normal and tumor tissues (group 1, n = 93). Normalized gene expression data for samples used in SAM are provided in Supplementary Table S3. Another SAM analysis was done with ERG-positive (n = 16) versus ERG-negative (n = 48) tumor samples only. The false discovery rate was set to zero in both analyses. Second, hierarchical clustering analysis was done to identify the genes coexpressed with ERG in prostate samples in an unsupervised manner. Third, we identified genes whose expression was most closely associated with that of ERG's by using a Perl implementation of Pearson's correlation (correlation factors >0.5 or
0.5). Finally, Gene Ontology analysis using the DAVID GO analysis tool (http://www.niaid.abcc.ncifcrf.gov) and gene set enrichment analysis (Broad Institute of MIT and Harvard) were done using the same expression data as in the SAM analyses. The data from the three different patient cohorts and the four different analysis methods were overlaid to define the most consistent alterations associated with ERG gene expression.
| Results and Discussion |
|---|
|
|
|---|
|
40% of advanced prostate cancers. Six of these samples (nos. 3, 4, 7, 13, 14, and 16) were strongly positive, whereas in one advanced prostate cancer sample (no. 10), the fusion transcript was expressed at very low levels (data not shown). Sequence analysis of the RT-PCR products indicate that the most common fusion transcript was a fusion of exon 1 of TMPRSS2 with the beginning of exon 4 of ERG (TMPRSS2:ERGa, ref. 1). In one tumor (no. 10), exon 1 of TMPRSS2 was fused to the beginning of exon 2 of ERG (TMPRSS2:ERGb, ref. 1). In sample no. 16, exon 3 of TMPRSS2 was fused to the beginning of exon 4 of ERG. In tumor no. 10, there was also a more complicated fusion consisting of exon 2 of TMPRSS2 fused to 95 nucleotides identical to a segment of ERG splice form 6 (ERG6 mRNA, AY204740.1; bp, 225-286), followed by the entire exon 4 of ERG. The schematic presentation of the novel TMPRSS2:ERG fusion transcript is presented in Supplementary Fig. S1. In summary, the sequence analysis of TMPRSS2:ERG fusion transcripts indicated significant variation in fusion breakpoints. Genomic rearrangements at the ERG locus by array-CGH analysis. Genomic rearrangements, either deletions at the ERG locus or interstitial deletions between the TMPRSS2 and ERG loci, were identified by array CGH in five out of the six samples displaying TMPRSS2:ERG gene fusions with high ERG expression (nos. 3, 4, 13, 14, and 16; Fig. 2 ). This indicates that the ERG activation is not caused by simple balanced translocation, but by a variety of unbalanced genetic rearrangements that bring together these two adjacent loci. The proximity of the two genes (with a genomic distance of only 2.8 Mb), and location in the same DNA strand, may facilitate the fusion gene formation and allow a simple intergenic deletion to activate the ERG gene by the fusion to TMPRSS2. This was the most prevalent genetic alteration by array-CGH. In one of the tumors (no. 4), there were two microdeletions both at the TMPRSS2 and ERG loci (measuring 911 and 159 kb, respectively), suggesting that complex genetic rearrangements led to fusion gene formation.
|
In silico analysis of ERG coexpressed genes in prostate tumors. The ERG gene is the most consistently overexpressed oncogene in malignant epithelial cells of the prostate (4). However, its functional role in prostate cancer development and progression has not yet been clearly determined. To identify potential ERG target genes and deregulated biological processes in vivo, we analyzed gene coexpression data from three different prostate data sets, altogether consisting of 410 human prostate tissue samples.
The largest data set (n = 284) consisted of nine previously published Affymetrix gene expression studies. These data were normalized to render them directly comparable.5 Expression data were analyzed by multiple statistical methods to characterize the most consistent ERG-associated genes. First, the results from SAM analysis indicated that 136 genes were significantly differentially expressed between ERG-positive and ERG-negative prostate samples; 92 genes were positively correlated (>1.5-fold change, positive hits), and 44 genes were negatively correlated with ERG (Supplementary Table S5). SAM was also done with prostate cancers only (Fig. 3 ). Second, the ERG clusters generated by K means and hierarchical clustering showed a 46% and 59% overlap, respectively, with the positive SAM hit list. Third, ERG correlation analyses across all prostate samples (n = 284) or all prostate cancer samples (n = 147) were done. For the top 200 genes correlating with ERG in prostate cancers, there was an 85% overlap with the SAM-positive hits using linear correlation and a 75% overlap for log-transformed correlation.
|
The third data set (n = 14) consisted of the advanced prostate cancers analyzed in this study. Although the number of samples was small and many samples were derived from hormone-refractory prostate cancers, three genes among the top ERG-correlating genes were also among the top 10 positive SAM hit gene list. The correlation results from the three different data sets overlaid with the SAM results are presented in Supplementary Table S5.
Histone deacetylase 1 (HDAC1) was the only gene among the top ERG coexpressed genes in all three data sets and therefore was the most consistent feature of ERG-overexpressing prostate cancers. By RT-PCR validation of the HDAC1 expression levels, all ERG-positive prostate cancers were strongly HDAC1-positive, whereas the ERG-negative tumors showed more variable, but significantly lower (P = 0.0001) expression levels (Supplementary Fig. S2). HDAC1 catalyzes the deacetylation of lysine residues of the core histones and other proteins, leading to epigenetic silencing of target genes. HDAC1 has been previously shown to be strongly expressed in hormone-refractory prostate cancers (13). Currently, it remains unclear if HDAC1 is a direct transcriptional target of ERG or whether HDAC1 up-regulation results from other changes occurring in ERG-overexpressing tumors. ERG has been shown to interact indirectly with HDAC1 via SETDB1 methyl transferase (14, 15).
Results from gene ontology analyses indicated that the "organogenesis" as well as "cell growth and maintenance" were the most significantly (P = 0.001 and 0.02) overrepresented gene ontology terms in ERG-positive tumors. Gene set enrichment analysis results, presented in Fig. 4 , indicated that the WNT and PITX2 pathways were among the most highly enriched pathways in ERG-overexpressing tumors (16, 17). The WNT pathway controls organogenesis by inducing, e.g., PITX2 transcription factor, which serves as an important modulator of growth control genes (17). HDAC1 itself is linked to these two pathways. The down-regulated gene sets in ERG-overexpressing tumors included CCR5, cell death, and tumor necrosis factor/FAS. Interestingly, the HDAC pathway with known HDAC target genes and regulators was also highlighted by this analysis, suggesting that the up-regulation of HDAC1 in ERG-positive tumors led, as could be expected, to the down-regulation of HDAC target genes (18, 19). The genes showing core enrichment in the identified pathways are presented in Supplementary Table S6.
|
| Acknowledgments |
|---|
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
| Footnotes |
|---|
4 M. Wolf et al., in preparation. ![]()
5 Autio et al., in preparation. ![]()
Received 5/31/06. Revised 9/ 5/06. Accepted 9/21/06.
| References |
|---|
|
|
|---|
Pitx2 pathway mediating cell-type-specific proliferation during development. Cell 2002;111:67385.[CrossRef][Medline]This article has been cited by other articles:
![]() |
G Attard, J E Ang, D Olmos, and J S de Bono Dissecting prostate carcinogenesis through ETS gene rearrangement studies: implications for anticancer drug development J. Clin. Pathol., August 1, 2008; 61(8): 891 - 896. [Abstract] [Full Text] [PDF] |
||||
![]() |
J.-M. Mosquera, S. Perner, E. M. Genega, M. Sanda, M. D. Hofer, K. D. Mertz, P. L. Paris, J. Simko, T. A. Bismar, G. Ayala, et al. Characterization of TMPRSS2-ERG Fusion High-Grade Prostatic Intraepithelial Neoplasia and Potential Clinical Implications Clin. Cancer Res., June 1, 2008; 14(11): 3380 - 3385. [Abstract] [Full Text] [PDF] |
||||
![]() |
O. R. Saramaki, A. E. Harjula, P. M. Martikainen, R. L. Vessella, T. L.J. Tammela, and T. Visakorpi TMPRSS2:ERG Fusion Identifies a Subgroup of Prostate Cancers with a Favorable Prognosis Clin. Cancer Res., June 1, 2008; 14(11): 3395 - 3400. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Mehra, S. A. Tomlins, J. Yu, X. Cao, L. Wang, A. Menon, M. A. Rubin, K. J. Pienta, R. B. Shah, and A. M. Chinnaiyan Characterization of TMPRSS2-ETS Gene Aberrations in Androgen-Independent Metastatic Prostate Cancer Cancer Res., May 15, 2008; 68(10): 3584 - 3590. [Abstract] [Full Text] [PDF] |
||||
![]() |
O. Klezovitch, M. Risk, I. Coleman, J. M. Lucas, M. Null, L. D. True, P. S. Nelson, and V. Vasioukhin A causal role for ERG in neoplastic transformation of prostate epithelium PNAS, February 12, 2008; 105(6): 2105 - 2110. [Abstract] [Full Text] [PDF] |
||||
![]() |
B. E. Helgeson, S. A. Tomlins, N. Shah, B. Laxman, Q. Cao, J. R. Prensner, X. Cao, N. Singla, J. E. Montie, S. Varambally, et al. Characterization of TMPRSS2:ETV5 and SLC45A3:ETV5 Gene Fusions in Prostate Cancer Cancer Res., January 1, 2008; 68(1): 73 - 80. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. Jhavar, A. Reid, J. Clark, Z. Kote-Jarai, T. Christmas, A. Thompson, C. Woodhouse, C. Ogden, C. Fisher, C. Corbishley, et al. Detection of TMPRSS2-ERG Translocations in Human Prostate Cancer by Expression Profiling Using GeneChip Human Exon 1.0 ST Arrays J. Mol. Diagn., January 1, 2008; 10(1): 50 - 57. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. Demichelis and M. A Rubin TMPRSS2-ETS fusion prostate cancer: biological and clinical implications J. Clin. Pathol., November 1, 2007; 60(11): 1185 - 1186. [Full Text] [PDF] |
||||
![]() |
A. B Rajput, M. A Miller, A. De Luca, N. Boyd, S. Leung, A. Hurtado-Coll, L. Fazli, E. C Jones, J. B Palmer, M. E Gleave, et al. Frequency of the TMPRSS2:ERG gene fusion is increased in moderate to poorly differentiated prostate cancers J. Clin. Pathol., November 1, 2007; 60(11): 1238 - 1243. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| Cancer Research | Clinical Cancer Research |
| Cancer Epidemiology Biomarkers & Prevention | Molecular Cancer Therapeutics |
| Molecular Cancer Research | Cancer Prevention Research |
| Cancer Prevention Journals Portal | Cancer Reviews Online |
| Annual Meeting Education Book | Meeting Abstracts Online |