ABSTRACT
BACKGROUND: Understanding the heterogeneous genotypes and phenotypes of prostate cancer is fundamental to improving the way we treat this disease. As yet, there are no validated descriptions of prostate cancer subgroups derived from integrated genomics linked with clinical outcome. METHODS: In a study of 482 tumour, benign and germline samples from 259 men with primary prostate cancer, we used integrative analysis of copy number alterations (CNA) and array transcriptomics to identify genomic loci that affect expression levels of mRNA in an expression quantitative trait loci (eQTL) approach, to stratify patients into subgroups that we then associated with future clinical behaviour, and compared with either CNA or transcriptomics alone. FINDINGS: We identified five separate patient subgroups with distinct genomic alterations and expression profiles based on 100 discriminating genes in our separate discovery and validation sets of 125 and 103 men. These subgroups were able to consistently predict biochemical relapse (p = 0.0017 and p = 0.016 respectively) and were further validated in a third cohort with long-term follow-up (p = 0.027). We show the relative contributions of gene expression and copy number data on phenotype, and demonstrate the improved power gained from integrative analyses. We confirm alterations in six genes previously associated with prostate cancer (MAP3K7, MELK, RCBTB2, ELAC2, TPD52, ZBTB4), and also identify 94 genes not previously linked to prostate cancer progression that would not have been detected using either transcript or copy number data alone. We confirm a number of previously published molecular changes associated with high risk disease, including MYC amplification, and NKX3-1, RB1 and PTEN deletions, as well as over-expression of PCA3 and AMACR, and loss of MSMB in tumour tissue. A subset of the 100 genes outperforms established clinical predictors of poor prognosis (PSA, Gleason score), as well as previously published gene signatures (p = 0.0001). We further show how our molecular profiles can be used for the early detection of aggressive cases in a clinical setting, and inform treatment decisions. INTERPRETATION: For the first time in prostate cancer this study demonstrates the importance of integrated genomic analyses incorporating both benign and tumour tissue data in identifying molecular alterations leading to the generation of robust gene sets that are predictive of clinical outcome in independent patient cohorts.
Subject(s)
Gene Dosage , Prostatic Neoplasms/genetics , Transcriptome/genetics , Adult , Aged , Aged, 80 and over , Cluster Analysis , Cohort Studies , Gene Expression Profiling , Gene Expression Regulation, Neoplastic , Genome, Human , Humans , Male , Middle Aged , Prognosis , RNA, Messenger/genetics , RNA, Messenger/metabolism , Recurrence , Reproducibility of Results , Risk FactorsABSTRACT
Illumina's GoldenGate technology is a two-channel microarray platform that allows for the simultaneous interrogation of about 1,500 locations in the genome. GoldenGate has proved a flexible platform not only in the choice of those 1,500 locations, but also in the choice of the property being measured at them. It retains the desirable properties of Illumina's BeadArrays in that the probes (in this case 'beads') are randomly arranged across the microarray, there are multiple instances of each probe and many samples can be processed simultaneously. As for other Illumina technologies, however, these properties are not exploited as they might be. Here we review the various common adaptations of the GoldenGate platform, review the analysis methods that are associated with each adaptation and then, with the aid of a number of example data sets we illustrate some of the improvements that can be made over the default analysis.
Subject(s)
Data Interpretation, Statistical , Genomics , Oligonucleotide Array Sequence Analysis/instrumentation , Oligonucleotide Array Sequence Analysis/methods , Alleles , Color , DNA Methylation , Gene Expression , Gene Expression Profiling , Genotype , HumansABSTRACT
SUMMARY: With their many replicates and their random layouts, Illumina BeadArrays provide greater scope fordetecting spatial artefacts than do other microarray technologies. They are also robust to artefact exclusion, yet there is a lack of tools that can perform these tasks for Illumina. We present BASH, a tool for this purpose. BASH adopts the concepts of Harshlight, but implements them in a manner that utilizes the unique characteristics of the Illumina technology. Using bead-level data, spatial artefacts of various kinds can thus be identified and excluded from further analyses. AVAILABILITY: The beadarray Bioconductor package (version 1.10 onwards), www.bioconductor.org
