Katelyn Mullen

The genomic features of pancreatic ductal adenocarcinoma (PDAC) have been well described, yet the evolutionary contexts within which these features occur remain unexplored. We studied genome landscapes, phylogenies, and clonal compositions of 91 PDACs in relation to clinicopathologic features. There was no difference in the number of driver mutations or evolutionary timing when each mutation occurred. High truncal density, a metric of the accumulation of somatic mutations in the lineage that gave rise to each PDAC, was significantly associated with worse overall survival. Polyclonal, monoclonal, or mixed polyclonal/monoclonal metastases were identified across the cohort, highlighting multiple forms of intertumoral heterogeneity. Advanced stage and treated PDACs had higher odds of being polyclonal, whereas oligometastatic PDACs had fewer driver alterations, a lower fractional allelic loss, and increased likelihood of being monoclonal. In sum, our findings reveal novel insights into the dynamic nature of the PDAC genome beyond established genetic paradigms. Significance: Although the pancreatic cancer genome has been described, it has not been explored with respect to stages of diagnosis or treatment bottlenecks. We now describe and quantify the genomic features of PDAC in the context of evolutionary metrics and in doing so have identified a novel prognostic biomarker.

Multiple large-scale genomic profiling efforts have been undertaken in osteosarcoma to define the genomic drivers of tumorigenesis, therapeutic response, and disease recurrence. The spatial and temporal intratumor heterogeneity could also play a role in promoting tumor growth and treatment resistance. We conducted longitudinal whole-genome sequencing of 37 tumor samples from 8 patients with relapsed or refractory osteosarcoma. Each patient had at least one sample from a primary site and a metastatic or relapse site. Subclonal copy-number alterations were identified in all patients except one. In 5 patients, subclones from the primary tumor emerged and dominated at subsequent relapses. MYC gain/amplification was enriched in the treatment-resistant clones in 6 of 7 patients with multiple clones. Amplifications in other potential driver genes, such as CCNE1, RAD21, VEGFA, and IGF1R, were also observed in the resistant copy-number clones. A chromosomal duplication timing analysis revealed that complex genomic rearrangements typically occurred prior to diagnosis, supporting a macroevolutionary model of evolution, where a large number of genomic aberrations are acquired over a short period of time followed by clonal selection, as opposed to ongoing evolution. A mutational signature analysis of recurrent tumors revealed that homologous repair deficiency (HRD)-related SBS3 increases at each time point in patients with recurrent disease, suggesting that HRD continues to be an active mutagenic process after diagnosis. Overall, by examining the clonal relationships between temporally and spatially separated samples from patients with relapsed/refractory osteosarcoma, this study sheds light on the intratumor heterogeneity and potential drivers of treatment resistance in this disease. Significance: The chemoresistant population in recurrent osteosarcoma is subclonal at diagnosis, emerges at the time of primary resection due to selective pressure from neoadjuvant chemotherapy, and is characterized by unique oncogenic amplifications.

Despite insights gained by bulk DNA sequencing of cancer it remains challenging to resolve the admixture of normal and tumor cells, and/or of distinct tumor subclones; high-throughput single-cell DNA sequencing circumvents these and brings cancer genomic studies to higher resolution. However, its application has been limited to liquid tumors or a small batch of solid tumors, mainly because of the lack of a scalable workflow to process solid tumor samples. Here we optimize a highly automated nuclei extraction workflow that achieves fast and reliable targeted single-nucleus DNA library preparation of 38 samples from 16 pancreatic ductal adenocarcinoma patients, with an average library yield per sample of 2867 single nuclei. We demonstrate that this workflow not only performs well using low cellularity or low tumor purity samples but reveals genomic evolution patterns of pancreatic ductal adenocarcinoma as well.

Multiple large-scale tumor genomic profiling efforts have been undertaken in osteosarcoma, however, little is known about the spatial and temporal intratumor heterogeneity and how it may drive treatment resistance. We performed whole-genome sequencing of 37 tumor samples from eight patients with relapsed or refractory osteosarcoma. Each patient had at least one sample from a primary site and a metastatic or relapse site. We identified subclonal copy number alterations in all but one patient. We observed that in five patients, a subclonal copy number clone from the primary tumor emerged and dominated at subsequent relapses. MYC gain/amplification was enriched in the treatment-resistant clone in 6 out of 7 patients with more than one clone. Amplifications in other potential driver genes, such as CCNE1 , RAD21 , VEGFA , and IGF1R , were also observed in the resistant copy number clones. Our study sheds light on intratumor heterogeneity and the potential drivers of treatment resistance in osteosarcoma. Significance: Subclonal copy number clones emerged and dominated in relapsed osteosarcoma, with MYC gain/amplification being the defining characteristic in our cohort. Selective pressure from neoadjuvant chemotherapy revealed this clone at the time of primary resection, highlighting that genomic profiling at this time may identify clones that are selected for, or determine innate resistance to primary chemotherapy.

Despite insights gained by bulk DNA sequencing of cancer it remains challenging to resolve the admixture of normal and tumor cells, and/or of distinct tumor subclones; high throughput single-cell DNA sequencing circumvents these and brings cancer genomic studies to higher resolution. However, its application has been limited to liquid tumors or a small batch of solid tumors, mainly because of the lack of a scalable workflow to process solid tumor samples. Here we optimized a highly automated nuclei extraction workflow that achieved fast and reliable targeted single-nucleus DNA library preparation of 38 samples from 16 pancreatic adenocarcinoma (PDAC) patients, with an average library yield per sample of 2867 single nuclei. We demonstrate that this workflow not only performs well using low cellularity or low tumor purity samples but reveals novel genomic evolution patterns of PDAC as well.

Despite insights gained by bulk DNA sequencing of cancer it remains challenging to resolve the admixture of normal and tumor cells, and/or of distinct tumor subclones; high throughput single-cell DNA sequencing circumvents these and brings cancer genomic studies to higher resolution. However, its application has been limited to liquid tumors or a small batch of solid tumors, mainly because of the lack of a scalable workflow to process solid tumor samples. Here we optimized a highly automated nuclei extraction workflow that achieved fast and reliable targeted single-nucleus DNA library preparation of 38 samples from 16 pancreatic adenocarcinoma (PDAC) patients, with an average library yield per sample of 2867 single nuclei. We demonstrate that this workflow not only performs well using low cellularity or low tumor purity samples but reveals novel genomic evolution patterns of PDAC as well.