Maley, Carlo C.

Abstract Pancreatic ductal adenocarcinoma (PDAC) is predicted to be the second cause of cancer death in the US by 2040. Five-year overall survival is only ~12% despite decades of research. The mainstay of treatment remains chemotherapy in the absence of highly recurrent, therapeutically actionable, genetic targets. For the minority of patients (~15-20%) who are candidates for curative-intent surgical resection, six months of intensive neoadjuvant chemotherapy with/without chemoradiation has become standard-of-care. However, despite evidence of radiographic, metabolic, and biochemical responses, significant viable residual disease on pathologic analysis is almost universally observed with approximately 40-60% risk of disease recurrence. Thus, there is a critical need to better understand why PDAC remains recalcitrant to current therapeutic modalities. Pathologic evaluation of PDAC reveals a complex microenvironment consisting of neoplastic cells dynamically interacting with both “normal” cellular (e.g. mesenchymal, immune) and non-cellular (e.g. collagen) stromal elements. While biologically important, experimental models and small-scale efforts to dissect the PDAC ecosystem are limited by focus on individual tumor compartments and generation of conflicting results. Consequently, there is a need for a more comprehensive understanding of the structure and dynamics of the PDAC TME through the theoretical lens of evolutionary ecology. Current therapies are not focused on the PDAC ecosystem. We hypothesize that cell-cell interactions of PDAC cells with its local niche are important for conferring treatment resistance. Thus, better understanding of said interactions could yield novel ecologically-informed therapies. Consequently, we developed a digital pathology library of approximately 1000 PDAC resection cases (whole-slide H&E images) performed at Mayo Clinic between 2013-2021. Approximately 400 cases are treatment naive while about 600 were treated with neoadjuvant therapy. While some qualitative features of pathologic analysis have prognostic value (e.g., pathologic treatment response, presence of perineural or lymphovascular invasion, extent of desmoplasia), we aim to identify prognostic spatial relationships using tools from landscape ecology. We utilized the HoVer-Net model for nuclei segmentation and classification and spatial statistics methods and landscape fragmentation analysis to understand the spatial organization within the PDAC tumor ecology. With these methodologies, we can assess the spatial context of cell interactions, habitat/niche/microenvironment formation, ecological relationships between cell populations, and the different tissue components. The different spatial statistics/metrics are used to quantify levels of immune cell infiltration, mixing, as well as the complexity and diversity of microenvironment phenotypes leading to high levels of tumor heterogeneity. Such measurements can help predict the likelihood of a PDAC tumor adaptive potential that could lead to fast evolutionary responses to therapy, progression, and poor prognosis. Citation Format: Merih D. Toruner, Luis Cisneros, Zafar Siddiqui, Chady Meroueh, Ritambhara Singh, Carlo Maley, Martin E. Fernandez-Zapico, Ryan M. Carr. Dissecting the pancreatic cancer microenvironment using landscape ecology [abstract]. In: Proceedings of the AACR Special Conference in Cancer Research: Translating Cancer Evolution and Data Science: The Next Frontier; 2023 Dec 3-6; Boston, Massachusetts. Philadelphia (PA): AACR; Cancer Res 2024;84(3 Suppl_2):Abstract nr B026.

Adaptive therapy, an ecologically inspired approach to cancer treatment, aims to overcome resistance and reduce toxicity by leveraging competitive interactions between drug-sensitive and drug-resistant subclones, prioritizing patient survival and quality of life instead of killing the maximum number of cancer cells. In preparation for a clinical trial, we used endocrine-resistant MCF7 breast cancer to stimulate second-line therapy and tested adaptive therapy using capecitabine, gemcitabine, or their combination in a mouse xenograft model. Dose modulation adaptive therapy with capecitabine alone increased survival time relative to MTD but not statistically significantly (HR = 0.22, 95% CI = 0.043–1.1, p = 0.065). However, when we alternated the drugs in both dose modulation (HR = 0.11, 95% CI = 0.024–0.55, p = 0.007) and intermittent adaptive therapies, the survival time was significantly increased compared to high-dose combination therapy (HR = 0.07, 95% CI = 0.013–0.42, p = 0.003). Overall, the survival time increased with reduced dose for both single drugs (p < 0.01) and combined drugs (p < 0.001), resulting in tumors with fewer proliferation cells (p = 0.0026) and more apoptotic cells (p = 0.045) compared to high-dose therapy. Adaptive therapy favors slower-growing tumors and shows promise in two-drug alternating regimens instead of being combined.

Abstract PURPOSE To examine incremental values of estrogen receptor (ER) status, body mass index (BMI), menopausal status, and a previously reported multi-gene classifier over commonly used clinical factors (i.e. age, tumor grade, comedonecrosis, surgical margins, and treatment) in predicting risk of any ipsilateral recurrence (IR) event within five years after DCIS diagnosis. METHODS A derivation cohort consisted of participants in the Translational Breast Cancer Research Consortium (TBCRC) 038, a retrospective multicenter cohort study in women undergoing surgical resection for DCIS between 01/01/1998 and 02/29/2016 (n=216). The validation cohort, the Repository of Archival Human Breast Tissue (RAHBT) at Washington University School of Medicine, provided cases meeting the same eligibility criteria as TBCRC038 (n=97). Participants in both cohorts had RNA-seq data and either developed IR 1-5y after initial DCIS diagnosis or were free from subsequent breast events for at least five years. The previously reported 812-gene classifier had been developed from a subset of the TBCRC038 samples using a negative-binomial regression model to identify differentially expressed genes in the primary tumor associated with subsequent recurrence events. This classifier has been shown to be highly correlated with 5-year invasive, DCIS, and all breast cancer events, and validated in the RAHBT cohort. Cox proportional hazards regression was used to estimate hazard ratios (HRs) of IR in the TBCRC038 cohort (76 with IR). The clinical score was developed using clinical predictors (aforementioned clinical factors and ER) and their regression coefficients from the model with the maximum predictive accuracy (e.g. c-index) and the minimum number of predictors; the summary score integrated the clinical score and multi-gene classifier. Predictive performance of both clinical and summary scores was validated in the RAHBT cohort (20 with IR). RESULTS In the TBCRC cohort derivation set, we used a multivariable model based on clinical factors alone (clinical score) and found that ER status, but not BMI or menopausal status, was independently associated with a higher IR risk (HR=2.06, 95% CI 1.18-3.58). Adding the multi-gene classifier to the clinical factors-based model (summary score) in the TBCRC038 test set increased predictive accuracy (c-index 0.68 to 0.70), with the genomic classifier-adjusted HR of 14.96 (95% CI 8.64-25.91). The summary score had higher predictive performance for IR risk than clinical score alone (c-index 0.82 vs. 0.70). In the RAHBT validation samples, model performance was similarly improved using summary scores clinical factors-based model plus multigene classifier as compared to clinical scores along (c-index 0.74 vs. 0.58). CONCLUSION Combining clinical factors and a multigene classifier provided more accurate risk estimates of IR within five years after excision of DCIS than clinical factors alone. Figure 1. Observed and predicted recurrence-free survival in the first five years after initial DCIS diagnosis in the RAHBT validation cohort, by risk groups defined by clinical scores (left) and clinical score plus multigene classifier (right). Citation Format: Ying Liu, Siri H. Strand, Lorraine King, Bryan Harmon, Fergus J. Couch, Kristalyn Gallagher, Mark Kilgore, Shi Wei, Angela DeMichele, Tari King, Priscilla F. McAuliffe, Jeffrey Marks, Carlo Maley, Robert West, E Shelley Hwang, Graham A. Colditz. Using clinical characteristics and molecular markers to predict the risk of subsequent ipsilateral breast events after excision of DCIS [abstract]. In: Proceedings of the 2022 San Antonio Breast Cancer Symposium; 2022 Dec 6-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2023;83(5 Suppl):Abstract nr P1-07-02.

AbstractUnderstanding intratumor heterogeneity is critical for studying tumorigenesis and designing personalized treatments. To decompose the mixed cell population in a tumor, subclones are inferred computationally based on variant allele frequency (VAF) from bulk sequencing data. In this study, we showed that sequencing depth, mean VAF, and variance of VAF of a subclone are confounded. Without considering this effect, current methods require deep-sequencing data (&gt;300× depth) to reliably infer subclones. Here, we present a novel algorithm that incorporates depth-variance and mean-variance dependencies in a clustering error model and successfully identifies subclones in tumors sequenced at depths of as low as 30×. We implemented the algorithm as a model-based adaptive grouping of subclones (MAGOS) method. Analyses of computer simulated data and empirical sequencing data showed that MAGOS outperformed existing methods on minimum sequencing depth, decomposition accuracy, and computation efficiency. The most prominent improvements were observed in analyzing tumors sequenced at depths between 30× and 200×, whereas the performance was comparable between MAGOS and existing methods on deeply sequenced tumors. MAGOS supports analysis of single-nucleotide variants and copy number variants from a single sample or multiple samples of a tumor. We applied MAGOS to whole-exome data of late-stage liver cancers and discovered that high subclone count in a tumor was a significant risk factor of poor prognosis. Lastly, our analysis suggested that sequencing multiple samples of the same tumor at standard depth is more cost-effective and robust for subclone characterization than deep sequencing a single sample. MAGOS is available at github (https://github.com/liliulab/magos).

Abstract Background. DCIS consists of a molecularly heterogeneous group of premalignant lesions, with variable risk of invasive progression. Understanding biomarkers for invasive progression could help individualize treatment recommendations based upon tumor biology. As part of the NCI Human Tumor Atlas Network (HTAN), we conducted comprehensive genomic analyses on two large DCIS case-control cohorts. Methods. We performed smart3-seq and low-pass whole genome sequencing on two independent, retrospective, longitudinally sampled DCIS case-control cohorts. TBCRC 038 was a multicenter cohort diagnosed with DCIS between 1998 and 2016 at one of the Translational Breast Cancer Research sites; the RAHBT (Resource of Archival Human Breast Tissue) cohort included women identified through the St. Louis Breast Tissue Repository, and the Women’s Health Repository diagnosed between 1997 and 2001. We studied the spectrum of molecular changes present and sought genomic predictors of subsequent ipsilateral breast events (iBEs: DCIS recurrence or invasive progression) in both DCIS epithelium and stroma in formalin fixed paraffin embedded tissue. We generated de novo tumor and stroma-centric subtypes for DCIS that represents fundamental transcriptomic organization. Copy number analysis was performed using low-pass DNA sequencing. Non-negative matrix factorization (NMF) was applied to the RNA expression of all coding genes to identify clusters. A negative-binomial regression model was used to identify differentially expressed genes. Results. We analyzed 677 DCIS samples from 481 patients with 7.1 years median follow-up. In TBCRC samples, we identified three clusters via NMF in TBCRC referred to as ER low, quiescent, and ER high. The ER-low cluster had significantly higher levels of ERBB2 and lower levels of ESR1 compared to quiescent and ER-high clusters. Quiescent cluster lesions were less proliferative and less metabolically active than ER high and ER low subtypes. These findings were replicated in the RAHBT cohort. Focusing on the stromal component of DCIS from laser capture microdissection in RAHBT samples, we identified four distinct DCIS-associated stromal clusters. A “normal-like” stromal cluster with ECM organization and PI3K-AKT signaling; a “collagen-rich” stromal cluster; a “desmoplastic” stromal cluster with high fibroblast and total myeloid abundance, mostly associated with macrophages and myeloid dendritic cells (mDC); and an “immune-dense” stromal cluster. Further, we compared differentially expressed genes in patients with or without subsequent iBEs within 5 years of diagnosis. Hypothesizing that the resulting 812 DE genes (DESeq2) represent multiple routes to subsequent iBEs, we leveraged NMF to identify paths to progression. In both TBCRC and RAHBT cohorts, poor outcome groups exhibited increased ER, MYC signaling, and oxidative phosphorylation, supporting that these pathways are important for DCIS recurrence and progression. Conclusion. Comprehensive genomic profiling in two independent DCIS cohorts with longitudinal outcomes shows distinct DCIS stromal expression patterns and immune cell composition. RNA expression profiles reveal underlying tumor biology that is associated with later iBEs in both cohorts. These studies provide new insight into DCIS biology and will guide the design of diagnostic strategies to prevent invasive progression. Citation Format: Siri H Strand, Belén Rivero-Gutiérrez, Kathleen E Houlahan, Jose A Seoane, Lorraine M King, Tyler Risom, Lunden Simpson, Sujay Vennam, Aziz Khan, Timothy Hardman, Bryan E Harmon, Fergus J Couch, Kristalyn Gallagher, Mark Kilgore, Shi Wei, Angela DeMichele, Tari King, Priscilla F McAuliffe, Julie Nangia, Joanna Lee, Jennifer Tseng, Anna Maria Storniolo, Alastair Thompson, Gaorav Gupta, Robyn Burns, Deborah J Veis, Katherine DeSchryver, Chunfang Zhu, Magdalena Matusiak, Jason Wang, Shirley X Zhu, Jen Tappenden, Daisy Yi Ding, Dadong Zhang, Jingqin Luo, Shu Jiang, Sushama Varma, Cody Straub, Sucheta Srivastava, Christina Curtis, Rob Tibshirani, Robert Michael Angelo, Allison Hall, Kouros Owzar, Kornelia Polyak, Carlo Maley, Jeffrey R Marks, Graham A Colditz, E Shelley Hwang, Robert B West. The Breast PreCancer Atlas DCIS genomic signatures define biology and correlate with clinical outcomes: An analysis of TBCRC 038 and RAHBT cohorts [abstract]. In: Proceedings of the 2021 San Antonio Breast Cancer Symposium; 2021 Dec 7-10; San Antonio, TX. Philadelphia (PA): AACR; Cancer Res 2022;82(4 Suppl):Abstract nr GS4-07.

We review and synthesize key findings and limitations of mathematical models for prostate cancer, both from theoretical work and data-validated approaches, especially concerning clinical applications. Our focus is on models of prostate cancer dynamics under treatment, particularly with a view toward optimizing hormone-based treatment schedules and estimating the onset of treatment resistance under various assumptions. Population models suggest that intermittent or adaptive therapy is more beneficial to delay cancer relapse as compared to the standard continuous therapy if treatment resistance comes at a competitive cost for cancer cells. Another consensus among existing work is that the standard biomarker for cancer growth, prostate-specific antigen, may not always correlate well with cancer progression. Instead, its doubling rate appears to be a better indicator of tumor growth. Much of the existing work utilizes simple ordinary differential equations due to difficulty in collecting spatial data and due to the early success of using prostate-specific antigen in mathematical modeling. However, a shift toward more complex and realistic models is taking place, which leaves many of the theoretical and mathematical questions unexplored. Furthermore, as adaptive therapy displays better potential than existing treatment protocols, an increasing number of studies incorporate this treatment into modeling efforts. Although existing modeling work has explored and yielded useful insights on the treatment of prostate cancer, the road to clinical application is still elusive. Among the pertinent issues needed to be addressed to bridge the gap from modeling work to clinical application are (1) real-time data validation and model identification, (2) sensitivity analysis and uncertainty quantification for model prediction, and (3) optimal treatment/schedule while considering drug properties, interactions, and toxicity. To address these issues, we suggest in-depth studies on various aspects of the parameters in dynamical models such as the evolution of parameters over time. We hope this review will assist future attempts at studying prostate cancer.

Abstract Motivation Functions of cancer driver genes vary substantially across tissues and organs. Distinguishing passenger genes, oncogenes (OGs) and tumor-suppressor genes (TSGs) for each cancer type is critical for understanding tumor biology and identifying clinically actionable targets. Although many computational tools are available to predict putative cancer driver genes, resources for context-aware classifications of OGs and TSGs are limited. Results We show that the direction and magnitude of somatic selection of protein-coding mutations are significantly different for passenger genes, OGs and TSGs. Based on these patterns, we develop a new method (genes under selection in tumors) to discover OGs and TSGs in a cancer-type specific manner. Genes under selection in tumors shows a high accuracy (92%) when evaluated via strict cross-validations. Its application to 10 172 tumor exomes found known and novel cancer drivers with high tissue-specificities. In 11 out of 13 OGs shared among multiple cancer types, we found functional domains selectively engaged in different cancers, suggesting differences in disease mechanisms. Availability and implementation An R implementation of the GUST algorithm is available at https://github.com/liliulab/gust. A database with pre-computed results is available at https://liliulab.shinyapps.io/gust. Supplementary information Supplementary data are available at Bioinformatics online.

AbstractThis article reviews progress in chemopreventive drug development, especially data and concepts that are new since the 2002 AACR report on treatment and prevention of intraepithelial neoplasia. Molecular biomarker expressions involved in mechanisms of carcinogenesis and genetic progression models of intraepithelial neoplasia are discussed and analyzed for how they can inform mechanism-based, molecularly targeted drug development as well as risk stratification, cohort selection, and end-point selection for clinical trials. We outline the concept of augmenting the risk, mechanistic, and disease data from histopathologic intraepithelial neoplasia assessments with molecular biomarker data. Updates of work in 10 clinical target organ sites include new data on molecular progression, significant completed trials, new agents of interest, and promising directions for future clinical studies. This overview concludes with strategies for accelerating chemopreventive drug development, such as integrating the best science into chemopreventive strategies and regulatory policy, providing incentives for industry to accelerate preventive drugs, fostering multisector cooperation in sharing clinical samples and data, and creating public-private partnerships to foster new regulatory policies and public education.