Risk SNP-Mediated Promoter-Enhancer Switching Drives Prostate Cancer through lncRNA PCAT19.

The prostate cancer (PCa) risk-associated SNP rs11672691 is positively associated with aggressive disease at diagnosis. We showed that rs11672691 maps to the promoter of a short isoform of long noncoding RNA PCAT19 (PCAT19-short), which is in the third intron of the long isoform (PCAT19-long). The risk variant is associated with decreased and increased levels of PCAT19-short and PCAT19-long, respectively. Mechanistically, the risk SNP region is bifunctional with both promoter and enhancer activity. The risk variants of rs11672691 and its LD SNP rs887391 decrease binding of transcription factors NKX3.1 and YY1 to the promoter of PCAT19-short, resulting in weaker promoter but stronger enhancer activity that subsequently activates PCAT19-long. PCAT19-long interacts with HNRNPAB to activate a subset of cell-cycle genes associated with PCa progression, thereby promoting PCa tumor growth and metastasis. Taken together, these findings reveal a risk SNP-mediated promoter-enhancer switching mechanism underlying both initiation and progression of aggressive PCa.

BPG: Seamless, automated and interactive visualization of scientific data.

BACKGROUND: We introduce BPG, a framework for generating publication-quality, highly-customizable plots in the R statistical environment. RESULTS: This open-source package includes multiple methods of displaying high-dimensional datasets and facilitates generation of complex multi-panel figures, making it suitable for complex datasets. A web-based interactive tool allows online figure customization, from which R code can be downloaded for integration with computational pipelines. CONCLUSION: BPG provides a new approach for linking interactive and scripted data visualization and is available at http://labs.oicr.on.ca/boutros-lab/software/bpg or via CRAN at https://cran.r-project.org/web/packages/BoutrosLab.plotting.general.

NanoStringNormCNV: pre-processing of NanoString CNV data.

Summary: The NanoString System is a well-established technology for measuring RNA and DNA abundance. Although it can estimate copy number variation, relatively few tools support analysis of these data. To address this gap, we created NanoStringNormCNV, an R package for pre-processing and copy number variant calling from NanoString data. This package implements algorithms for pre-processing, quality-control, normalization and copy number variation detection. A series of reporting and data visualization methods support exploratory analyses. To demonstrate its utility, we apply it to a new dataset of 96 genes profiled on 41 prostate tumour and 24 matched normal samples. Availability and implementation: NanoStringNormCNV is implemented in R and is freely available at http://labs.oicr.on.ca/boutros-lab/software/nanostringnormcnv. Contact: paul.boutros@oicr.on.ca. Supplementary information: Supplementary data are available at Bioinformatics online.

Comparative toxicoproteogenomics of mouse and rat liver identifies TCDD-resistance genes.

The aryl hydrocarbon receptor (AHR) mediates many toxic effects of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD). However, the AHR alone does not explain the widely different outcomes among organisms. To identify the other factors involved, we evaluated three transgenic mouse lines, each expressing a different rat AHR isoform (rWT, DEL, and INS) providing widely different resistance to TCDD toxicity, as well as C57BL/6 and DBA/2 mice which exhibit a ~ tenfold divergence in TCDD sensitivity (exposures of 5-1000 µg/kg TCDD). We supplement these with whole-genome sequencing, together with transcriptomic and proteomic analyses of the corresponding rat models, Long-Evans (L-E) and Han/Wistar (H/W) rats (having a ~ 1000-fold difference in their TCDD sensitivities; 100 µg/kg TCDD), to identify genes associated with TCDD-response phenotypes. Overall, we identified up to 50% of genes with altered mRNA abundance following TCDD exposure are associated with a single AHR isoform (33.8%, 11.7%, 5.2% and 0.3% of 3076 genes altered unique to rWT, DEL, C57BL/6 and INS respectively following 1000 µg/kg TCDD). Hepatic Pxdc1 was significantly repressed in all three TCDD-sensitive animal models (C57BL/6 and rWT mice, and L-E rat) after TCDD exposure. Three genes, including Cxxc5, Sugp1 and Hgfac, demonstrated different AHRE-1 (full) motif occurrences within their promoter regions between rat strains, as well as different patterns of mRNA abundance. Several hepatic proteins showed parallel up- or downward alterations with their RNAs, with three genes (SNRK, IGTP and IMPA2) showing consistent, strain-dependent changes. These data show the value of integrating genomic, transcriptomic and proteomic evidence across multi-species models in toxicologic studies.

Compendium of TCDD-mediated transcriptomic response datasets in mammalian model systems.

BACKGROUND: 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) is the most potent congener of the dioxin class of environmental contaminants. Exposure to TCDD causes a wide range of toxic outcomes, ranging from chloracne to acute lethality. The severity of toxicity is highly dependent on the aryl hydrocarbon receptor (AHR). Binding of TCDD to the AHR leads to changes in transcription of numerous genes. Studies evaluating the transcriptional changes brought on by TCDD may provide valuable insight into the role of the AHR in human health and disease. We therefore compiled a collection of transcriptomic datasets that can be used to aid the scientific community in better understanding the transcriptional effects of ligand-activated AHR. RESULTS: Specifically, we have created a datasets package - TCDD.Transcriptomics - for the R statistical environment, consisting of 63 unique experiments comprising 377 samples, including various combinations of 3 species (human derived cell lines, mouse and rat), 4 tissue types (liver, kidney, white adipose tissue and hypothalamus) and a wide range of TCDD exposure times and doses. These datasets have been fully standardized using consistent preprocessing and annotation packages (available as of September 14, 2015). To demonstrate the utility of this R package, a subset of "AHR-core" genes were evaluated across the included datasets. Ahrr, Nqo1 and members of the Cyp family were significantly induced following exposure to TCDD across the studies as expected while Aldh3a1 was induced specifically in rat liver. Inmt was altered only in liver tissue and primarily by rat-AHR. CONCLUSIONS: Analysis of the "AHR-core" genes demonstrates a continued need for studies surrounding the impact of AHR-activity on the transcriptome; genes believed to be consistently regulated by ligand-activated AHR show surprisingly little overlap across species and tissues. Until now, a comprehensive assessment of the transcriptome across these studies was challenging due to differences in array platforms, processing methods and annotation versions. We believe that this package, which is freely available for download ( http://labs.oicr.on.ca/boutros-lab/tcdd-transcriptomics ) will prove to be a highly beneficial resource to the scientific community evaluating the effects of TCDD exposure as well as the variety of functions of the AHR.

Downregulation of histone H2A and H2B pathways is associated with anthracycline sensitivity in breast cancer.

BACKGROUND: Drug resistance in breast cancer is the major obstacle to effective treatment with chemotherapy. While upregulation of multidrug resistance genes is an important component of drug resistance mechanisms in vitro, their clinical relevance remains to be determined. Therefore, identifying pathways that could be targeted in the clinic to eliminate anthracycline-resistant breast cancer remains a major challenge. METHODS: We generated paired native and epirubicin-resistant MDA-MB-231, MCF7, SKBR3 and ZR-75-1 epirubicin-resistant breast cancer cell lines to identify pathways contributing to anthracycline resistance. Native cell lines were exposed to increasing concentrations of epirubicin until resistant cells were generated. To identify mechanisms driving epirubicin resistance, we used a complementary approach including gene expression analyses to identify molecular pathways involved in resistance, and small-molecule inhibitors to reverse resistance. In addition, we tested its clinical relevance in a BR9601 adjuvant clinical trial. RESULTS: Characterisation of epirubicin-resistant cells revealed that they were cross-resistant to doxorubicin and SN-38 and had alterations in apoptosis and cell-cycle profiles. Gene expression analysis identified deregulation of histone H2A and H2B genes in all four cell lines. Histone deacetylase small-molecule inhibitors reversed resistance and were cytotoxic for epirubicin-resistant cell lines, confirming that histone pathways are associated with epirubicin resistance. Gene expression of a novel 18-gene histone pathway module analysis of the BR9601 adjuvant clinical trial revealed that patients with low expression of the 18-gene histone module benefited from anthracycline treatment more than those with high expression (hazard ratio 0.35, 95 % confidence interval 0.13-0.96, p = 0.042). CONCLUSIONS: This study revealed a key pathway that contributes to anthracycline resistance and established model systems for investigating drug resistance in all four major breast cancer subtypes. As the histone modification can be targeted with small-molecule inhibitors, it represents a possible means of reversing clinical anthracycline resistance. TRIAL REGISTRATION: ClinicalTrials.gov identifier NCT00003012 . Registered on 1 November 1999.

Tumour genomic and microenvironmental heterogeneity for integrated prediction of 5-year biochemical recurrence of prostate cancer: a retrospective cohort study.

BACKGROUND: Clinical prognostic groupings for localised prostate cancers are imprecise, with 30-50% of patients recurring after image-guided radiotherapy or radical prostatectomy. We aimed to test combined genomic and microenvironmental indices in prostate cancer to improve risk stratification and complement clinical prognostic factors. METHODS: We used DNA-based indices alone or in combination with intra-prostatic hypoxia measurements to develop four prognostic indices in 126 low-risk to intermediate-risk patients (Toronto cohort) who will receive image-guided radiotherapy. We validated these indices in two independent cohorts of 154 (Memorial Sloan Kettering Cancer Center cohort [MSKCC] cohort) and 117 (Cambridge cohort) radical prostatectomy specimens from low-risk to high-risk patients. We applied unsupervised and supervised machine learning techniques to the copy-number profiles of 126 pre-image-guided radiotherapy diagnostic biopsies to develop prognostic signatures. Our primary endpoint was the development of a set of prognostic measures capable of stratifying patients for risk of biochemical relapse 5 years after primary treatment. FINDINGS: Biochemical relapse was associated with indices of tumour hypoxia, genomic instability, and genomic subtypes based on multivariate analyses. We identified four genomic subtypes for prostate cancer, which had different 5-year biochemical relapse-free survival. Genomic instability is prognostic for relapse in both image-guided radiotherapy (multivariate analysis hazard ratio [HR] 4·5 [95% CI 2·1-9·8]; p=0·00013; area under the receiver operator curve [AUC] 0·70 [95% CI 0·65-0·76]) and radical prostatectomy (4·0 [1·6-9·7]; p=0·0024; AUC 0·57 [0·52-0·61]) patients with prostate cancer, and its effect is magnified by intratumoral hypoxia (3·8 [1·2-12]; p=0·019; AUC 0·67 [0·61-0·73]). A novel 100-loci DNA signature accurately classified treatment outcome in the MSKCC low-risk to intermediate-risk cohort (multivariate analysis HR 6·1 [95% CI 2·0-19]; p=0·0015; AUC 0·74 [95% CI 0·65-0·83]). In the independent MSKCC and Cambridge cohorts, this signature identified low-risk to high-risk patients who were most likely to fail treatment within 18 months (combined cohorts multivariate analysis HR 2·9 [95% CI 1·4-6·0]; p=0·0039; AUC 0·68 [95% CI 0·63-0·73]), and was better at predicting biochemical relapse than 23 previously published RNA signatures. INTERPRETATION: This is the first study of cancer outcome to integrate DNA-based and microenvironment-based failure indices to predict patient outcome. Patients exhibiting these aggressive features after biopsy should be entered into treatment intensification trials. FUNDING: Movember Foundation, Prostate Cancer Canada, Ontario Institute for Cancer Research, Canadian Institute for Health Research, NIHR Cambridge Biomedical Research Centre, The University of Cambridge, Cancer Research UK, Cambridge Cancer Charity, Prostate Cancer UK, Hutchison Whampoa Limited, Terry Fox Research Institute, Princess Margaret Cancer Centre Foundation, PMH-Radiation Medicine Program Academic Enrichment Fund, Motorcycle Ride for Dad (Durham), Canadian Cancer Society.

Identification of differentially expressed proteins in direct expressed prostatic secretions of men with organ-confined versus extracapsular prostate cancer.

Current protocols for the screening of prostate cancer cannot accurately discriminate clinically indolent tumors from more aggressive ones. One reliable indicator of outcome has been the determination of organ-confined versus nonorgan-confined disease but even this determination is often only made following prostatectomy. This underscores the need to explore alternate avenues to enhance outcome prediction of prostate cancer patients. Fluids that are proximal to the prostate, such as expressed prostatic secretions (EPS), are attractive sources of potential prostate cancer biomarkers as these fluids likely bathe the tumor. Direct-EPS samples from 16 individuals with extracapsular (n = 8) or organ-confined (n = 8) prostate cancer were used as a discovery cohort, and were analyzed in duplicate by a nine-step MudPIT on a LTQ-Orbitrap XL mass spectrometer. A total of 624 unique proteins were identified by at least two unique peptides with a 0.2% false discovery rate. A semiquantitative spectral counting algorithm identified 133 significantly differentially expressed proteins in the discovery cohort. Integrative data mining prioritized 14 candidates, including two known prostate cancer biomarkers: prostate-specific antigen and prostatic acid phosphatase, which were significantly elevated in the direct-EPS from the organ-confined cancer group. These and five other candidates (SFN, MME, PARK7, TIMP1, and TGM4) were verified by Western blotting in an independent set of direct-EPS from patients with biochemically recurrent disease (n = 5) versus patients with no evidence of recurrence upon follow-up (n = 10). Lastly, we performed proof-of-concept SRM-MS-based relative quantification of the five candidates using unpurified heavy isotope-labeled synthetic peptides spiked into pools of EPS-urines from men with extracapsular and organ-confined prostate tumors. This study represents the first efforts to define the direct-EPS proteome from two major subclasses of prostate cancer using shotgun proteomics and verification in EPS-urine by SRM-MS.

Inter-strain heterogeneity in rat hepatic transcriptomic responses to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD).

The biochemical and toxic effects of 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) have been the subject of intense study for decades. It is now clear that essentially all TCDD-induced toxicities are mediated by DNA-protein interactions involving the Aryl Hydrocarbon Receptor (AHR). Nevertheless, it remains unknown which AHR target genes cause TCDD toxicities. Several groups, including our own, have developed rodent model systems to probe these questions. mRNA expression profiling of these model systems has revealed significant inter-species heterogeneity in rodent hepatic responses to TCDD. It has remained unclear if this variability also exists within a species, amongst rodent strains. To resolve this question, we profiled the hepatic transcriptomic response to TCDD of diverse rat strains (L-E, H/W, F344 and Wistar rats) and two lines derived from L-E×H/W crosses, at consistent age, sex, and dosing (100 µg/kg TCDD for 19 h). Using this uniquely consistent dataset, we show that the majority of TCDD-induced alterations in mRNA abundance are strain/line-specific: only 11 genes were affected by TCDD across all strains, including well-known dioxin-responsive genes such as Cyp1a1 and Nqo1. Our analysis identified two novel universally dioxin-responsive genes as well as 4 genes induced by TCDD in dioxin-sensitive rats only. These 6 genes are strong candidates to explain TCDD-related toxicities, so we validated them using 152 animals in time-course (0 to 384 h) and dose-response (0 to 3000 µg/kg) experiments. This study reveals that different rat strains exhibit dramatic transcriptional heterogeneity in their hepatic responses to TCDD and that inter-strain comparisons can help identify candidate toxicity-related genes.

Hepatic transcriptomic responses to TCDD in dioxin-sensitive and dioxin-resistant rats during the onset of toxicity.

The dioxin congener 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) causes a wide range of toxic effects in rodent species, all of which are mediated by a ligand-dependent transcription-factor, the aryl hydrocarbon receptor (AHR). The Han/Wistar (Kuopio) (H/W) strain shows exceptional resistance to many TCDD-induced toxicities; the LD50 of > 9600 µg/kg for H/W rats is higher than for any other wild-type mammal known. We previously showed that this resistance primarily results from H/W rats expressing a variant AHR isoform that has a substantial portion of the AHR transactivation domain deleted. Despite this large deletion, H/W rats are not entirely refractory to the effects of TCDD; the variant AHR in these animals remains fully competent to up-regulate well-known dioxin-inducible genes. TCDD-sensitive (Long-Evans, L-E) and resistant (H/W) rats were treated with either corn-oil (with or without feed-restriction) or 100 µg/kg TCDD for either four or ten days. Hepatic transcriptional profiling was done using microarrays, and was validated by RT-PCR analysis of 41 genes. A core set of genes was altered in both strains at all time points tested, including CYP1A1, CYP1A2, CYP1B1, Nqo1, Aldh3a1, Tiparp, Exoc3, and Inmt. Outside this core, the strains differed significantly in the breadth of response: three-fold more genes were altered in L-E than H/W rats. At ten days almost all expressed genes were dysregulated in L-E rats, likely reflecting emerging toxic responses. Far fewer genes were affected by feed-restriction, suggesting that only a minority of the TCDD-induced changes are secondary to the wasting syndrome.

Comparative survival analysis of multiparametric tests-when molecular tests disagree-A TEAM Pathology study.

Multiparametric assays for risk stratification are widely used in the management of both node negative and node positive hormone receptor positive invasive breast cancer. Recent data from multiple sources suggests that different tests may provide different risk estimates at the individual patient level. The TEAM pathology study consists of 3284 postmenopausal ER+ve breast cancers treated with endocrine therapy Using genes comprising the following multi-parametric tests OncotypeDx®, Prosigna™ and MammaPrint® signatures were trained to recapitulate true assay results. Patients were then classified into risk groups and survival assessed. Whilst likelihood χ2 ratios suggested limited value for combining tests, Kaplan-Meier and LogRank tests within risk groups suggested combinations of tests provided statistically significant stratification of potential clinical value. Paradoxically whilst Prosigna-trained results stratified Oncotype-trained subgroups across low and intermediate risk categories, only intermediate risk Prosigna-trained cases were further stratified by Oncotype-trained results. Both Oncotype-trained and Prosigna-trained results further stratified MammaPrint-trained low risk cases, and MammaPrint-trained results also stratified Oncotype-trained low and intermediate risk groups but not Prosigna-trained results. Comparisons between existing multiparametric tests are challenging, and evidence on discordance between tests in risk stratification presents further dilemmas. Detailed analysis of the TEAM pathology study suggests a complex inter-relationship between test results in the same patient cohorts which requires careful evaluation regarding test utility. Further prognostic improvement appears both desirable and achievable.

Computational approaches to support comparative analysis of multiparametric tests: Modelling versus Training.

Multiparametric assays for risk stratification are widely used in the management of breast cancer, with applications being developed for a number of other cancer settings. Recent data from multiple sources suggests that different tests may provide different risk estimates at the individual patient level. There is an increasing need for robust methods to support cost effective comparisons of test performance in multiple settings. The derivation of similar risk classifications using genes comprising the following multi-parametric tests Oncotype DX® (Genomic Health.), Prosigna™ (NanoString Technologies, Inc.), MammaPrint® (Agendia Inc.) was performed using different computational approaches. Results were compared to the actual test results. Two widely used approaches were applied, firstly computational "modelling" of test results using published algorithms and secondly a "training" approach which used reference results from the commercially supplied tests. We demonstrate the potential for errors to arise when using a "modelling" approach without reference to real world test results. Simultaneously we show that a "training" approach can provide a highly cost-effective solution to the development of real-world comparisons between different multigene signatures. Comparisons between existing multiparametric tests is challenging, and evidence on discordance between tests in risk stratification presents further dilemmas. We present an approach, modelled in breast cancer, which can provide health care providers and researchers with the potential to perform robust and meaningful comparisons between multigene tests in a cost-effective manner. We demonstrate that whilst viable estimates of gene signatures can be derived from modelling approaches, in our study using a training approach allowed a close approximation to true signature results.

Temporal Stability and Prognostic Biomarker Potential of the Prostate Cancer Urine miRNA Transcriptome.

BACKGROUND: The development of noninvasive tests for the early detection of aggressive prostate tumors is a major unmet clinical need. miRNAs are promising noninvasive biomarkers: they play essential roles in tumorigenesis, are stable under diverse analytical conditions, and can be detected in body fluids. METHODS: We measured the longitudinal stability of 673 miRNAs by collecting serial urine samples from 10 patients with localized prostate cancer. We then measured temporally stable miRNAs in an independent training cohort (n = 99) and created a biomarker predictive of Gleason grade using machine-learning techniques. Finally, we validated this biomarker in an independent validation cohort (n = 40). RESULTS: We found that each individual has a specific urine miRNA fingerprint. These fingerprints are temporally stable and associated with specific biological functions. We identified seven miRNAs that were stable over time within individual patients and integrated them with machine-learning techniques to create a novel biomarker for prostate cancer that overcomes interindividual variability. Our urine biomarker robustly identified high-risk patients and achieved similar accuracy as tissue-based prognostic markers (area under the receiver operating characteristic = 0.72, 95% confidence interval = 0.69 to 0.76 in the training cohort, and area under the receiver operating characteristic curve = 0.74, 95% confidence interval = 0.55 to 0.92 in the validation cohort). CONCLUSIONS: These data highlight the importance of quantifying intra- and intertumoral heterogeneity in biomarker development. This noninvasive biomarker may usefully supplement invasive or expensive radiologic- and tissue-based assays.

2,3,7,8-Tetrachlorodibenzo-p-dioxin modifies alternative splicing in mouse liver.

Alternative splicing is a co-transcriptional mechanism that generates protein diversity by including or excluding exons in different combinations, thereby expanding the diversity of protein isoforms of a single gene. Abnormalities in this process can result in deleterious effects to human health, and several xenobiotics are known to interfere with splicing regulation through multiple mechanisms. These changes could lead to human diseases such as cancer, neurological disorders, autoimmune diseases, and developmental disorders. 2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) is an environmental contaminant generated as a byproduct of various industrial activities. Exposure to this dioxin has been linked to a wide range of pathologies through the alteration of multiple cellular processes. However, the effects of TCDD exposure on alternative splicing have not yet been studied. Here, we investigated whether a single po. dose of 5 µg/kg or 500 µg/kg TCDD influence hepatic alternative splicing in adult male C57BL/6Kou mouse. We identified several genes whose alternative splicing of precursor messenger RNAs was modified following TCDD exposure. In particular, we demonstrated that alternative splicing of Cyp1a1, Ahrr, and Actn1 was significantly altered after TCDD treatment. These findings show that the exposure to TCDD has an impact on alternative-splicing, and suggest a new avenue for understanding TCDD-mediated toxicity and pathogenesis.

Identification of Distinct Prognostic Groups: Implications for Patient Selection to Targeted Therapies Among Anti-Endocrine Therapy-Resistant Early Breast Cancers.

PURPOSE: Hormone receptor-positive breast cancer remains an ongoing therapeutic challenge, despite optimal anti-endocrine therapies. In this study, we assessed the prognostic ability of genomic signatures to identify patients at risk for recurrence after endocrine therapy. Analysis was performed on the basis of an a priori hypothesis related to molecular pathways, which might predict response to existing targeted therapies. PATIENTS AND METHODS: A subset of patients from the Tamoxifen Versus Exemestane Adjuvant Multinational trial (ClinicalTrials.gov identifiers: NCT00279448 and NCT00032136, and NCT00036270) pathology cohort were analyzed to determine the prognostic ability of mutational and copy number aberration biomarkers that represent the cyclin D/cyclin-dependent kinase (CCND/CDK), fibroblast growth factor receptor/fibroblast growth factor (FGFR/FGF), and phosphatidylinositol 3-kinase/protein kinase B (PI3K/ATK) pathways to inform the potential choice of additional therapies to standard endocrine treatment. Copy number analysis and targeted sequencing was performed. Pathways were identified as aberrant if there were copy number aberrations and/or mutations in any of the predetermined pathway genes: CCND1/CCND2/CCND3/CDK4/CDK6, FGFR1/FGFR2/FGFR2/FGFR4, and AKT1/AKT2/PIK3CA/PTEN. RESULTS: The 390 of 420 samples that passed quality control were analyzed for distant metastasis-free survival between groups. Patients with no changes in the CCND/CDK pathway experienced a better distant metastasis-free survival (hazard ratio, 1.94; 95% CI, 1.45 to 2.61; P < .001) than those who possessed aberrations. In the FGFR/FGF and PI3K/AKT pathways, a similar outcome was observed (hazard ratio, 1.43 [95% CI, 1.07 to 1.92; P = .017] and 1.34 [95% CI, 1.00 to 1.81; P = .053], respectively). CONCLUSION: We show that aberrations of genes in these pathways are independently linked to a higher risk of relapse after endocrine treatment. Improvement of the clinical management of early breast cancers could be made by identifying those for whom current endocrine therapies are sufficient, thus reducing unnecessary treatment, and secondly, by identifying those who are at high risk for recurrence and linking molecular features that drive these cancers to treatment with targeted therapies.

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.

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.

Genome-wide germline correlates of the epigenetic landscape of prostate cancer.

Oncogenesis is driven by germline, environmental and stochastic factors. It is unknown how these interact to produce the molecular phenotypes of tumors. We therefore quantified the influence of germline polymorphisms on the somatic epigenome of 589 localized prostate tumors. Predisposition risk loci influence a tumor's epigenome, uncovering a mechanism for cancer susceptibility. We identified and validated 1,178 loci associated with altered methylation in tumoral but not nonmalignant tissue. These tumor methylation quantitative trait loci influence chromatin structure, as well as RNA and protein abundance. One prominent tumor methylation quantitative trait locus is associated with AKT1 expression and is predictive of relapse after definitive local therapy in both discovery and validation cohorts. These data reveal intricate crosstalk between the germ line and the epigenome of primary tumors, which may help identify germline biomarkers of aggressive disease to aid patient triage and optimize the use of more invasive or expensive diagnostic assays.

Deregulation of the spindle assembly checkpoint is associated with paclitaxel resistance in ovarian cancer.

BACKGROUND: Ovarian cancer is the leading gynecologic cancer diagnosed in North America and because related symptoms are not disease specific, this often leads to late detection, an advanced disease state, and the need for chemotherapy. Ovarian cancer is frequently sensitive to chemotherapy at diagnosis but rapid development of drug resistance leads to disease progression and ultimately death in the majority of patients. RESULTS: We have generated paclitaxel resistant ovarian cell lines from their corresponding native cell lines to determine driver mechanisms of drug resistance using gene expression arrays. These paclitaxel resistant ovarian cells demonstrate: (1) Increased IC50 for paclitaxel and docetaxel (10 to 75-fold) and cross-resistance to anthracyclines (2) Reduced cell apoptosis in the presence of paclitaxel (3) Gene depletion involving mitotic regulators BUB1 mitotic checkpoint serine/threonine kinase, cyclin BI (CCNB1), centromere protein E (CENPE), and centromere protein F (CENPF), and (4) Functional data validating gene depletion among mitotic regulators. CONCLUSIONS: We have generated model systems to explore drug resistance in ovarian cancer, which have revealed a key pathway related to the spindle assembly checkpoint underlying paclitaxel resistance in ovarian cell lines.

Pathway-based subnetworks enable cross-disease biomarker discovery.

Biomarkers lie at the heart of precision medicine. Surprisingly, while rapid genomic profiling is becoming ubiquitous, the development of biomarkers usually involves the application of bespoke techniques that cannot be directly applied to other datasets. There is an urgent need for a systematic methodology to create biologically-interpretable molecular models that robustly predict key phenotypes. Here we present SIMMS (Subnetwork Integration for Multi-Modal Signatures): an algorithm that fragments pathways into functional modules and uses these to predict phenotypes. We apply SIMMS to multiple data types across five diseases, and in each it reproducibly identifies known and novel subtypes, and makes superior predictions to the best bespoke approaches. To demonstrate its ability on a new dataset, we profile 33 genes/nodes of the PI3K pathway in 1734 FFPE breast tumors and create a four-subnetwork prediction model. This model out-performs a clinically-validated molecular test in an independent cohort of 1742 patients. SIMMS is generic and enables systematic data integration for robust biomarker discovery.