Search

Background: Metastasis is the leading cause of cancer related deaths in breast cancer patients. Lymphovascular invasion represents one of the earliest stages of metastasis wherein the cells are introduced to a very different and distinct microenvironment. Methods: We leveraged spatial techniques developed for limited specimens in archival tissue to study patient matched cross-sectional tumor samples from different stages of breast neoplasia including normal breast, ductal carcinoma in situ (DCIS), primary invasive carcinoma (IBC), lymphovascular invasion (LVI) and regional lymph node metastasis. We selected a set of 21 patients with ER+ breast cancer to generate cross-sectional samples of each of these stages, for a total of 331 samples. The areas of LVI were identified by a combination of H&E review and immunohistochemistry for podoplanin. We performed smart-3SEQ for gene expression profiling and light pass whole genome sequencing for DNA copy number alterations. Results: We profiled the spectrum of neoplasia for transcriptome-wide gene expression. Principal component analysis of all 252 DCIS, LVI, IBC, or metastasis samples using the top 500 genes with the highest variance demonstrated that clustering was roughly based on the diagnostic stage (i.e. DCIS, LVI, IBC, or metastasis). Differential gene expression profiling identified thousands of genes increased or decreased in expression across the transitional stages with the largest change in gene expression being the transition from normal breast to DCIS, dominated by gene expression down regulation. We next performed NMF clustering on 62 samples of LVI from 18 cases and identified two patterns of gene expression which define two subgroups. Gene ontology analysis revealed that one cluster was associated with increased proliferation and metabolism, whereas the second cluster was dominated by an immune response. When we analyzed the immune and proliferative LVI subgroups separately, we found that the immune profiles in the patient matched IBC and LVI samples from the LVI Immune cluster were similar, whereas the immune profiles in the patient matched IBC and LVI samples from the Proliferative cluster were significantly different. At the LVI stage, all immune cell populations estimated by CibersortX were decreased in the Proliferative LVI cluster. These changes were validated using immunofluorescence for proliferation (Ki67), T cells (CD3) and macrophages (CD68) on the same samples. Using the LVI centroids, we built a model that could predict the same clusters in the METABRIC IBC. Kaplan-Meier analysis showed a significant difference between groups, with the Proliferative-like IBC group having a worse prognosis than the Immune-like IBC group. Conclusions: We observed a dichotomy at the LVI stage with a more proliferative cluster that may escape the immune response and an immune cluster which has a microenvironment with a similar pattern to its primary IBC. The recognition of two groups of LVI, differing in immune association and proliferation, raises the possibility that the risk of metastasis could be different in these two groups, leading to different biological pathways of progression. Citation Format: Belén Rivero-Gutiérrez, Diego Mallo, Almudena Espín-Pérez, Sujay Vennam, Chunfang Zhu, Sushama Varma, Greg Scott, Joseph Foley, E Shelley Hwang, Carlo Maley, Robert West. Characterization of the lymphovascular invasion microenvironment reveals immune response dichotomy [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 PD2-09.

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.

The invasive leading edge represents a potential gateway for tumor metastasis. The role of fibroblasts from the tumor edge in promoting cancer invasion and metastasis has not been comprehensively elucidated. We hypothesize that cross-talk between tumor and stromal cells within the tumor microenvironment results in activation of key biological pathways depending on their position in the tumor (edge vs. core). Here we highlight phenotypic differences between tumor-adjacent-fibroblasts (TAF) from the invasive edge and tumor core fibroblasts from the tumor core, established from human lung adenocarcinomas. A multiomics approach that includes genomics, proteomics, and O-glycoproteomics was used to characterize cross-talk between TAFs and cancer cells. These analyses showed that O-glycosylation, an essential posttranslational modification resulting from sugar metabolism, alters key biological pathways including the cyclin-dependent kinase 4 (CDK4) and phosphorylated retinoblastoma protein axis in the stroma and indirectly modulates proinvasive features of cancer cells. In summary, the O-glycoproteome represents a new consideration for important biological processes involved in tumor–stroma cross-talk and a potential avenue to improve the anticancer efficacy of CDK4 inhibitors. Significance: A multiomics analysis of spatially distinct fibroblasts establishes the importance of the stromal O-glycoproteome in tumor–stroma interactions at the leading edge and provides potential strategies to improve cancer treatment. See related commentary by De Wever, p. 537

Ductal Carcinoma in Situ (DCIS) is a pre-invasive lesion that accounts for nearly 20% of new breast cancer diagnoses. Of these cases, about half will progress to invasive breast cancer (IBC) within ten years. However, because diagnostic criteria for delineating low risk lesions from those likely to progress to IBC have not been identified, many patients are receiving unnecessary chemotherapy and surgery that can result in therapy-related morbidity and death. With this in mind, we used Multiplexed Ion Beam Imaging by time of flight (MIBI-TOF) and RNA-seq laser-capture microdissection (SMART-3SEQ) to construct a comprehensive spatial atlas describing the structure, function, and cellular composition of DCIS. MIBI-TOF and SMART-3SEQ were used to compare lesions from patients that later developed IBC with those from age- and history matched DCIS controls without recurrence. Using a 37-marker staining panel to interrogate 137 lesions, we identified 30 distinct cell populations of the epithelial, stromal, and immune lineages that were arranged in recurrent cellular microenvironments specific for DCIS or invasive disease. We observe a coordinated shift in the immune and stromal compartments as invasive disease arises, including an expansion of immune cell diversity and transition to reactive stromal phenotypes in synchronous DCIS + IBC, which was distinct from the macrophage-dominant microenvironment of recurrent IBC. Single-cell segmentation using a deep learning model was combined with pixel-level coexpression analysis to determine how thickness, continuity, and phenotype of ductal myoepithelium changes as tumors progress from a pre-invasive state. These data were incorporated into a comprehensive model which was subsequently used to identify a subset of features that correlate with disease-free survival following DCIS tumor resection. Taken together, these features represent important prognostic metrics that can be used to separate pre-invasive from indolent DCIS tumors, and allow for tailored therapy that improves patient outcomes and quality of life in this disease. Citation Format: Tyler Risom, Belen Rivero, Candace Liu, Alex Baranski, Siri Strand, Noah Greenwald, Erin McCaffrey, Sushama Varma, Leeat Keren, Sucheta Srivastava, Chunfang Zhu, Sujay Vennam, Shelley Hwang, Graham Colditz, Sean Bendall, Robert West, Michael Angelo. Mapping the tumor and microenvironmental evolution underlying DCIS progression through multiplexed ion beam imaging [abstract]. In: Proceedings of the AACR Virtual Special Conference on Tumor Heterogeneity: From Single Cells to Clinical Impact; 2020 Sep 17-18. Philadelphia (PA): AACR; Cancer Res 2020;80(21 Suppl):Abstract nr PR05.

Currently, only patients with HER2-positive tumors are candidates for HER2-targeted therapies. However, recent clinical observations suggest that the survival of patients with HER2-low breast cancers, who lack HER2 amplification, may benefit from adjuvant therapy that targets HER2. In this study, we explored a mechanism through which these benefits may be obtained. Prompted by the hypothesis that HER2/HER3 signaling in breast tumor-initiating cells (TIC) promotes self-renewal and survival, we obtained evidence that neuregulin 1 (NRG1) produced by TICs promotes their proliferation and self-renewal in HER2-low tumors, including in triple-negative breast tumors. Pharmacologic inhibition of EGFR, HER2, or both receptors reduced breast TIC survival and self-renewal in vitro and in vivo and increased TIC sensitivity to ionizing radiation. Through a tissue microarray analysis, we found that NRG1 expression and associated HER2 activation occurred in a subset of HER2-low breast cancers. Our results offer an explanation for why HER2 inhibition blocks the growth of HER2-low breast tumors. Moreover, they argue that dual inhibition of EGFR and HER2 may offer a useful therapeutic strategy to target TICs in these tumors. In generating a mechanistic rationale to apply HER2-targeting therapies in patients with HER2-low tumors, this work shows why these therapies could benefit a considerably larger number of patients with breast cancer than they currently reach. Cancer Res; 74(1); 341–52. ©2013 AACR.

Cancer evolution involves cycles of genomic damage, epigenetic deregulation, and increased cellular proliferation that eventually culminate in the carcinoma phenotype. Early breast neoplasias include usual ductal hyperplasia, columnar cell lesions, and flat epithelial atypia and some are thought to represent precursor stages. To elucidate their role in cancer evolution we performed comparative whole genome sequencing of early neoplasias, matched normal tissue, and carcinomas from six patients. The identified somatic mutations served as lineage markers to build trees that relate the tissue samples. On the basis of the lineage trees we inferred the order, timing, and rates of mutational events. We find that in a subset of cases, the last common ancestor of an early neoplasm and a carcinoma was hypermutated and had several aneuploidies, and that evolution further accelerated in the carcinoma lineage. In contrast to highly advanced tumors that are the focus of much of current cancer genome sequencing, the early neoplasia genomes, and the carcinomas as well, harbor a striking paucity of potentially functional somatic point mutations. The earliest significant events we could detect, those changes that occurred in common ancestors of neoplastic and tumor cells, are aneuploidies. Several aneuploidies are recurrent, suggesting that they are among the earliest genomic events that predispose breast tissue to eventual development of invasive carcinoma. Citation Information: Cancer Res 2012;72(24 Suppl):Abstract nr PD05-09.

Objective: Nasopharyngeal carcinoma (NPC) is a common Epstein-barr virus-associated epithelial malignancy in several parts of the world including Southeast Asia. While the majority of patients are treated uniformly with a combination of chemo and radiation therapy, 20% of patients experience recurrence, most commonly in the form of distant metastasis. NPC subtypes based on underlying differences in biology are unexplored and are likely responsible for the heterogeneous clinical response. As NPC biopsies are small and have significant stromal infiltrate, obtaining pure epithelial cells for genomic profiling is a challenge. We aim to overcome these limitations to identify biological subtypes and dysregulated molecular pathways in NPC. Methods: We first evaluated 217 whole exomes sequences, including 102 microdissected tumors, from three previous studies for single nucleotide variants and copy number changes, following GATK standards and using next generation sequencing copy number callers. We then applied a novel method for gene expression profiling developed in our lab to a new cohort of EBV-positive primary NPC cases from our institution. We performed laser capture microdissection, separately dissecting tumor, normal and microenvironment for each case. We applied Smart-3SEQ, a novel 3' end RNA-Seq technique which allows for the accurate quantification of transcript abundance in dissected FFPE samples comprising only a few hundred cells. Results: We achieved a per-base concordance of 80.4% between copy number profiles by SNP array and whole exome sequencing. Unsupervised clustering identified three distinct copy-number groups of NPC tumors, with a low copy-number group demonstrating an 18.7% better 5-year disease-specific survival, not attributable to stage. Apart from broad cytogenetic changes, narrow regions of amplifications (e.g. 1q21, 11q13) and deletions (e.g. 9p21, 11q22) were important for defining copy-number subtypes. Preliminary differential gene expression analysis showed that genes involved in cell cycle and cellular differentiation were significantly dysregulated in tumor cells (p < 0.001 and p = 0.03 respectively), while genes involved in cilia assembly and flagella transport were upregulated in normal cells (p < 0.001 for both). Our further analysis includes defining tumor subtypes based on gene expression, identifying key driver pathways, and correlating with EBV-latent gene expression, the immune environment, as well as clinical outcome. Conclusion: NPC tumors are biologically heterogeneous and can be classified based on their mutational and gene expression profiles. This provides an important basis for the consideration of escalation and de-escalation of therapy in selected patient groups. Citation Format: Kai Xun Joshua Tay, Chunfang Zhu, Sujay Vennam, Sushama Varma, Quynh-Thu Le, John Sunwoo, Robert West. Biological subtypes of nasopharyngeal carcinoma by genomic profiling [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2018; 2018 Apr 14-18; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2018;78(13 Suppl):Abstract nr 3411.

Sarcomas of the uterus are derived from uterine smooth muscle (leiomyosarcomas) or from endometrial stromal cells. The latter group includes three subtypes: the low-grade (LG) and high-grade (HG) endometrial stromal sarcomas (ESS), each defined by characteristic chromosomal translocations, and undifferentiated uterine sarcoma (UUS), the most aggressive form. The majority of UUS patients present with high-stage disease and even patients with stage I tumors often die within 2 years from diagnosis. Adjuvant radio- and chemotherapy do not improve the clinical outcome. No prior comprehensive molecular analysis of all three types of endometrial stromal tumors has been reported. We performed genomic analysis (whole exome sequencing and aCGH), gene expression profiling (RNA-seq and microarrays), immunohistochemistry (IHC) and FISH on 31 endometrial stromal tumors, including 17 UUS and 14 LG and HG ESS. All ESS cases carried characteristic JAZF1, PHF1, MBTD1 or YWHAE rearrangements. Selected findings were validated by IHC, qRT-PCR and Sanger sequencing. Genomic analysis revealed that UUS are characterized by complex chromosomal aberrations. The most frequent somatic mutations (in TP53, KRAS, FBXW7, PIK3CA and ERBB3) and copy number changes (e.g. gains of 19q, 8q11.1-q24.3 and 3q26.2-q29) in UUS resemble those seen in uterine carcinosarcomas and carcinomas. UUS over-expressed numerous genes encoding M2 macrophage-specific markers, including CD163, CCL18, mannose receptor, stabilin-1, and macrophage galactose-type C-type lectin 2. Immunohistochemistry for CD163 and CD68 confirmed high tumor-associated macrophages (TAM) counts in UUS tissues compared to the ESS. In UUS, we also observed significantly higher mRNA expression for genes implicated in angiogenesis (CD34, MMP9), immunosuppression and tumor invasion (e.g. CCR2, IL10RA, IRAK3, CCL13, CXCL9, CXCL10). TAMs are associated with progression, resistance to cytotoxic therapy and metastatic spread in most human tumor types. In animal models, TAMs were shown to induce increased angiogenesis. CSF1 is a major regulator of macrophage function and is implicated in these tumor-promoting effects. Clinical trials have shown an effect of inhibitors of the CSF1 pathway in tenosynovial giant cell tumors and our data indicate that these inhibitors may be considered for the treatment of UUS. Our study also showed significantly elevated expression in UUS of two additional therapeutic targets involving macrophage pathways: CCR2 (C-C chemokine receptor type 2) that can be targeted by monoclonal antibodies, and BTK (Bruton's tyrosine kinase) that can be targeted by ibrutinib. In conclusion, our findings demonstrate abundant presence of TAMs and over-expression of TAM-associated markers in UUS compared to the other tumors derived from endometrial stroma. Moreover, our results indicate novel potential therapeutic targets for UUS patients management. Citation Format: Joanna Przybyl, Magdalena Kowalewska, Anna Quattrone, Barbara Dewaele, Vanessa Vanspauwen, Sushama Varma, Robert T. Sweeney, Michal Swierniak, Elwira Bakula-Zalewska, Janusz A. Siedlecki, Mariusz Bidzinski, Jan Cools, Matt van de Rijn, Maria Debiec-Rychter. Tumor associated macrophages in undifferentiated uterine sarcoma: association with angiogenesis and therapeutic implications. [abstract]. In: Proceedings of the 107th Annual Meeting of the American Association for Cancer Research; 2016 Apr 16-20; New Orleans, LA. Philadelphia (PA): AACR; Cancer Res 2016;76(14 Suppl):Abstract nr 3182.

Rare cancer types are not only understudied, but are typically represented by formalin-fixed paraffin-embedded (FFPE) (rather than freshly-frozen) specimens that are suboptimal for genomic analysis. Ameloblastoma is one such rare tumor type, thought to arise from ameloblasts, the cells that deposit enamel during tooth development. Though typically benign, ameloblastomas are locally destructive to the jaw and face, and new non-surgical interventions are needed. To discover novel driver mutations and therapeutic targets, we optimized methods and performed whole-transcriptome sequencing and/or targeted exon sequencing (TruSeq Cancer Panel) of 8 FFPE cases. Identified mutations were verified, and then evaluated on a larger, independent set of 21 FFPE cases by PCR and Sanger sequencing. From the analysis, we identified recurrent somatic mutations in three key developmental or signaling pathways, including Hedgehog, fibroblast growth factor, and MAP kinase pathways. Functional interrogation of a novel Hedgehog pathway mutation confirmed increased basal pathway activity, and defined the response profile to various pharmacologic Hedgehog inhibitors. Together, our results define new ameloblastoma drivers and nominate new molecularly-directed therapies for this rare but disfiguring disease. More generally, our findings validate a robust approach for discovering driver mutations in rare cancers. Citation Format: Andrew C. McClary, Robert T. Sweeney, Jewison Biscocho, Benjamin R. Myers, Lila Neahring, Kevin A. Kwei, Kunbin Qu, Xue Gong, Tony Ng, Carol D. Jones, Sushama Varma, Justin I. Odegaard, Brian Rubin, Megan L. Troxell, Robert J. Pelham, James L. Zehnder, Philip A. Beachy, Jonathan R. Pollack, Robert B. West. Ameloblastoma driver mutations revealed by next-generation sequencing of formalin-fixed paraffin-embedded specimens. [abstract]. In: Proceedings of the 105th Annual Meeting of the American Association for Cancer Research; 2014 Apr 5-9; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2014;74(19 Suppl):Abstract nr 3436. doi:10.1158/1538-7445.AM2014-3436