ABSTRACT
Advanced prostate cancers comprise distinct phenotypes, but tumor classification remains clinically challenging. Here, we harnessed circulating tumor DNA (ctDNA) to study tumor phenotypes by ascertaining nucleosome positioning patterns associated with transcription regulation. We sequenced plasma ctDNA whole genomes from patient-derived xenografts representing a spectrum of androgen receptor active (ARPC) and neuroendocrine (NEPC) prostate cancers. Nucleosome patterns associated with transcriptional activity were reflected in ctDNA at regions of genes, promoters, histone modifications, transcription factor binding, and accessible chromatin. We identified the activity of key phenotype-defining transcriptional regulators from ctDNA, including AR, ASCL1, HOXB13, HNF4G, and GATA2. To distinguish NEPC and ARPC in patient plasma samples, we developed prediction models that achieved accuracies of 97% for dominant phenotypes and 87% for mixed clinical phenotypes. Although phenotype classification is typically assessed by IHC or transcriptome profiling from tumor biopsies, we demonstrate that ctDNA provides comparable results with diagnostic advantages for precision oncology. SIGNIFICANCE: This study provides insights into the dynamics of nucleosome positioning and gene regulation associated with cancer phenotypes that can be ascertained from ctDNA. New methods for classification in phenotype mixtures extend the utility of ctDNA beyond assessments of somatic DNA alterations with important implications for molecular classification and precision oncology. This article is highlighted in the In This Issue feature, p. 517.
Subject(s)
Circulating Tumor DNA , Prostatic Neoplasms , Male , Humans , Circulating Tumor DNA/genetics , Nucleosomes/genetics , Precision Medicine , Prostatic Neoplasms/pathology , Gene Expression Regulation, Neoplastic , PhenotypeABSTRACT
Cancers with homology-directed DNA repair (HRR) deficiency exhibit high response rates to poly(ADP-ribose) polymerase inhibitors (PARPi) and platinum chemotherapy. Though mutations disrupting BRCA1 and BRCA2 associate with HRR deficiency (HRRd), patterns of genomic aberrations and mutation signatures may be more sensitive and specific indicators of compromised repair. Here, we evaluated whole-exome sequences from 418 metastatic prostate cancers (mPCs) and determined that one-fifth exhibited genomic characteristics of HRRd that included Catalogue Of Somatic Mutations In Cancer mutation signature 3. Notably, a substantial fraction of tumors with genomic features of HRRd lacked biallelic loss of a core HRR-associated gene, such as BRCA2. In this subset, HRRd associated with loss of chromodomain helicase DNA binding protein 1 but not with mutations in serine-protein kinase ATM, cyclin dependent kinase 12, or checkpoint kinase 2. HRRd genomic status was strongly correlated with responses to PARPi and platinum chemotherapy, a finding that supports evaluating biomarkers reflecting functional HRRd for treatment allocation.