Accurate detection of benign and malignant renal tumor subtypes with MethylBoostER: An epigenetic marker-driven learning framework.

Current gold standard diagnostic strategies are unable to accurately differentiate malignant from benign small renal masses preoperatively; consequently, 20% of patients undergo unnecessary surgery. Devising a more confident presurgical diagnosis is key to improving treatment decision-making. We therefore developed MethylBoostER, a machine learning model leveraging DNA methylation data from 1228 tissue samples, to classify pathological subtypes of renal tumors (benign oncocytoma, clear cell, papillary, and chromophobe RCC) and normal kidney. The prediction accuracy in the testing set was 0.960, with class-wise ROC AUCs >0.988 for all classes. External validation was performed on >500 samples from four independent datasets, achieving AUCs >0.89 for all classes and average accuracies of 0.824, 0.703, 0.875, and 0.894 for the four datasets. Furthermore, consistent classification of multiregion samples (N = 185) from the same patient demonstrates that methylation heterogeneity does not limit model applicability. Following further clinical studies, MethylBoostER could facilitate a more confident presurgical diagnosis to guide treatment decision-making in the future.

Comprehensive characterization of cell-free tumor DNA in plasma and urine of patients with renal tumors.

BACKGROUND: Cell-free tumor-derived DNA (ctDNA) allows non-invasive monitoring of cancers, but its utility in renal cell cancer (RCC) has not been established. METHODS: Here, a combination of untargeted and targeted sequencing methods, applied to two independent cohorts of patients (n = 91) with various renal tumor subtypes, were used to determine ctDNA content in plasma and urine. RESULTS: Our data revealed lower plasma ctDNA levels in RCC relative to other cancers of similar size and stage, with untargeted detection in 27.5% of patients from both cohorts. A sensitive personalized approach, applied to plasma and urine from select patients (n = 22) improved detection to ~ 50%, including in patients with early-stage disease and even benign lesions. Detection in plasma, but not urine, was more frequent amongst patients with larger tumors and in those patients with venous tumor thrombus. With data from one extensively characterized patient, we observed that plasma and, for the first time, urine ctDNA may better represent tumor heterogeneity than a single tissue biopsy. Furthermore, in a subset of patients (n = 16), longitudinal sampling revealed that ctDNA can track disease course and may pre-empt radiological identification of minimal residual disease or disease progression on systemic therapy. Additional datasets will be required to validate these findings. CONCLUSIONS: These data highlight RCC as a ctDNA-low malignancy. The biological reasons for this are yet to be determined. Nonetheless, our findings indicate potential clinical utility in the management of patients with renal tumors, provided improvement in isolation and detection approaches.

Measurement of Plasma Cell-Free Mitochondrial Tumor DNA Improves Detection of Glioblastoma in Patient-Derived Orthotopic Xenograft Models.

The factors responsible for the low detection rate of cell-free tumor DNA (ctDNA) in the plasma of patients with glioblastoma (GBM) are currently unknown. In this study, we measured circulating nucleic acids in patient-derived orthotopically implanted xenograft (PDOX) models of GBM (n = 64) and show that tumor size and cell proliferation, but not the integrity of the blood-brain barrier or cell death, affect the release of ctDNA in treatment-naïve GBM PDOX. Analysis of fragment length profiles by shallow genome-wide sequencing (<0.2× coverage) of host (rat) and tumor (human) circulating DNA identified a peak at 145 bp in the human DNA fragments, indicating a difference in the origin or processing of the ctDNA. The concentration of ctDNA correlated with cell death only after treatment with temozolomide and radiotherapy. Digital PCR detection of plasma tumor mitochondrial DNA (tmtDNA), an alternative to detection of nuclear ctDNA, improved plasma DNA detection rate (82% vs. 24%) and allowed detection in cerebrospinal fluid and urine. Mitochondrial mutations are prevalent across all cancers and can be detected with high sensitivity, at low cost, and without prior knowledge of tumor mutations via capture-panel sequencing. Coupled with the observation that mitochondrial copy number increases in glioma, these data suggest analyzing tmtDNA as a more sensitive method to detect and monitor tumor burden in cancer, specifically in GBM, where current methods have largely failed. SIGNIFICANCE: These findings show that detection of tumor mitochondrial DNA is more sensitive than circulating tumor DNA analysis to detect and monitor tumor burden in patient-derived orthotopic xenografts of glioblastoma.

Enhanced detection of circulating tumor DNA by fragment size analysis.

Existing methods to improve detection of circulating tumor DNA (ctDNA) have focused on genomic alterations but have rarely considered the biological properties of plasma cell-free DNA (cfDNA). We hypothesized that differences in fragment lengths of circulating DNA could be exploited to enhance sensitivity for detecting the presence of ctDNA and for noninvasive genomic analysis of cancer. We surveyed ctDNA fragment sizes in 344 plasma samples from 200 patients with cancer using low-pass whole-genome sequencing (0.4×). To establish the size distribution of mutant ctDNA, tumor-guided personalized deep sequencing was performed in 19 patients. We detected enrichment of ctDNA in fragment sizes between 90 and 150 bp and developed methods for in vitro and in silico size selection of these fragments. Selecting fragments between 90 and 150 bp improved detection of tumor DNA, with more than twofold median enrichment in >95% of cases and more than fourfold enrichment in >10% of cases. Analysis of size-selected cfDNA identified clinically actionable mutations and copy number alterations that were otherwise not detected. Identification of plasma samples from patients with advanced cancer was improved by predictive models integrating fragment length and copy number analysis of cfDNA, with area under the curve (AUC) >0.99 compared to AUC <0.80 without fragmentation features. Increased identification of cfDNA from patients with glioma, renal, and pancreatic cancer was achieved with AUC > 0.91 compared to AUC < 0.5 without fragmentation features. Fragment size analysis and selective sequencing of specific fragment sizes can boost ctDNA detection and could complement or provide an alternative to deeper sequencing of cfDNA.

Synthetic lethality between androgen receptor signalling and the PARP pathway in prostate cancer.

Emerging data demonstrate homologous recombination (HR) defects in castration-resistant prostate cancers, rendering these tumours sensitive to PARP inhibition. Here we demonstrate a direct requirement for the androgen receptor (AR) to maintain HR gene expression and HR activity in prostate cancer. We show that PARP-mediated repair pathways are upregulated in prostate cancer following androgen-deprivation therapy (ADT). Furthermore, upregulation of PARP activity is essential for the survival of prostate cancer cells and we demonstrate a synthetic lethality between ADT and PARP inhibition in vivo. Our data suggest that ADT can functionally impair HR prior to the development of castration resistance and that, this potentially could be exploited therapeutically using PARP inhibitors in combination with androgen-deprivation therapy upfront in advanced or high-risk prostate cancer.Tumours with homologous recombination (HR) defects become sensitive to PARPi. Here, the authors show that androgen receptor (AR) regulates HR and AR inhibition activates the PARP pathway in vivo, thus inhibition of both AR and PARP is required for effective treatment of high risk prostate cancer.

The importance of DNA methylation in prostate cancer development.

After briefly reviewing the nature of DNA methylation, its general role in cancer and the tools available to interrogate it, we consider the literature surrounding DNA methylation as relating to prostate cancer. Specific consideration is given to recurrent alterations. A list of frequently reported genes is synthesized from 17 studies that have reported on methylation changes in malignant prostate tissue, and we chart the timing of those changes in the diseases history through amalgamation of several previously published data sets. We also review associations with genetic alterations and hormone signalling, before the practicalities of investigating prostate cancer methylation using cell lines are assessed. We conclude by outlining the interplay between DNA methylation and prostate cancer metabolism and their regulation by androgen receptor, with a specific discussion of the mitochondria and their associations with DNA methylation.

HNF1B variants associate with promoter methylation and regulate gene networks activated in prostate and ovarian cancer.

Two independent regions within HNF1B are consistently identified in prostate and ovarian cancer genome-wide association studies (GWAS); their functional roles are unclear. We link prostate cancer (PC) risk SNPs rs11649743 and rs3760511 with elevated HNF1B gene expression and allele-specific epigenetic silencing, and outline a mechanism by which common risk variants could effect functional changes that increase disease risk: functional assays suggest that HNF1B is a pro-differentiation factor that suppresses epithelial-to-mesenchymal transition (EMT) in unmethylated, healthy tissues. This tumor-suppressor activity is lost when HNF1B is silenced by promoter methylation in the progression to PC. Epigenetic inactivation of HNF1B in ovarian cancer also associates with known risk SNPs, with a similar impact on EMT. This represents one of the first comprehensive studies into the pleiotropic role of a GWAS-associated transcription factor across distinct cancer types, and is the first to describe a conserved role for a multi-cancer genetic risk factor.

Methods to Identify Chromatin-Bound Protein Complexes: From Genome-Wide to Locus-Specific Approaches.

High-throughput sequencing approaches coupled with functional genomics experiments have facilitated a rapid growth in our understanding of chromatin biology, from genome-wide maps of transcription factor binding and histone modifications to insights into higher order chromatin organization under specific cellular conditions. However in most cases these methods require a prior knowledge of the system of interest (e.g., targets for immunoprecipitation or modulation) and therefore are limited in their utility to identify novel components of pathways or for the study of uncharacterized pathways. Several orthologous proteomics approaches have been developed recently that bridge this gap, allowing the identification of protein complexes globally or at specific genomic loci. In this chapter the relative advantages of each approach will be explored and a detailed protocol given for DNA pull-down of a specific androgen receptor (AR) genomic target.

Choline Kinase Alpha as an Androgen Receptor Chaperone and Prostate Cancer Therapeutic Target.

BACKGROUND: The androgen receptor (AR) is a major drug target in prostate cancer (PCa). We profiled the AR-regulated kinome to identify clinically relevant and druggable effectors of AR signaling. METHODS: Using genome-wide approaches, we interrogated all AR regulated kinases. Among these, choline kinase alpha (CHKA) expression was evaluated in benign (n = 195), prostatic intraepithelial neoplasia (PIN) (n = 153) and prostate cancer (PCa) lesions (n = 359). We interrogated how CHKA regulates AR signaling using biochemical assays and investigated androgen regulation of CHKA expression in men with PCa, both untreated (n = 20) and treated with an androgen biosynthesis inhibitor degarelix (n = 27). We studied the effect of CHKA inhibition on the PCa transcriptome using RNA sequencing and tested the effect of CHKA inhibition on cell growth, clonogenic survival and invasion. Tumor xenografts (n = 6 per group) were generated in mice using genetically engineered prostate cancer cells with inducible CHKA knockdown. Data were analyzed with χ(2) tests, Cox regression analysis, and Kaplan-Meier methods. All statistical tests were two-sided. RESULTS: CHKA expression was shown to be androgen regulated in cell lines, xenografts, and human tissue (log fold change from 6.75 to 6.59, P = .002) and was positively associated with tumor stage. CHKA binds directly to the ligand-binding domain (LBD) of AR, enhancing its stability. As such, CHKA is the first kinase identified as an AR chaperone. Inhibition of CHKA repressed the AR transcriptional program including pathways enriched for regulation of protein folding, decreased AR protein levels, and inhibited the growth of PCa cell lines, human PCa explants, and tumor xenografts. CONCLUSIONS: CHKA can act as an AR chaperone, providing, to our knowledge, the first evidence for kinases as molecular chaperones, making CHKA both a marker of tumor progression and a potential therapeutic target for PCa.

Epigenetic and oncogenic regulation of SLC16A7 (MCT2) results in protein over-expression, impacting on signalling and cellular phenotypes in prostate cancer.

Monocarboxylate Transporter 2 (MCT2) is a major pyruvate transporter encoded by the SLC16A7 gene. Recent studies pointed to a consistent overexpression of MCT2 in prostate cancer (PCa) suggesting MCT2 as a putative biomarker and molecular target. Despite the importance of this observation the mechanisms involved in MCT2 regulation are unknown. Through an integrative analysis we have discovered that selective demethylation of an internal SLC16A7/MCT2 promoter is a recurrent event in independent PCa cohorts. This demethylation is associated with expression of isoforms differing only in 5'-UTR translational control motifs, providing one contributing mechanism for MCT2 protein overexpression in PCa. Genes co-expressed with SLC16A7/MCT2 also clustered in oncogenic-related pathways and effectors of these signalling pathways were found to bind at the SLC16A7/MCT2 gene locus. Finally, MCT2 knock-down attenuated the growth of PCa cells. The present study unveils an unexpected epigenetic regulation of SLC16A7/MCT2 isoforms and identifies a link between SLC16A7/MCT2, Androgen Receptor (AR), ETS-related gene (ERG) and other oncogenic pathways in PCa. These results underscore the importance of combining data from epigenetic, transcriptomic and protein level changes to allow more comprehensive insights into the mechanisms underlying protein expression, that in our case provide additional weight to MCT2 as a candidate biomarker and molecular target in PCa.

Effect of mutation order on myeloproliferative neoplasms.

BACKGROUND: Cancers result from the accumulation of somatic mutations, and their properties are thought to reflect the sum of these mutations. However, little is known about the effect of the order in which mutations are acquired. METHODS: We determined mutation order in patients with myeloproliferative neoplasms by genotyping hematopoietic colonies or by means of next-generation sequencing. Stem cells and progenitor cells were isolated to study the effect of mutation order on mature and immature hematopoietic cells. RESULTS: The age at which a patient presented with a myeloproliferative neoplasm, acquisition of JAK2 V617F homozygosity, and the balance of immature progenitors were all influenced by mutation order. As compared with patients in whom the TET2 mutation was acquired first (hereafter referred to as "TET2-first patients"), patients in whom the Janus kinase 2 (JAK2) mutation was acquired first ("JAK2-first patients") had a greater likelihood of presenting with polycythemia vera than with essential thrombocythemia, an increased risk of thrombosis, and an increased sensitivity of JAK2-mutant progenitors to ruxolitinib in vitro. Mutation order influenced the proliferative response to JAK2 V617F and the capacity of double-mutant hematopoietic cells and progenitor cells to generate colony-forming cells. Moreover, the hematopoietic stem-and-progenitor-cell compartment was dominated by TET2 single-mutant cells in TET2-first patients but by JAK2-TET2 double-mutant cells in JAK2-first patients. Prior mutation of TET2 altered the transcriptional consequences of JAK2 V617F in a cell-intrinsic manner and prevented JAK2 V617F from up-regulating genes associated with proliferation. CONCLUSIONS: The order in which JAK2 and TET2 mutations were acquired influenced clinical features, the response to targeted therapy, the biology of stem and progenitor cells, and clonal evolution in patients with myeloproliferative neoplasms. (Funded by Leukemia and Lymphoma Research and others.).

HES5 silencing is an early and recurrent change in prostate tumourigenesis.

Prostate cancer is the most common cancer in men, resulting in over 10 000 deaths/year in the UK. Sequencing and copy number analysis of primary tumours has revealed heterogeneity within tumours and an absence of recurrent founder mutations, consistent with non-genetic disease initiating events. Using methylation profiling in a series of multi-focal prostate tumours, we identify promoter methylation of the transcription factor HES5 as an early event in prostate tumourigenesis. We confirm that this epigenetic alteration occurs in 86-97% of cases in two independent prostate cancer cohorts (n=49 and n=39 tumour-normal pairs). Treatment of prostate cancer cells with the demethylating agent 5-aza-2'-deoxycytidine increased HES5 expression and downregulated its transcriptional target HES6, consistent with functional silencing of the HES5 gene in prostate cancer. Finally, we identify and test a transcriptional module involving the AR, ERG, HES1 and HES6 and propose a model for the impact of HES5 silencing on tumourigenesis as a starting point for future functional studies.

Origins and functional consequences of somatic mitochondrial DNA mutations in human cancer.

Recent sequencing studies have extensively explored the somatic alterations present in the nuclear genomes of cancers. Although mitochondria control energy metabolism and apoptosis, the origins and impact of cancer-associated mutations in mtDNA are unclear. In this study, we analyzed somatic alterations in mtDNA from 1675 tumors. We identified 1907 somatic substitutions, which exhibited dramatic replicative strand bias, predominantly C > T and A > G on the mitochondrial heavy strand. This strand-asymmetric signature differs from those found in nuclear cancer genomes but matches the inferred germline process shaping primate mtDNA sequence content. A number of mtDNA mutations showed considerable heterogeneity across tumor types. Missense mutations were selectively neutral and often gradually drifted towards homoplasmy over time. In contrast, mutations resulting in protein truncation undergo negative selection and were almost exclusively heteroplasmic. Our findings indicate that the endogenous mutational mechanism has far greater impact than any other external mutagens in mitochondria and is fundamentally linked to mtDNA replication.

Inactivating CUX1 mutations promote tumorigenesis.

A major challenge in cancer genetics is to determine which low-frequency somatic mutations are drivers of tumorigenesis. Here we interrogate the genomes of 7,651 diverse human cancers and find inactivating mutations in the homeodomain transcription factor gene CUX1 (cut-like homeobox 1) in ~1-5% of various tumors. Meta-analysis of CUX1 mutational status in 2,519 cases of myeloid malignancies reveals disruptive mutations associated with poor survival, highlighting the clinical significance of CUX1 loss. In parallel, we validate CUX1 as a bona fide tumor suppressor using mouse transposon-mediated insertional mutagenesis and Drosophila cancer models. We demonstrate that CUX1 deficiency activates phosphoinositide 3-kinase (PI3K) signaling through direct transcriptional downregulation of the PI3K inhibitor PIK3IP1 (phosphoinositide-3-kinase interacting protein 1), leading to increased tumor growth and susceptibility to PI3K-AKT inhibition. Thus, our complementary approaches identify CUX1 as a pan-driver of tumorigenesis and uncover a potential strategy for treating CUX1-mutant tumors.

Mapping protein-DNA interactions using ChIP-sequencing.

Chromatin immunoprecipitation (ChIP) allows enrichment of genomic regions which are associated with specific transcription factors, histone modifications, and indeed any other epitopes which are present on chromatin. The original ChIP methods used site-specific PCR and Southern blotting to confirm which regions of the genome were enriched, on a candidate basis. The combination of ChIP with genomic tiling arrays (ChIP-chip) allowed a more unbiased approach to map ChIP-enriched sites. However, limitations of microarray probe design and probe number have a detrimental impact on the coverage, resolution, sensitivity, and cost of whole-genome tiling microarray sets for higher eukaryotes with large genomes. The combination of ChIP with high-throughput sequencing technology has allowed more comprehensive surveys of genome occupancy, greater resolution, and lower cost for whole genome coverage. Herein, we provide a comparison of high-throughput sequencing platforms and a survey of ChIP-seq analysis tools, discuss experimental design, and describe a detailed ChIP-seq method.Chromatin immunoprecipitation (ChIP) allows enrichment of genomic regions which are associated with specific transcription factors, histone modifications, and indeed any other epitopes which are present on chromatin. The original ChIP methods used site-specific PCR and Southern blotting to confirm which regions of the genome were enriched, on a candidate basis. The combination of ChIP with genomic tiling arrays (ChIP-chip) allowed a more unbiased approach to map ChIP-enriched sites. However, limitations of microarray probe design and probe number have a detrimental impact on the coverage, resolution, sensitivity, and cost of whole-genome tiling microarray sets for higher eukaryotes with large genomes. The combination of ChIP with high-throughput sequencing technology has allowed more comprehensive surveys of genome occupancy, greater resolution, and lower cost for whole genome coverage. Herein, we provide a comparison of high-throughput sequencing platforms and a survey of ChIP-seq analysis tools, discuss experimental design, and describe a detailed ChIP-seq method.

Global identification of androgen response elements.

Chromatin immunoprecipitation (ChIP) is an invaluable tool in the study of transcriptional regulation. ChIP methods require both a priori knowledge of the transcriptional regulators which are important for a given biological system and high-quality specific antibodies for these targets. The androgen receptor (AR) is known to play essential roles in male sexual development, in prostate cancer and in the function of many other AR-expressing cell types (e.g. neurons and myocytes). As a ligand-activated transcription factor the AR also represents an endogenous, inducible system to study transcriptional biology. Therefore, ChIP studies of the AR can make use of treatment contrast experiments to define its transcriptional targets. To date several studies have mapped AR binding sites using ChIP in combination with genome tiling microarrays (ChIP-chip) or direct sequencing (ChIP-seq), mainly in prostate cancer cell lines and with varying degrees of genomic coverage. These studies have provided new insights into the DNA sequences to which the AR can bind, identified AR cooperating transcription factors, mapped thousands of potential AR regulated genes and provided insights into the biological processes regulated by the AR. However, further ChIP studies will be required to fully characterise the dynamics of the AR-regulated transcriptional programme, to map the occupancy of different AR transcriptional complexes which result in different transcriptional output and to delineate the transcriptional networks downstream of the AR.

The androgen receptor fuels prostate cancer by regulating central metabolism and biosynthesis.

The androgen receptor (AR) is a key regulator of prostate growth and the principal drug target for the treatment of prostate cancer. Previous studies have mapped AR targets and identified some candidates which may contribute to cancer progression, but did not characterize AR biology in an integrated manner. In this study, we took an interdisciplinary approach, integrating detailed genomic studies with metabolomic profiling and identify an anabolic transcriptional network involving AR as the core regulator. Restricting flux through anabolic pathways is an attractive approach to deprive tumours of the building blocks needed to sustain tumour growth. Therefore, we searched for targets of the AR that may contribute to these anabolic processes and could be amenable to therapeutic intervention by virtue of differential expression in prostate tumours. This highlighted calcium/calmodulin-dependent protein kinase kinase 2, which we show is overexpressed in prostate cancer and regulates cancer cell growth via its unexpected role as a hormone-dependent modulator of anabolic metabolism. In conclusion, it is possible to progress from transcriptional studies to a promising therapeutic target by taking an unbiased interdisciplinary approach.