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.

A collection of breast cancer cell lines for the study of functionally distinct cancer subtypes.

Recent studies suggest that thousands of genes may contribute to breast cancer pathophysiologies when deregulated by genomic or epigenomic events. Here, we describe a model "system" to appraise the functional contributions of these genes to breast cancer subsets. In general, the recurrent genomic and transcriptional characteristics of 51 breast cancer cell lines mirror those of 145 primary breast tumors, although some significant differences are documented. The cell lines that comprise the system also exhibit the substantial genomic, transcriptional, and biological heterogeneity found in primary tumors. We show, using Trastuzumab (Herceptin) monotherapy as an example, that the system can be used to identify molecular features that predict or indicate response to targeted therapies or other physiological perturbations.

Genomic and transcriptional aberrations linked to breast cancer pathophysiologies.

This study explores the roles of genome copy number abnormalities (CNAs) in breast cancer pathophysiology by identifying associations between recurrent CNAs, gene expression, and clinical outcome in a set of aggressively treated early-stage breast tumors. It shows that the recurrent CNAs differ between tumor subtypes defined by expression pattern and that stratification of patients according to outcome can be improved by measuring both expression and copy number, especially high-level amplification. Sixty-six genes deregulated by the high-level amplifications are potential therapeutic targets. Nine of these (FGFR1, IKBKB, ERBB2, PROCC, ADAM9, FNTA, ACACA, PNMT, and NR1D1) are considered druggable. Low-level CNAs appear to contribute to cancer progression by altering RNA and cellular metabolism.

Integrated genome and transcriptome sequencing identifies a novel form of hybrid and aggressive prostate cancer.

Next-generation sequencing is making sequence-based molecular pathology and personalized oncology viable. We selected an individual initially diagnosed with conventional but aggressive prostate adenocarcinoma and sequenced the genome and transcriptome from primary and metastatic tissues collected prior to hormone therapy. The histology-pathology and copy number profiles were remarkably homogeneous, yet it was possible to propose the quadrant of the prostate tumour that likely seeded the metastatic diaspora. Despite a homogeneous cell type, our transcriptome analysis revealed signatures of both luminal and neuroendocrine cell types. Remarkably, the repertoire of expressed but apparently private gene fusions, including C15orf21:MYC, recapitulated this biology. We hypothesize that the amplification and over-expression of the stem cell gene MSI2 may have contributed to the stable hybrid cellular identity. This hybrid luminal-neuroendocrine tumour appears to represent a novel and highly aggressive case of prostate cancer with unique biological features and, conceivably, a propensity for rapid progression to castrate-resistance. Overall, this work highlights the importance of integrated analyses of genome, exome and transcriptome sequences for basic tumour biology, sequence-based molecular pathology and personalized oncology.

From sequence to molecular pathology, and a mechanism driving the neuroendocrine phenotype in prostate cancer.

The current paradigm of cancer care relies on predictive nomograms which integrate detailed histopathology with clinical data. However, when predictions fail, the consequences for patients are often catastrophic, especially in prostate cancer where nomograms influence the decision to therapeutically intervene. We hypothesized that the high dimensional data afforded by massively parallel sequencing (MPS) is not only capable of providing biological insights, but may aid molecular pathology of prostate tumours. We assembled a cohort of six patients with high-risk disease, and performed deep RNA and shallow DNA sequencing in primary tumours and matched metastases where available. Our analysis identified copy number abnormalities, accurately profiled gene expression levels, and detected both differential splicing and expressed fusion genes. We revealed occult and potentially dormant metastases, unambiguously supporting the patients' clinical history, and implicated the REST transcriptional complex in the development of neuroendocrine prostate cancer, validating this finding in a large independent cohort. We massively expand on the number of novel fusion genes described in prostate cancer; provide fresh evidence for the growing link between fusion gene aetiology and gene expression profiles; and show the utility of fusion genes for molecular pathology. Finally, we identified chromothripsis in a patient with chronic prostatitis. Our results provide a strong foundation for further development of MPS-based molecular pathology.

Poly-gene fusion transcripts and chromothripsis in prostate cancer.

Complex genome rearrangements are frequently observed in cancer but their impact on tumor molecular biology is largely unknown. Recent studies have identified a new phenomenon involving the simultaneous generation of tens to hundreds of genomic rearrangements, called chromothripsis. To understand the molecular consequences of these events, we sequenced the genomes and transcriptomes of two prostate tumors exhibiting evidence of chromothripsis. We identified several complex fusion transcripts, each containing sequence from three different genes, originating from different parts of the genome. One such poly-gene fusion transcript appeared to be expressed from a chain of small genomic fragments. Furthermore, we detected poly-gene fusion transcripts in the prostate cancer cell line LNCaP, suggesting they may represent a common phenomenon. Finally in one tumor with chromothripsis, we identified multiple mutations in the p53 signaling pathway, expanding on recent work associating aberrant DNA damage response mechanisms with chromothripsis. Overall, our data show that chromothripsis can manifest as massively rearranged transcriptomes. The implication that multigenic changes can give rise to poly-gene fusion transcripts is potentially of great significance to cancer genetics.

Optimally discriminative subnetwork markers predict response to chemotherapy.

MOTIVATION: Molecular profiles of tumour samples have been widely and successfully used for classification problems. A number of algorithms have been proposed to predict classes of tumor samples based on expression profiles with relatively high performance. However, prediction of response to cancer treatment has proved to be more challenging and novel approaches with improved generalizability are still highly needed. Recent studies have clearly demonstrated the advantages of integrating protein-protein interaction (PPI) data with gene expression profiles for the development of subnetwork markers in classification problems. RESULTS: We describe a novel network-based classification algorithm (OptDis) using color coding technique to identify optimally discriminative subnetwork markers. Focusing on PPI networks, we apply our algorithm to drug response studies: we evaluate our algorithm using published cohorts of breast cancer patients treated with combination chemotherapy. We show that our OptDis method improves over previously published subnetwork methods and provides better and more stable performance compared with other subnetwork and single gene methods. We also show that our subnetwork method produces predictive markers that are more reproducible across independent cohorts and offer valuable insight into biological processes underlying response to therapy. AVAILABILITY: The implementation is available at: http://www.cs.sfu.ca/~pdao/personal/OptDis.html CONTACT: cenk@cs.sfu.ca; alapuk@prostatecentre.com; ccollins@prostatecentre.com.

Comrad: detection of expressed rearrangements by integrated analysis of RNA-Seq and low coverage genome sequence data.

MOTIVATION: Comrad is a novel algorithmic framework for the integrated analysis of RNA-Seq and whole genome shotgun sequencing (WGSS) data for the purposes of discovering genomic rearrangements and aberrant transcripts. The Comrad framework leverages the advantages of both RNA-Seq and WGSS data, providing accurate classification of rearrangements as expressed or not expressed and accurate classification of the genomic or non-genomic origin of aberrant transcripts. A major benefit of Comrad is its ability to accurately identify aberrant transcripts and associated rearrangements using low coverage genome data. As a result, a Comrad analysis can be performed at a cost comparable to that of two RNA-Seq experiments, significantly lower than an analysis requiring high coverage genome data. RESULTS: We have applied Comrad to the discovery of gene fusions and read-throughs in prostate cancer cell line C4-2, a derivative of the LNCaP cell line with androgen-independent characteristics. As a proof of concept, we have rediscovered in the C4-2 data 4 of the 6 fusions previously identified in LNCaP. We also identified six novel fusion transcripts and associated genomic breakpoints, and verified their existence in LNCaP, suggesting that Comrad may be more sensitive than previous methods that have been applied to fusion discovery in LNCaP. We show that many of the gene fusions discovered using Comrad would be difficult to identify using currently available techniques. AVAILABILITY: A C++ and Perl implementation of the method demonstrated in this article is available at http://compbio.cs.sfu.ca/.

The RAB25 small GTPase determines aggressiveness of ovarian and breast cancers.

High-density array comparative genomic hybridization (CGH) showed amplification of chromosome 1q22 centered on the RAB25 small GTPase, which is implicated in apical vesicle trafficking, in approximately half of ovarian and breast cancers. RAB25 mRNA levels were selectively increased in stage III and IV serous epithelial ovarian cancers compared to other genes within the amplified region, implicating RAB25 as a driving event in the development of the amplicon. Increased DNA copy number or RNA level of RAB25 was associated with markedly decreased disease-free survival or overall survival in ovarian and breast cancers, respectively. Forced expression of RAB25 markedly increased anchorage-dependent and anchorage-independent cell proliferation, prevented apoptosis and anoikis, including that induced by chemotherapy, and increased aggressiveness of cancer cells in vivo. The inhibition of apoptosis was associated with a decrease in expression of the proapoptotic molecules, BAK and BAX, and activation of the antiapoptotic phosphatidylinositol 3 kinase (PI3K) and AKT pathway, providing potential mechanisms for the effects of RAB25 on tumor aggressiveness. Overall, these studies implicate RAB25, and thus the RAB family of small G proteins, in aggressiveness of epithelial cancers.

Prevalence of the EGFR T790M and other resistance mutations in the Australian population and histopathological correlation in a small subset of cases.

We sought to review the prevalence of EGFR T790M and other EGFR mutations associated with either proven or probable tyrosine kinase inhibitor (TKI) resistance in the Australasian lung cancer population and to perform histopathological correlation in a subset of cases. Retrospective statistical analysis was performed on a set of targeted lung cancer gene mutation tests (FIND IT gene panel) performed at Sonic Healthcare during 2018 and early 2019. A total of 1833 lung adenocarcinoma tumour samples underwent somatic mutation testing. EGFR mutations were found in 28% (n=514) of patients, in whom 9.3% (n=48) T790M mutations were present (always combined with other EGFR mutations) and 4.8% (n=25) exon 20 insertions were found. We also compared the prevalence of EGFR mutations identified in our population with that of the four largest publicly available lung cancer cohorts (total n=576 samples). Finally, a subset of 38 samples of primary/and or metastatic lung adenocarcinomas from 23 patients, including five with serial biopsies, underwent detailed morphological analysis. No reproducible morphological correlates were found to be associated with T790M, exon 20 resistance mutations or rarer co-occurring EGFR mutations. Although this may be subject to referral bias towards patients with resistant disease, the incidence of EGFR and T790M mutations is higher in this series from an Australasian population than in other similar publicly available lung adenocarcinoma cohorts. We conclude that histopathological features cannot be used to predict the acquisition of EGFR resistance.

The Proteogenomic Landscape of Curable Prostate Cancer.

DNA sequencing has identified recurrent mutations that drive the aggressiveness of prostate cancers. Surprisingly, the influence of genomic, epigenomic, and transcriptomic dysregulation on the tumor proteome remains poorly understood. We profiled the genomes, epigenomes, transcriptomes, and proteomes of 76 localized, intermediate-risk prostate cancers. We discovered that the genomic subtypes of prostate cancer converge on five proteomic subtypes, with distinct clinical trajectories. ETS fusions, the most common alteration in prostate tumors, affect different genes and pathways in the proteome and transcriptome. Globally, mRNA abundance changes explain only ∼10% of protein abundance variability. As a result, prognostic biomarkers combining genomic or epigenomic features with proteomic ones significantly outperform biomarkers comprised of a single data type.

Amplification of PVT1 contributes to the pathophysiology of ovarian and breast cancer.

PURPOSE: This study was designed to elucidate the role of amplification at 8q24 in the pathophysiology of ovarian and breast cancer because increased copy number at this locus is one of the most frequent genomic abnormalities in these cancers. EXPERIMENTAL DESIGN: To accomplish this, we assessed the association of amplification at 8q24 with outcome in ovarian cancers using fluorescence in situ hybridization to tissue microarrays and measured responses of ovarian and breast cancer cell lines to specific small interfering RNAs against the oncogene MYC and a putative noncoding RNA, PVT1, both of which map to 8q24. RESULTS: Amplification of 8q24 was associated with significantly reduced survival duration. In addition, small interfering RNA-mediated reduction in either PVT1 or MYC expression inhibited proliferation in breast and ovarian cancer cell lines in which they were both amplified and overexpressed but not in lines in which they were not amplified/overexpressed. Inhibition of PVT1 expression also induced a strong apoptotic response in cell lines in which it was overexpressed but not in lines in which it was not amplified/overexpressed. Inhibition of MYC, on the other hand, did not induce an apoptotic response in cell lines in which MYC was amplified and overexpressed. CONCLUSIONS: These results suggest that MYC and PVT1 contribute independently to ovarian and breast pathogenesis when overexpressed because of genomic abnormalities. They also suggest that PVT1-mediated inhibition of apoptosis may explain why amplification of 8q24 is associated with reduced survival duration in patients treated with agents that act through apoptotic mechanisms.

Formalin fixation increases deamination mutation signature but should not lead to false positive mutations in clinical practice.

Genomic analysis of cancer tissues is an essential aspect of personalized oncology treatment. Though it has been suggested that formalin fixation of patient tissues may be suboptimal for molecular studies, this tissue processing approach remains the industry standard. Therefore clinical molecular laboratories must be able to work with formalin fixed, paraffin embedded (FFPE) material. This study examines the effects of pre-analytic variables introduced by routine pathology processing on specimens used for clinical reports produced by next-generation sequencing technology. Tissue resected from three colorectal cancer patients was subjected to 2, 15, 24, and 48 hour fixation times in neutral buffered formalin. DNA was extracted from all tissues twice, once with uracil-N-glycosylase (UNG) treatment to counter deamination effects, and once without. Of note, deamination events at methylated cytosine, as found at CpG sites, remains unaffected by UNG. After extraction a two-step PCR targeted sequencing method was performed using the Illumina MiSeq and the data was analyzed via a custom-built bioinformatics pipeline, including filtration of reads with mapping quality <30. A larger baseline group of samples (n = 20) was examined to establish if there was a sample performance difference between the two DNA extraction methods, with/without UNG treatment. There was no statistical difference between sequencing performance of the two extraction methods when comparing read counts (raw, mapped, and filtered) and read quality (% mapped, % filtered). Analyzing mutation type, there was no significant difference between mutation calls until the 48 hour fixation treatment. At 48 hours there is a significant increase in C/G->T/A mutations that is not represented in DNA treated with UNG. This suggests these errors may be due to deamination events triggered by a longer fixation time. However the allelic frequency of these events remained below the limit of detection for reportable mutations in this assay (<2%). We do however recommend that suspected intratumoral heterogeneity events be verified by re-sequencing the same FFPE block.

The role of mRNA splicing in prostate cancer.

Alternative splicing (AS) is a crucial step in gene expression. It is subject to intricate regulation, and its deregulation in cancer can lead to a wide array of neoplastic phenotypes. A large body of evidence implicates splice isoforms in most if not all hallmarks of cancer, including growth, apoptosis, invasion and metastasis, angiogenesis, and metabolism. AS has important clinical implications since it can be manipulated therapeutically to treat cancer and represents a mechanism of resistance to therapy. In prostate cancer (PCa) AS also plays a prominent role and this review will summarize the current knowledge of alternatively spliced genes with important functional consequences. We will highlight accumulating evidence on AS of the components of the two critical pathways in PCa: androgen receptor (AR) and phosphoinositide 3-kinase (PI3K). These observations together with data on dysregulation of splice factors in PCa suggest that AR and PI3K pathways may be interconnected with previously unappreciated splicing regulatory networks. In addition, we will discuss several lines of evidence implicating splicing regulation in the development of the castration resistance.

A meta-analysis approach for characterizing pan-cancer mechanisms of drug sensitivity in cell lines.

Understanding the heterogeneous drug response of cancer patients is essential to precision oncology. Pioneering genomic analyses of individual cancer subtypes have begun to identify key determinants of resistance, including up-regulation of multi-drug resistance (MDR) genes and mutational alterations of drug targets. However, these alterations are sufficient to explain only a minority of the population, and additional mechanisms of drug resistance or sensitivity are required to explain the remaining spectrum of patient responses to ultimately achieve the goal of precision oncology. We hypothesized that a pan-cancer analysis of in vitro drug sensitivities across numerous cancer lineages will improve the detection of statistical associations and yield more robust and, importantly, recurrent determinants of response. In this study, we developed a statistical framework based on the meta-analysis of expression profiles to identify pan-cancer markers and mechanisms of drug response. Using the Cancer Cell Line Encyclopaedia (CCLE), a large panel of several hundred cancer cell lines from numerous distinct lineages, we characterized both known and novel mechanisms of response to cytotoxic drugs including inhibitors of Topoisomerase 1 (TOP1; Topotecan, Irinotecan) and targeted therapies including inhibitors of histone deacetylases (HDAC; Panobinostat) and MAP/ERK kinases (MEK; PD-0325901, AZD6244). Notably, our analysis implicated reduced replication and transcriptional rates, as well as deficiency in DNA damage repair genes in resistance to TOP1 inhibitors. The constitutive activation of several signaling pathways including the interferon/STAT-1 pathway was implicated in resistance to the pan-HDAC inhibitor. Finally, a number of dysregulations upstream of MEK were identified as compensatory mechanisms of resistance to the MEK inhibitors. In comparison to alternative pan-cancer analysis strategies, our approach can better elucidate relevant drug response mechanisms. Moreover, the compendium of putative markers and mechanisms identified through our analysis can serve as a foundation for future studies into these drugs.

Heterogeneity in the inter-tumor transcriptome of high risk prostate cancer.

BACKGROUND: Genomic analyses of hundreds of prostate tumors have defined a diverse landscape of mutations and genome rearrangements, but the transcriptomic effect of this complexity is less well understood, particularly at the individual tumor level. We selected a cohort of 25 high-risk prostate tumors, representing the lethal phenotype, and applied deep RNA-sequencing and matched whole genome sequencing, followed by detailed molecular characterization. RESULTS: Ten tumors were exposed to neo-adjuvant hormone therapy and expressed marked evidence of therapy response in all except one extreme case, which demonstrated early resistance via apparent neuroendocrine transdifferentiation. We observe high inter-tumor heterogeneity, including unique sets of outlier transcripts in each tumor. Interestingly, outlier expression converged on druggable cellular pathways associated with cell cycle progression, translational control or immune regulation, suggesting distinct contemporary pathway affinity and a mechanism of tumor stratification. We characterize hundreds of novel fusion transcripts, including a high frequency of ETS fusions associated with complex genome rearrangements and the disruption of tumor suppressors. Remarkably, several tumors express unique but potentially-oncogenic non-ETS fusions, which may contribute to the phenotype of individual tumors, and have significance for disease progression. Finally, one ETS-negative tumor has a striking tandem duplication genotype which appears to be highly aggressive and present at low recurrence in ETS-negative prostate cancer, suggestive of a novel molecular subtype. CONCLUSIONS: The multitude of rare genomic and transcriptomic events detected in a high-risk tumor cohort offer novel opportunities for personalized oncology and their convergence on key pathways and functions has broad implications for precision medicine.

Next generation sequencing of prostate cancer from a patient identifies a deficiency of methylthioadenosine phosphorylase, an exploitable tumor target.

Castrate-resistant prostate cancer (CRPC) and neuroendocrine carcinoma of the prostate are invariably fatal diseases for which only palliative therapies exist. As part of a prostate tumor sequencing program, a patient tumor was analyzed using Illumina genome sequencing and a matched renal capsule tumor xenograft was generated. Both tumor and xenograft had a homozygous 9p21 deletion spanning the MTAP, CDKN2, and ARF genes. It is rare for this deletion to occur in primary prostate tumors, yet approximately 10% express decreased levels of methylthioadenosine phosphorylase (MTAP) mRNA. Decreased MTAP expression is a prognosticator for poor outcome. Moreover, it seems that this deletion is more common in CRPC than in primary prostate cancer. We show for the first time that treatment with methylthioadenosine and high dose 6-thioguanine causes marked inhibition of a patient-derived neuroendocrine xenograft growth while protecting the host from 6-thioguanine toxicity. This therapeutic approach can be applied to other MTAP-deficient human cancers as deletion or hypermethylation of the MTAP gene occurs in a broad spectrum of tumors at high frequency. The combination of genome sequencing and patient-derived xenografts can identify candidate therapeutic agents and evaluate them for personalized oncology.

Exon-level microarray analyses identify alternative splicing programs in breast cancer.

Protein isoforms produced by alternative splicing (AS) of many genes have been implicated in several aspects of cancer genesis and progression. These observations motivated a genome-wide assessment of AS in breast cancer. We accomplished this by measuring exon level expression in 31 breast cancer and nonmalignant immortalized cell lines representing luminal, basal, and claudin-low breast cancer subtypes using Affymetrix Human Junction Arrays. We analyzed these data using a computational pipeline specifically designed to detect AS with a low false-positive rate. This identified 181 splice events representing 156 genes as candidates for AS. Reverse transcription-PCR validation of a subset of predicted AS events confirmed 90%. Approximately half of the AS events were associated with basal, luminal, or claudin-low breast cancer subtypes. Exons involved in claudin-low subtype-specific AS were significantly associated with the presence of evolutionarily conserved binding motifs for the tissue-specific Fox2 splicing factor. Small interfering RNA knockdown of Fox2 confirmed the involvement of this splicing factor in subtype-specific AS. The subtype-specific AS detected in this study likely reflects the splicing pattern in the breast cancer progenitor cells in which the tumor arose and suggests the utility of assays for Fox-mediated AS in cancer subtype definition and early detection. These data also suggest the possibility of reducing the toxicity of protein-targeted breast cancer treatments by targeting protein isoforms that are not present in limiting normal tissues.

Organic cation transporters are determinants of oxaliplatin cytotoxicity.

Although the platinum-based anticancer drugs cisplatin, carboplatin, and oxaliplatin have similar DNA-binding properties, only oxaliplatin is active against colorectal tumors. The mechanisms for this tumor specificity of platinum-based compounds are poorly understood but could be related to differences in uptake. This study shows that the human organic cation transporters (OCT) 1 and 2 (SLC22A1 and SLC22A2) markedly increase oxaliplatin, but not cisplatin or carboplatin, accumulation and cytotoxicity in transfected cells, indicating that oxaliplatin is an excellent substrate of these transporters. The cytotoxicity of oxaliplatin was greater than that of cisplatin in six colon cancer cell lines [mean +/- SE of IC(50) in the six cell lines, 3.9 +/- 1.4 micromol/L (oxaliplatin) versus 11 +/- 2.0 micromol/L (cisplatin)] but was reduced by an OCT inhibitor, cimetidine, to a level similar to, or even lower than that of, cisplatin (29 +/- 11 micromol/L for oxaliplatin versus 19 +/- 4.3 micromol/L for cisplatin). Structure-activity studies indicated that organic functionalities on nonleaving groups coordinated to platinum are critical for selective uptake by OCTs. These results indicate that OCT1 and OCT2 are major determinants of the anticancer activity of oxaliplatin and may contribute to its antitumor specificity. They also strongly suggest that expression of OCTs in tumors should be investigated as markers for selecting specific platinum-based therapies in individual patients. The development of new anticancer drugs, specifically targeted to OCTs, represents a novel strategy for targeted drug therapy. The results of the present structure-activity studies indicate specific tactics for realizing this goal.