Genomic evolution shapes prostate cancer disease type.

The development of cancer is an evolutionary process involving the sequential acquisition of genetic alterations that disrupt normal biological processes, enabling tumor cells to rapidly proliferate and eventually invade and metastasize to other tissues. We investigated the genomic evolution of prostate cancer through the application of three separate classification methods, each designed to investigate a different aspect of tumor evolution. Integrating the results revealed the existence of two distinct types of prostate cancer that arise from divergent evolutionary trajectories, designated as the Canonical and Alternative evolutionary disease types. We therefore propose the evotype model for prostate cancer evolution wherein Alternative-evotype tumors diverge from those of the Canonical-evotype through the stochastic accumulation of genetic alterations associated with disruptions to androgen receptor DNA binding. Our model unifies many previous molecular observations, providing a powerful new framework to investigate prostate cancer disease progression.

Estimation of tumor cell total mRNA expression in 15 cancer types predicts disease progression.

Single-cell RNA sequencing studies have suggested that total mRNA content correlates with tumor phenotypes. Technical and analytical challenges, however, have so far impeded at-scale pan-cancer examination of total mRNA content. Here we present a method to quantify tumor-specific total mRNA expression (TmS) from bulk sequencing data, taking into account tumor transcript proportion, purity and ploidy, which are estimated through transcriptomic/genomic deconvolution. We estimate and validate TmS in 6,590 patient tumors across 15 cancer types, identifying significant inter-tumor variability. Across cancers, high TmS is associated with increased risk of disease progression and death. TmS is influenced by cancer-specific patterns of gene alteration and intra-tumor genetic heterogeneity as well as by pan-cancer trends in metabolic dysregulation. Taken together, our results indicate that measuring cell-type-specific total mRNA expression in tumor cells predicts tumor phenotypes and clinical outcomes.

Quantifying the influence of mutation detection on tumour subclonal reconstruction.

Whole-genome sequencing can be used to estimate subclonal populations in tumours and this intra-tumoural heterogeneity is linked to clinical outcomes. Many algorithms have been developed for subclonal reconstruction, but their variabilities and consistencies are largely unknown. We evaluate sixteen pipelines for reconstructing the evolutionary histories of 293 localized prostate cancers from single samples, and eighteen pipelines for the reconstruction of 10 tumours with multi-region sampling. We show that predictions of subclonal architecture and timing of somatic mutations vary extensively across pipelines. Pipelines show consistent types of biases, with those incorporating SomaticSniper and Battenberg preferentially predicting homogenous cancer cell populations and those using MuTect tending to predict multiple populations of cancer cells. Subclonal reconstructions using multi-region sampling confirm that single-sample reconstructions systematically underestimate intra-tumoural heterogeneity, predicting on average fewer than half of the cancer cell populations identified by multi-region sequencing. Overall, these biases suggest caution in interpreting specific architectures and subclonal variants.

Divergent mutational processes distinguish hypoxic and normoxic tumours.

Many primary tumours have low levels of molecular oxygen (hypoxia), and hypoxic tumours respond poorly to therapy. Pan-cancer molecular hallmarks of tumour hypoxia remain poorly understood, with limited comprehension of its associations with specific mutational processes, non-coding driver genes and evolutionary features. Here, as part of the ICGC/TCGA Pan-Cancer Analysis of Whole Genomes (PCAWG) Consortium, which aggregated whole genome sequencing data from 2658 cancers across 38 tumour types, we quantify hypoxia in 1188 tumours spanning 27 cancer types. Elevated hypoxia associates with increased mutational load across cancer types, irrespective of underlying mutational class. The proportion of mutations attributed to several mutational signatures of unknown aetiology directly associates with the level of hypoxia, suggesting underlying mutational processes for these signatures. At the gene level, driver mutations in TP53, MYC and PTEN are enriched in hypoxic tumours, and mutations in PTEN interact with hypoxia to direct tumour evolutionary trajectories. Overall, hypoxia plays a critical role in shaping the genomic and evolutionary landscapes of cancer.

Widespread and Functional RNA Circularization in Localized Prostate Cancer.

The cancer transcriptome is remarkably complex, including low-abundance transcripts, many not polyadenylated. To fully characterize the transcriptome of localized prostate cancer, we performed ultra-deep total RNA-seq on 144 tumors with rich clinical annotation. This revealed a linear transcriptomic subtype associated with the aggressive intraductal carcinoma sub-histology and a fusion profile that differentiates localized from metastatic disease. Analysis of back-splicing events showed widespread RNA circularization, with the average tumor expressing 7,232 circular RNAs (circRNAs). The degree of circRNA production was correlated to disease progression in multiple patient cohorts. Loss-of-function screening identified 11.3% of highly abundant circRNAs as essential for cell proliferation; for ∼90% of these, their parental linear transcripts were not essential. Individual circRNAs can have distinct functions, with circCSNK1G3 promoting cell growth by interacting with miR-181. These data advocate for adoption of ultra-deep RNA-seq without poly-A selection to interrogate both linear and circular transcriptomes.

Molecular landmarks of tumor hypoxia across cancer types.

Many primary-tumor subregions have low levels of molecular oxygen, termed hypoxia. Hypoxic tumors are at elevated risk for local failure and distant metastasis, but the molecular hallmarks of tumor hypoxia remain poorly defined. To fill this gap, we quantified hypoxia in 8,006 tumors across 19 tumor types. In ten tumor types, hypoxia was associated with elevated genomic instability. In all 19 tumor types, hypoxic tumors exhibited characteristic driver-mutation signatures. We observed widespread hypoxia-associated dysregulation of microRNAs (miRNAs) across cancers and functionally validated miR-133a-3p as a hypoxia-modulated miRNA. In localized prostate cancer, hypoxia was associated with elevated rates of chromothripsis, allelic loss of PTEN and shorter telomeres. These associations are particularly enriched in polyclonal tumors, representing a constellation of features resembling tumor nimbosus, an aggressive cellular phenotype. Overall, this work establishes that tumor hypoxia may drive aggressive molecular features across cancers and shape the clinical trajectory of individual tumors.

Development and Validation of a 28-gene Hypoxia-related Prognostic Signature for Localized Prostate Cancer.

BACKGROUND: Hypoxia is associated with a poor prognosis in prostate cancer. This work aimed to derive and validate a hypoxia-related mRNA signature for localized prostate cancer. METHOD: Hypoxia genes were identified in vitro via RNA-sequencing and combined with in vivo gene co-expression analysis to generate a signature. The signature was independently validated in eleven prostate cancer cohorts and a bladder cancer phase III randomized trial of radiotherapy alone or with carbogen and nicotinamide (CON). RESULTS: A 28-gene signature was derived. Patients with high signature scores had poorer biochemical recurrence free survivals in six of eight independent cohorts of prostatectomy-treated patients (Log rank test P < .05), with borderline significances achieved in the other two (P < .1). The signature also predicted biochemical recurrence in patients receiving post-prostatectomy radiotherapy (n = 130, P = .007) or definitive radiotherapy alone (n = 248, P = .035). Lastly, the signature predicted metastasis events in a pooled cohort (n = 631, P = .002). Prognostic significance remained after adjusting for clinic-pathological factors and commercially available prognostic signatures. The signature predicted benefit from hypoxia-modifying therapy in bladder cancer patients (intervention-by-signature interaction test P = .0026), where carbogen and nicotinamide was associated with improved survival only in hypoxic tumours. CONCLUSION: A 28-gene hypoxia signature has strong and independent prognostic value for prostate cancer patients.

The Evolutionary Landscape of Localized Prostate Cancers Drives Clinical Aggression.

The majority of newly diagnosed prostate cancers are slow growing, with a long natural life history. Yet a subset can metastasize with lethal consequences. We reconstructed the phylogenies of 293 localized prostate tumors linked to clinical outcome data. Multiple subclones were detected in 59% of patients, and specific subclonal architectures associate with adverse clinicopathological features. Early tumor development is characterized by point mutations and deletions followed by later subclonal amplifications and changes in trinucleotide mutational signatures. Specific genes are selectively mutated prior to or following subclonal diversification, including MTOR, NKX3-1, and RB1. Patients with low-risk monoclonal tumors rarely relapse after primary therapy (7%), while those with high-risk polyclonal tumors frequently do (61%). The presence of multiple subclones in an index biopsy may be necessary, but not sufficient, for relapse of localized prostate cancer, suggesting that evolution-aware biomarkers should be studied in prospective studies of low-risk tumors suitable for active surveillance.

A Prostate Cancer "Nimbosus": Genomic Instability and SChLAP1 Dysregulation Underpin Aggression of Intraductal and Cribriform Subpathologies.

BACKGROUND: Intraductal carcinoma (IDC) and cribriform architecture (CA) represent unfavorable subpathologies in localized prostate cancer. We recently showed that IDC shares a clonal ancestry with the adjacent glandular adenocarcinoma. OBJECTIVE: We investigated for the co-occurrence of "aggression" factors, genomic instability and hypoxia, and performed gene expression profiling of these tumors. DESIGN, SETTING, AND PARTICIPANTS: A total of 1325 men were treated for localized prostate cancer from four academic institutions (University Health Network, CHU de Québec-Université Laval, Memorial Sloan Kettering Cancer Center [MSKCC], and Erasmus Medical Center). Pathological specimens were centrally reviewed. Gene copy number and expression, and intraprostatic oxygenation were assessed. OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS: IDC/CA was separately assessed for biochemical relapse risk in the Canadian and MSKCC cohorts. Both cohorts were pooled for analyses on metastasis. RESULTS AND LIMITATION: Presence of IDC/CA independently predicted for increased risks of biochemical relapse (HRCanadian 2.17, p<0.001; HRMSKCC 2.32, p=0.0035) and metastasis (HRpooled 3.31, p<0.001). IDC/CA+ cancers were associated with an increased percentage of genome alteration (PGA [median] 7.2 vs 3.0, p<0.001), and hypoxia (64.0% vs 45.5%, p=0.17). Combinatorial genomic-pathological indices offered the strongest discrimination for metastasis (C-index 0.805 [clinical+IDC/CA+PGA] vs 0.786 [clinical+IDC/CA] vs 0.761 [clinical]). Profiling of mRNA abundance revealed that long noncoding RNA, SChLAP1, was the only gene expressed at >3-fold higher (p<0.0001) in IDC/CA+ than in IDC/CA- tumors, independently corroborated by increased SChLAP1 RNA in situ hybridization signal. Optimal treatment intensification for IDC/CA+ prostate cancer requires prospective testing. CONCLUSIONS: The poor outcome associated with IDC and CA subpathologies is associated with a constellation of genomic instability, SChLAP1 expression, and hypoxia. We posit a novel concept in IDC/CA+ prostate cancer, "nimbosus" (gathering of stormy clouds, Latin), which manifests as increased metastatic capacity and lethality. PATIENT SUMMARY: A constellation of unfavorable molecular characteristics co-occur with intraductal and cribriform subpathologies in prostate cancer. Modern imaging for surveillance and treatment intensification trials should be considered in this adverse subgroup.

Comparing continuous and discrete analyses of breast cancer survival information.

Treatment of cancer is becoming increasingly personalized and biomarkers continue to be developed to refine treatment decisions. Tumour mRNA abundance data is commonly used to develop such biomarkers, often to predict patient survival. However, survival analyses present unique challenges and it is unknown whether analysing mRNA abundance information in a discrete or continuous manner yields different results. To address this, we analysed 1988 primary breast tumour transcriptomes. When compared univariately, approximately 60% of all genes showed differences between the discrete and continuous Cox proportional hazards models with q-value differences spanning four orders of magnitude for some genes. Further, hybrid models using both continuous and discrete data used to classify poor prognosis via random forest outperformed models using a single type of information. Thus some genes appear to continuously contribute to poor prognosis while others display threshold effects, and incorporating this into biomarker development is a key unexplored avenue.