Multiregion sampling of de novo metastatic prostate cancer reveals complex polyclonality and augments clinical genotyping

De novo metastatic prostate cancer is highly aggressive, but the paucity of routinely collected tissue has hindered genomic stratification and precision oncology. Here, we leveraged a rare study of surgical intervention in 43 de novo metastatic prostate cancers to assess somatic genotypes across 607 synchronous primary and metastatic tissue regions plus circulating tumor DNA. Intra-prostate heterogeneity was pervasive and impacted clinically relevant genes, resulting in discordant genotypes between select primary restricted regions and synchronous metastases. Additional complexity was driven by polyclonal metastatic seeding from phylogenetically related primary populations. When simulating clinical practice relying on a single tissue region, genomic heterogeneity plus variable tumor fraction across samples caused inaccurate genotyping of dominant disease; however, pooling extracted DNA from multiple biopsy cores before sequencing can rescue misassigned somatic genotypes. Our results define the relationship between synchronous treatment-sensitive primary and metastatic lesions in men with de novo metastatic prostate cancer and provide a framework for implementing genomics-guided patient management.

The diverse effects of phenotypic dominance on hybrid fitness

When divergent populations interbreed, their alleles are brought together in hybrids. In the initial F1 cross, most divergent loci are heterozygous. Therefore, F1 fitness can be influenced by dominance effects that could not have been selected to function well together. We present a systematic study of these F1 dominance effects by introducing variable phenotypic dominance into Fisher's geometric model. We show that dominance often reduces hybrid fitness, which can generate optimal outbreeding followed by a steady decline in F1 fitness, as is often observed. We also show that "lucky" beneficial effects sometimes arise by chance, which might be important when hybrids can access novel environments. We then show that dominance can lead to violations of Haldane's Rule (reduced fitness of the heterogametic F1) but strengthens Darwin's Corollary (F1 fitness differences between cross directions). Taken together, results show that the effects of dominance on hybrid fitness can be surprisingly difficult to isolate, because they often resemble the effects of uniparental inheritance or expression. Nevertheless, we identify a pattern of environment-dependent heterosis that only dominance can explain, and for which there is some suggestive evidence. Our results also show how existing data set upper bounds on the size of dominance effects. These bounds could explain why additive models often provide good predictions for later-generation recombinant hybrids, even when dominance qualitatively changes outcomes for the F1.

Deep whole-genome ctDNA chronology of treatment-resistant prostate cancer

Circulating tumour DNA (ctDNA) in blood plasma is an emerging tool for clinical cancer genotyping and longitudinal disease monitoring¹. However, owing to past emphasis on targeted and low-resolution profiling approaches, our understanding of the distinct populations that comprise bulk ctDNA is incomplete^2-12. Here we perform deep whole-genome sequencing of serial plasma and synchronous metastases in patients with aggressive prostate cancer. We comprehensively assess all classes of genomic alterations and show that ctDNA contains multiple dominant populations, the evolutionary histories of which frequently indicate whole-genome doubling and shifts in mutational processes. Although tissue and ctDNA showed concordant clonally expanded cancer driver alterations, most individual metastases contributed only a minor share of total ctDNA. By comparing serial ctDNA before and after clinical progression on potent inhibitors of the androgen receptor (AR) pathway, we reveal population restructuring converging solely on AR augmentation as the dominant genomic driver of acquired treatment resistance. Finally, we leverage nucleosome footprints in ctDNA to infer mRNA expression in synchronously biopsied metastases, including treatment-induced changes in AR transcription factor signalling activity. Our results provide insights into cancer biology and show that liquid biopsy can be used as a tool for comprehensive multi-omic discovery.

Abstract 3625: Clonal architecture and evolution of treatment-resistant prostate cancer via deep whole-genome ctDNA sequencing

Background: Circulating tumor DNA (ctDNA) in blood plasma is an emerging tool for clinical cancer genotyping and longitudinal disease monitoring. However, integration of ctDNA tests into clinical management is critically impeded by the poor understanding of the distinct somatic populations comprising bulk ctDNA—including their relationship to synchronous metastatic tissue, and their temporal dynamics during standard-of-care treatment. Prior approaches relying on targeted and/or low-resolution techniques (e.g. targeted exon sequencing, low-pass (shallow) whole-genome sequencing; WGS) do not permit comprehensive dissection of clonal architecture and unbiased analysis of putative resistance mechanisms.

Methods: We performed deep WGS on serial plasma ctDNA (median depth: 185×) and synchronous metastatic tissue biopsies with high tumor purity from 35 patients with metastatic castration-resistant prostate cancer. We developed a subclonal reconstruction algorithm optimized for our data enabling resolution of per-patient evolutionary histories and ctDNA clonal composition. ctDNA nucleosome footprinting was used to infer mRNA abundance in synchronously biopsied metastases and androgen receptor (AR) transcription factor activity at 3224 AR binding sites (ARBS).

Results: We comprehensively assess all classes of genomic alterations and demonstrate that ctDNA harbors greater populational heterogeneity than metastatic tissue (p<0.001). The evolutionary histories of ctDNA populations indicate frequent whole-genome doubling and attenuation of C>T aging-associated mutation signature during subclonal differentiation. Although driver alterations were largely concordant between tissue and ctDNA, each individual metastasis contributed only a minor share of total ctDNA (average ctDNA contribution: 17%). By comparing serial ctDNA before and after clinical progression on potent AR pathway inhibitors, we reveal population restructuring converging solely on AR copy augmentation as the dominant genomic driver of acquired treatment-resistance. Nucleosome depletion at transcription start-sites is highly correlated with same-patient metastatic tissue mRNA abundance, indicating that ctDNA fragmentomics can recapitulate transcriptomic patterns in metastatic lesions. Most ctDNA samples exhibited strong ARBS nucleosome depletion which correlated with AR gene copy number (R=0.36, p=0.003). Finally, serial ctDNA nucleosome profiling at ARBS revealed adaptive transcriptomic resistance to AR pathway inhibitors, including lineage switch to a neuroendocrine-like (AR-low) state.

Conclusions: We show that the populations comprising ctDNA are typically complex and more heterogeneous than those found in bulk WGS of a synchronous metastasis. Our work advocates for liquid biopsy as a comprehensive multi-omic discovery tool for cancers with high ctDNA fractions.

Citation Format: Cameron Herberts, Matti Annala, Joonatan Sipola, Sarah W. Ng, Xinyi E. Chen, Anssi Nurminen, Olga Korhonen, Aslı D. Munzur, Kevin Beja, Elena Schönlau, Cecily Q. Bernales, Elie Ritch, Jack V. Bacon, Nathan A. Lack, Matti Nykter, Rahul Aggarwal, Eric J. Small, Martin E. Gleave, David A. Quigley, Felix Y. Feng, Kim N. Chi, Alexander W. Wyatt. Clonal architecture and evolution of treatment-resistant prostate cancer via deep whole-genome ctDNA sequencing [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 3625.

Adult-type granulosa cell tumor of the ovary: a FOXL2-centric disease

Adult-type granulosa cell tumors (aGCTs) account for 90% of malignant ovarian sex cord-stromal tumors and 2-5% of all ovarian cancers. These tumors are usually diagnosed at an early stage and are treated with surgery. However, one-third of patients relapse between 4 and 8 years after initial diagnosis, and there are currently no effective treatments other than surgery for these relapsed patients. As the majority of aGCTs (>95%) harbor a somatic mutation in FOXL2 (c.C402G; p.C134W), the aim of this study was to identify genetic mutations besides FOXL2 C402G in aGCTs that could explain the clinical diversity of this disease. Whole-genome sequencing of 10 aGCTs and their matched normal blood was performed to identify somatic mutations. From this analysis, a custom amplicon-based panel was designed to sequence 39 genes of interest in a validation cohort of 83 aGCTs collected internationally. KMT2D inactivating mutations were present in 10 of 93 aGCTs (10.8%), and the frequency of these mutations was similar between primary and recurrent aGCTs. Inactivating mutations, including a splice site mutation in candidate tumor suppressor WNK2 and nonsense mutations in PIK3R1 and NLRC5, were identified at a low frequency in our cohort. Missense mutations were identified in cell cycle-related genes TP53, CDKN2D, and CDK1. From these data, we conclude that aGCTs are comparatively a homogeneous group of tumors that arise from a limited set of genetic events and are characterized by the FOXL2 C402G mutation. Secondary mutations occur in a subset of patients but do not explain the diverse clinical behavior of this disease. As the FOXL2 C402G mutation remains the main driver of this disease, progress in the development of therapeutics for aGCT would likely come from understanding the functional consequences of the FOXL2 C402G mutation.

Single cell transcriptomes of normal endometrial derived organoids uncover novel cell type markers and cryptic differentiation of primary tumours

Endometrial carcinoma, the most common gynaecological cancer, develops from endometrial epithelium which is composed of secretory and ciliated cells. Pathologic classification is unreliable and there is a need for prognostic tools. We used single cell sequencing to study organoid model systems derived from normal endometrial endometrium to discover novel markers specific for endometrial ciliated or secretory cells. A marker of secretory cells (MPST) and several markers of ciliated cells (FAM92B, WDR16, and DYDC2) were validated by immunohistochemistry on organoids and tissue sections. We performed single cell sequencing on endometrial and ovarian tumours and found both secretory-like and ciliated-like tumour cells. We found that ciliated cell markers (DYDC2, CTH, FOXJ1, and p73) and the secretory cell marker MPST were expressed in endometrial tumours and positively correlated with disease-specific and overall survival of endometrial cancer patients. These findings suggest that expression of differentiation markers in tumours correlates with less aggressive disease, as would be expected for tumours that retain differentiation capacity, albeit cryptic in the case of ciliated cells. These markers could be used to improve the risk stratification of endometrial cancer patients, thereby improving their management. We further assessed whether consideration of MPST expression could refine the ProMiSE molecular classification system for endometrial tumours. We found that higher expression levels of MPST could be used to refine stratification of three of the four ProMiSE molecular subgroups, and that any level of MPST expression was able to significantly refine risk stratification of the copy number high subgroup which has the worst prognosis. Taken together, this shows that single cell sequencing of putative cells of origin has the potential to uncover novel biomarkers that could be used to guide management of cancers. © 2020 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.