Cindy Q. Yao

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.

Hypoxia promotes stem cell phenotypes and poor prognosis through epigenetic regulation of DICER

MicroRNAs are small regulatory RNAs that post-transcriptionally control gene expression. Reduced expression of DICER, the enzyme involved in microRNA processing, is frequently observed in cancer and is associated with poor clinical outcome in various malignancies. Yet the underlying mechanisms are not well understood. Here, we identify tumor hypoxia as a regulator of DICER expression in large cohorts of breast cancer patients. We show that DICER expression is suppressed by hypoxia through an epigenetic mechanism that involves inhibition of oxygen-dependent H3K27me3 demethylases KDM6A/B and results in silencing of the DICER promoter. Subsequently, reduced miRNA processing leads to derepression of the miR-200 target ZEB1, stimulates the epithelial to mesenchymal transition and ultimately results in the acquisition of stem cell phenotypes in human mammary epithelial cells. Our study uncovers a previously unknown relationship between oxygen-sensitive epigenetic regulators, miRNA biogenesis and tumor stem cell phenotypes that may underlie poor outcome in breast cancer.

Modulation of long noncoding RNAs by risk SNPs underlying genetic predispositions to prostate cancer

Long noncoding RNAs (lncRNAs) represent an attractive class of candidates to mediate cancer risk. Through integrative analysis of the lncRNA transcriptome with genomic data and SNP data from prostate cancer genome-wide association studies (GWAS), we identified 45 candidate lncRNAs associated with risk to prostate cancer. We further evaluated the mechanism underlying the top hit, PCAT1, and found that a risk-associated variant at rs7463708 increases binding of ONECUT2, a novel androgen receptor (AR)-interacting transcription factor, at a distal enhancer that loops to the PCAT1 promoter, resulting in upregulation of PCAT1 upon prolonged androgen treatment. In addition, PCAT1 interacts with AR and LSD1 and is required for their recruitment to the enhancers of GNMT and DHCR24, two androgen late-response genes implicated in prostate cancer development and progression. PCAT1 promotes prostate cancer cell proliferation and tumor growth in vitro and in vivo. These findings suggest that modulating lncRNA expression is an important mechanism for risk-associated SNPs in promoting prostate transformation.

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 LD₅₀ 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.

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.

Targeted proteomics identifies liquid-biopsy signatures for extracapsular prostate cancer

Biomarkers are rapidly gaining importance in personalized medicine. Although numerous molecular signatures have been developed over the past decade, there is a lack of overlap and many biomarkers fail to validate in independent patient cohorts and hence are not useful for clinical application. For these reasons, identification of novel and robust biomarkers remains a formidable challenge. We combine targeted proteomics with computational biology to discover robust proteomic signatures for prostate cancer. Quantitative proteomics conducted in expressed prostatic secretions from men with extraprostatic and organ-confined prostate cancers identified 133 differentially expressed proteins. Using synthetic peptides, we evaluate them by targeted proteomics in a 74-patient cohort of expressed prostatic secretions in urine. We quantify a panel of 34 candidates in an independent 207-patient cohort. We apply machine-learning approaches to develop clinical predictive models for prostate cancer diagnosis and prognosis. Our results demonstrate that computationally guided proteomics can discover highly accurate non-invasive biomarkers.

Translating a Prognostic DNA Genomic Classifier into the Clinic: Retrospective Validation in 563 Localized Prostate Tumors

BACKGROUND Localized prostate cancer is clinically heterogeneous, despite clinical risk groups that represent relative prostate cancer-specific mortality. We previously developed a 100-locus DNA classifier capable of substratifying patients at risk of biochemical relapse within clinical risk groups. OBJECTIVE The 100-locus genomic classifier was refined to 31 functional loci and tested with standard clinical variables for the ability to predict biochemical recurrence (BCR) and metastasis. DESIGN, SETTING, AND PARTICIPANTS Four retrospective cohorts of radical prostatectomy specimens from patients with localized disease were pooled, and an additional 102-patient cohort used to measure the 31-locus genomic classifier with the NanoString platform. OUTCOME MEASUREMENTS AND STATISTICAL ANALYSIS The genomic classifier scores were tested for their ability to predict BCR (n=563) and metastasis (n=154), and compared with clinical risk stratification schemes. RESULTS AND LIMITATIONS The 31-locus genomic classifier performs similarly to the 100-locus classifier. It identifies patients with elevated BCR rates (hazard ratio=2.73, p<0.001) and patients that eventually develop metastasis (hazard ratio=7.79, p<0.001). Combining the genomic classifier with standard clinical variables outperforms clinical models. Finally, the 31-locus genomic classifier was implemented using a NanoString assay. The study is limited to retrospective cohorts. CONCLUSIONS The 100-locus and 31-locus genomic classifiers reliably identify a cohort of men with localized disease who have an elevated risk of failure. The NanoString assay will be useful for selecting patients for treatment deescalation or escalation in prospective clinical trials based on clinico-genomic scores from pretreatment biopsies. PATIENT SUMMARY It is challenging to determine whether tumors confined to the prostate are aggressive, leading to significant undertreatment and overtreatment. We validated a test based on prostate tumor DNA that improves estimations of relapse risk, and that can help guide treatment planning.

Molecular characterisation of isogenic taxane resistant cell lines identify novel drivers of drug resistance.

BackgroundTaxanes such as paclitaxel and docetaxel are used successfully to treat breast cancer, usually in combination with other agents. They interfere with microtubules causing cell cycle arrest; however, the mechanisms underlying the clinical effects of taxanes are yet to be fully elucidated.MethodsIsogenic paclitaxel resistant (PACR) MDA‒MB‒231, paclitaxel resistant ZR75‒1 and docetaxel resistant (DOCR) ZR75‒1 cell lines were generated by incrementally increasing taxane dose in native cell lines in vitro. We used aCGH analysis to identify mechanisms driving taxane resistance.ResultsTaxane resistant cell lines exhibited an 18-170 fold increased resistance to taxanes, with the ZR75-1 resistant cell lines also demonstrating cross resistance to anthracyclines. Paclitaxel treatment of native cells resulted in a G2/M block and a decrease in the G1 phase of the cell cycle. However, in the resistant cell lines, minimal changes were present. Functional network analysis revealed that the mitotic prometaphase was lost in the resistant cell lines.ConclusionThis study established a model system for examining taxane resistance and demonstrated that both MDR and mitosis represent common mechanism of taxane resistance.

BPG: Seamless, Automated and Interactive Visualization of Scientific Data

We introduce BPG, an easy-to-use framework for generating publication-quality, highly-customizable plots in the R statistical environment. This open-source package includes novel methods of displaying high-dimensional datasets and facilitates generation of complex multi-panel figures, making it ideal for complex datasets. A web-based interactive tool allows online figure customization, from which R code can be downloaded for seamless integration with computational pipelines. BPG is available at http://labs.oicr.on.ca/boutros-lab/software/bpg

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

Background2,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.ResultsSpecifically, 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.ConclusionsAnalysis 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.