Mutation order in acute myeloid leukemia identifies uncommon patterns of evolution and illuminates phenotypic heterogeneity

Acute myeloid leukemia (AML) has a poor prognosis and a heterogeneous mutation landscape. Although common mutations are well-studied, little research has characterized how the sequence of mutations relates to clinical features. Using published, single-cell DNA sequencing data from three institutions, we compared clonal evolution patterns in AML to patient characteristics, disease phenotype, and outcomes. Mutation trees, which represent the order of select mutations, were created for 207 patients from targeted panel sequencing data using 1 639 162 cells, 823 mutations, and 275 samples. In 224 distinct orderings of mutated genes, mutations related to DNA methylation typically preceded those related to cell signaling, but signaling-first cases did occur, and had higher peripheral cell counts, increased signaling mutation homozygosity, and younger patient age. Serial sample analysis suggested that NPM1 and DNA methylation mutations provide an advantage to signaling mutations in AML. Interestingly, WT1 mutation evolution shared features with signaling mutations, such as WT1-early being proliferative and occurring in younger individuals, trends that remained in multivariable regression. Some mutation orderings had a worse prognosis, but this was mediated by unfavorable mutations, not mutation order. These findings add a dimension to the mutation landscape of AML, identifying uncommon patterns of leukemogenesis and shedding light on heterogeneous phenotypes.

Loss of p53-DREAM-mediated repression of cell cycle genes as a driver of lymph node metastasis in head and neck cancer

BACKGROUND

The prognosis for patients with head and neck cancer (HNC) is poor and has improved little in recent decades, partially due to lack of therapeutic options. To identify effective therapeutic targets, we sought to identify molecular pathways that drive metastasis and HNC progression, through large-scale systematic analyses of transcriptomic data.

METHODS

We performed meta-analysis across 29 gene expression studies including 2074 primary HNC biopsies to identify genes and transcriptional pathways associated with survival and lymph node metastasis (LNM). To understand the biological roles of these genes in HNC, we identified their associated cancer pathways, as well as the cell types that express them within HNC tumor microenvironments, by integrating single-cell RNA-seq and bulk RNA-seq from sorted cell populations.

RESULTS

Patient survival-associated genes were heterogenous and included drivers of diverse tumor biological processes: these included tumor-intrinsic processes such as epithelial dedifferentiation and epithelial to mesenchymal transition, as well as tumor microenvironmental factors such as T cell-mediated immunity and cancer-associated fibroblast activity. Unexpectedly, LNM-associated genes were almost universally associated with epithelial dedifferentiation within malignant cells. Genes negatively associated with LNM consisted of regulators of squamous epithelial differentiation that are expressed within well-differentiated malignant cells, while those positively associated with LNM represented cell cycle regulators that are normally repressed by the p53-DREAM pathway. These pro-LNM genes are overexpressed in proliferating malignant cells of TP53 mutated and HPV + ve HNCs and are strongly associated with stemness, suggesting that they represent markers of pre-metastatic cancer stem-like cells. LNM-associated genes are deregulated in high-grade oral precancerous lesions, and deregulated further in primary HNCs with advancing tumor grade and deregulated further still in lymph node metastases.

CONCLUSIONS

In HNC, patient survival is affected by multiple biological processes and is strongly influenced by the tumor immune and stromal microenvironments. In contrast, LNM appears to be driven primarily by malignant cell plasticity, characterized by epithelial dedifferentiation coupled with EMT-independent proliferation and stemness. Our findings postulate that LNM is initially caused by loss of p53-DREAM-mediated repression of cell cycle genes during early tumorigenesis.

Mutation order in acute myeloid leukemia identifies uncommon patterns of evolution and illuminates phenotypic heterogeneity

Acute myeloid leukemia (AML) has a poor prognosis and a heterogeneous mutation landscape. Although common mutations are well-studied, little research has characterized how the sequence of mutations relates to clinical features. Using published, single-cell DNA sequencing data from three institutions, we compared clonal evolution patterns in AML to patient characteristics, disease phenotype, and outcomes. Mutation trees, which represent the order of select mutations, were created for 207 patients from targeted panel sequencing data using 1 639 162 cells, 823 mutations, and 275 samples. In 224 distinct orderings of mutated genes, mutations related to DNA methylation typically preceded those related to cell signaling, but signaling-first cases did occur, and had higher peripheral cell counts, increased signaling mutation homozygosity, and younger patient age. Serial sample analysis suggested that NPM1 and DNA methylation mutations provide an advantage to signaling mutations in AML. Interestingly, WT1 mutation evolution shared features with signaling mutations, such as WT1-early being proliferative and occurring in younger individuals, trends that remained in multivariable regression. Some mutation orderings had a worse prognosis, but this was mediated by unfavorable mutations, not mutation order. These findings add a dimension to the mutation landscape of AML, identifying uncommon patterns of leukemogenesis and shedding light on heterogenous phenotypes.

High-resolution alignment of single-cell and spatial transcriptomes with CytoSPACE

Recent studies have emphasized the importance of single-cell spatial biology, yet available assays for spatial transcriptomics have limited gene recovery or low spatial resolution. Here we introduce CytoSPACE, an optimization method for mapping individual cells from a single-cell RNA sequencing atlas to spatial expression profiles. Across diverse platforms and tissue types, we show that CytoSPACE outperforms previous methods with respect to noise tolerance and accuracy, enabling tissue cartography at single-cell resolution.

AML/T cell interactomics uncover correlates of patient outcomes and the key role of ICAM1 in T cell killing of AML

T cells are important for the control of acute myeloid leukemia (AML), a common and often deadly malignancy. We observed that some AML patient samples are resistant to killing by human engineered cytotoxic CD4 + T cells. Single-cell RNA-seq of primary AML samples and CD4 + T cells before and after their interaction uncovered transcriptional programs that correlate with AML sensitivity or resistance to CD4 + T cell killing. Resistance-associated AML programs were enriched in AML patients with poor survival, and killing-resistant AML cells did not engage T cells in vitro . Killing-sensitive AML potently activated T cells before being killed, and upregulated ICAM1 , a key component of the immune synapse with T cells. Without ICAM1, killing-sensitive AML became resistant to killing to primary ex vivo -isolated CD8 + T cells in vitro , and engineered CD4 + T cells in vitro and in vivo . Thus, ICAM1 on AML acts as an immune trigger, allowing T cell killing, and could affect AML patient survival in vivo .

SIGNIFICANCE

AML is a common leukemia with sub-optimal outcomes. We show that AML transcriptional programs correlate with susceptibility to T cell killing. Killing resistance-associated AML programs are enriched in patients with poor survival. Killing-sensitive, but not resistant AML activate T cells and upregulate ICAM1 that binds to LFA-1 on T cells, allowing immune synapse formation which is critical for AML elimination.

GRAPHICAL ABSTRACT

Translatome analysis reveals microglia and astrocytes to be distinct regulators of inflammation in the hyperacute and acute phases after stroke

Neuroinflammation is a hallmark of ischemic stroke, which is a leading cause of death and long-term disability. Understanding the exact cellular signaling pathways that initiate and propagate neuroinflammation after stroke will be critical for developing immunomodulatory stroke therapies. In particular, the precise mechanisms of inflammatory signaling in the clinically relevant hyperacute period, hours after stroke, have not been elucidated. We used the RiboTag technique to obtain microglia and astrocyte-derived mRNA transcripts in a hyperacute (4 h) and acute (3 days) period after stroke, as these two cell types are key modulators of acute neuroinflammation. Microglia initiated a rapid response to stroke at 4 h by adopting an inflammatory profile associated with the recruitment of immune cells. The hyperacute astrocyte profile was marked by stress response genes and transcription factors, such as Fos and Jun, involved in pro-inflammatory pathways such as TNF-α. By 3 days, microglia shift to a proliferative state and astrocytes strengthen their inflammatory response. The astrocyte pro-inflammatory response at 3 days is partially driven by the upregulation of the transcription factors C/EBPβ, Spi1, and Rel, which comprise 25% of upregulated transcription factor-target interactions. Surprisingly, few sex differences across all groups were observed. Expression and log₂ fold data for all sequenced genes are available on a user-friendly website for researchers to examine gene changes and generate hypotheses for stroke targets. Taken together, our data comprehensively describe the microglia and astrocyte-specific translatome response in the hyperacute and acute period after stroke and identify pathways critical for initiating neuroinflammation.

Abstract 3469: Lymph node colonization promotes distant tumor metastasis through the induction of tumor-specific immune tolerance

The majority of cancer-associated deaths result from distant organ metastasis, yet the mechanisms that enable this process remain poorly understood. For most solid tumors, colonization of regional or distant lymph nodes (LNs) typically precedes the formation of distant organ metastases, yet it remains unclear whether LN metastasis plays a functional role in disease progression. LNs are major sites of anti-tumor lymphocyte education, including in the context of immunotherapy, yet LN metastasis frequently correlates with further disease progression. Here, we find that LN metastasis represents a critical step in tumor progression through the capacity of such metastases to induce tumor-specific immune tolerance in a manner that promotes further dissemination of tumors to distant organs. Using an in vivo passaging approach of a non-metastatic syngeneic melanoma, we generated 300 unique cell lines exhibiting varying degrees of LN metastatic capacity. We show that the presence of these LN metastases enables distant organ seeding of metastases in a manner that the parental tumor cannot, and this effect is eliminated in mice lacking an adaptive immune response. Furthermore, this promotion of distant seeding by LN metastases is tumor specific. Using flow cytometry and single-cell sequencing to perform comprehensive immune profiling, we identify multiple cellular mediators of tolerance. In particular, we find that LN metastases have the capacity to both resist NK cell cytotoxicity and induce regulatory T cells (Tregs). Furthermore, depletion of NK cells in vivo enables non-metastatic tumors to disseminate to LNs, and ablation of Tregs using FoxP3-DTR mice eliminates the occurrence of lymphatic metastases. Adoptive transfer of Tregs from the LNs of mice bearing LN metastasis to naïve mice facilitates metastasis in a manner that Tregs from mice without LN metastases cannot, and we find that these Tregs are induced in an antigen-specific manner. Whole exome sequencing revealed that neither the metastatic proclivity nor immunosuppression evolve through the acquisition of driver mutations, loss of neoantigens, loss of MHC class I presentation, or decreases in melanoma antigen expression. Rather, by RNA-seq and ATAC-seq, we show that a conserved interferon signaling axis is upregulated in LN metastases and is rendered stable through epigenetic reprogramming of chromatin accessibility resulting from chronic exposure to interferons in vivo. Furthermore, using CRISPR/Cas9, we find that these pathways are required for LN metastatic seeding, and validate their conserved significance in additional mouse models of pancreatic ductal adenocarcinoma and head and neck squamous cell carcinoma and humans with LN metastatic disease. Together, these findings demonstrate a critical role for LN metastasis in promoting tumor-specific immunosuppression.

Citation Format: Nathan E. Reticker-Flynn, Weiruo Zhang, Julia A. Belk, Pamela A. Basto, Andrew J. Gentles, John B. Sunwoo, Ansuman T. Satpathy, Sylvia K. Plevritis, Edgar G. Engleman. Lymph node colonization promotes distant tumor metastasis through the induction of tumor-specific immune tolerance [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 3469.

Multimodal data fusion for cancer biomarker discovery with deep learning

Technological advances now make it possible to study a patient from multiple angles with high-dimensional, high-throughput multi-scale biomedical data. In oncology, massive amounts of data are being generated ranging from molecular, histopathology, radiology to clinical records. The introduction of deep learning has significantly advanced the analysis of biomedical data. However, most approaches focus on single data modalities leading to slow progress in methods to integrate complementary data types. Development of effective multimodal fusion approaches is becoming increasingly important as a single modality might not be consistent and sufficient to capture the heterogeneity of complex diseases to tailor medical care and improve personalised medicine. Many initiatives now focus on integrating these disparate modalities to unravel the biological processes involved in multifactorial diseases such as cancer. However, many obstacles remain, including lack of usable data as well as methods for clinical validation and interpretation. Here, we cover these current challenges and reflect on opportunities through deep learning to tackle data sparsity and scarcity, multimodal interpretability, and standardisation of datasets.

Translatome analysis reveals microglia and astrocytes to be distinct regulators of inflammation in the hyperacute and acute phases after stroke

Neuroinflammation is a hallmark of ischemic stroke, which is a leading cause of death and long-term disability. Understanding the exact cellular signaling pathways that initiate and propagate neuroinflammation after stroke will be critical for developing immunomodulatory stroke therapies. In particular, the precise mechanisms of inflammatory signaling in the clinically relevant hyperacute period, hours after stroke, have not been elucidated. We used the RiboTag technique to obtain astrocyte and microglia-derived mRNA transcripts in a hyperacute (4 hours) and acute (3 days) period after stroke, as these two cell types are key modulators of acute neuroinflammation. Microglia initiated a rapid response to stroke at 4 hours by adopting an inflammatory profile associated with the recruitment of immune cells. The hyperacute astrocyte profile was marked by stress response genes and transcription factors, such as Fos and Jun , involved in pro-inflammatory pathways such as TNF-α. By 3 days, microglia shift to a proliferative state and astrocytes strengthen their inflammatory response. The astrocyte pro-inflammatory response at 3 days is partially driven by the upregulation of the transcription factors C/EBPβ, Spi1 , and Rel , which comprise 25% of upregulated transcription factor-target interactions. Surprisingly, few sex differences across all groups were observed. Expression and log ₂ fold data for all sequenced genes are available on a user-friendly website for researchers to examine gene changes and generate hypotheses for stroke targets. Taken together our data comprehensively describe the astrocyte and microglia-specific translatome response in the hyperacute and acute period after stroke and identify pathways critical for initiating neuroinflammation.

Single-cell genomics in AML: extending the frontiers of AML research

The era of genomic medicine has allowed acute myeloid leukemia (AML) researchers to improve disease characterization, optimize risk-stratification systems, and develop new treatments. Although there has been significant progress, AML remains a lethal cancer because of its remarkably complex and plastic cellular architecture. This degree of heterogeneity continues to pose a major challenge, because it limits the ability to identify and therefore eradicate the cells responsible for leukemogenesis and treatment failure. In recent years, the field of single-cell genomics has led to unprecedented strides in the ability to characterize cellular heterogeneity, and it holds promise for the study of AML. In this review, we highlight advancements in single-cell technologies, outline important shortcomings in our understanding of AML biology and clinical management, and discuss how single-cell genomics can address these shortcomings as well as provide unique opportunities in basic and translational AML research.

Abstract PR013: Lymph node colonization promotes distant tumor metastasis through the induction of tumor-specific immune tolerance

The majority of cancer-associated deaths result from distant organ metastasis, yet the mechanisms that enable this process remain poorly understood. For most solid tumors, colonization of regional or distant lymph nodes (LNs) typically precedes the formation of distant organ metastases, yet it remains unclear whether LN metastasis plays a functional role in disease progression. LNs are major sites of anti-tumor lymphocyte education, including in the context of immunotherapy, yet LN metastasis frequently correlates with further disease progression. Here, we find that LN metastasis represents a critical step in tumor progression through the capacity of such metastases to induce tumor-specific immune tolerance in a manner that promotes further dissemination of tumors to distant organs. Using an in vivo passaging approach of a non-metastatic syngeneic melanoma, we generated 300 unique cell lines exhibiting varying degrees of LN metastatic capacity. We show that the presence of these LN metastases enables distant organ seeding of metastases in a manner that the parental tumor cannot, and this effect is eliminated in mice lacking an adaptive immune response. Furthermore, this promotion of distant seeding by LN metastases is tumor specific. Using flow cytometry and single-cell sequencing to perform comprehensive immune profiling, we identify multiple cellular mediators of tolerance. In particular, we find that LN metastases have the capacity to both resist NK cell cytotoxicity and induce regulatory T cells (Tregs). Furthermore, depletion of NK cells in vivo enables non-metastatic tumors to disseminate to LNs, and ablation of Tregs using FoxP3-DTR mice eliminates the occurrence of lymphatic metastases. Adoptive transfer of Tregs from the LNs of mice bearing LN metastasis to naïve mice facilitates metastasis in a manner that Tregs from mice without LN metastases cannot, and we find that these Tregs are induced in an antigen-specific manner. Whole exome sequencing revealed that neither the metastatic proclivity nor immunosuppression evolve through the acquisition of driver mutations, loss of neoantigens, loss of MHC class I presentation, or decreases in melanoma antigen expression. Rather, by RNA-seq and ATAC-seq, we show that a conserved interferon signaling axis is upregulated in LN metastases and is rendered stable through epigenetic reprogramming of chromatin accessibility resulting from chronic exposure to interferons in vivo. Furthermore, using CRISPR/Cas9, we find that these pathways are required for LN metastatic seeding, and validate their conserved significance in additional mouse models of pancreatic ductal adenocarcinoma and head and neck squamous cell carcinoma and humans with LN metastatic disease. Together, these findings demonstrate a critical role for LN metastasis in promoting tumor-specific immunosuppression.

Citation Format: Nathan E. Reticker-Flynn, Weiruo Zhang, Julia A. Belk, Pamela A. Basto, Andrew J. Gentles, John B. Sunwoo, Ansuman K. Satpathy, Sylvia K. Plevritis, Edgar G. Engleman. Lymph node colonization promotes distant tumor metastasis through the induction of tumor-specific immune tolerance [abstract]. In: Proceedings of the AACR Special Conference: Cancer Metastasis; 2022 Nov 14-17; Portland, OR. Philadelphia (PA): AACR; Cancer Res 2022;83(2 Suppl_2):Abstract nr PR013.

Profiling Cellular Ecosystems at Single-Cell Resolution and at Scale with EcoTyper

Tissues are composed of diverse cell types and cellular states that organize into distinct ecosystems with specialized functions. EcoTyper is a collection of machine learning tools for the large-scale delineation of cellular ecosystems and their constituent cell states from bulk, single-cell, and spatially resolved gene expression data. In this chapter, we provide a primer on EcoTyper and demonstrate its use for the discovery and recovery of cell states and ecosystems from healthy and diseased tissue specimens.

Peripheral blood DNA methylation profiles predict future development of B-cell Non-Hodgkin Lymphoma

Lack of accurate methods for early lymphoma detection limits the ability to cure patients. Since patients with Non-Hodgkin lymphomas (NHL) who present with advanced disease have worse outcomes, accurate and sensitive methods for early detection are needed to improve patient care. We developed a DNA methylation-based prediction tool for NHL, based on blood samples collected prospectively from 278 apparently healthy patients who were followed for up to 16 years to monitor for NHL development. A predictive score was developed using machine learning methods in a robust training/validation framework. Our predictive score incorporates CpG DNA methylation at 135 genomic positions, with higher scores predicting higher risk. It was 85% and 78% accurate for identifying patients at risk of developing future NHL, in patients with high or low epigenetic mitotic clock respectively, in a validation cohort. It was also sensitive at detecting active NHL (96.3% accuracy) and healthy status (95.6% accuracy) in additional independent cohorts. Scores optimized for specific NHL subtypes showed significant but lower accuracy for predicting other subtypes. Our score incorporates hyper-methylation of Polycomb and HOX genes, which have roles in NHL development, as well as PAX5 - a master transcriptional regulator of B-cell fate. Subjects with higher risk scores showed higher regulatory T-cells, memory B-cells, but lower naïve T helper lymphocytes fractions in the blood. Future prospective studies will be required to confirm the utility of our signature for managing patients who are at high risk for developing future NHL.

NSD1 mutations deregulate transcription and DNA methylation of bivalent developmental genes in Sotos syndrome

Sotos syndrome (SS), the most common overgrowth with intellectual disability (OGID) disorder, is caused by inactivating germline mutations of NSD1, which encodes a histone H3 lysine 36 methyltransferase. To understand how NSD1 inactivation deregulates transcription and DNA methylation (DNAm), and to explore how these abnormalities affect human development, we profiled transcription and DNAm in SS patients and healthy control individuals. We identified a transcriptional signature that distinguishes individuals with SS from controls and was also deregulated in NSD1-mutated cancers. Most abnormally expressed genes displayed reduced expression in SS; these downregulated genes consisted mostly of bivalent genes and were enriched for regulators of development and neural synapse function. DNA hypomethylation was strongly enriched within promoters of transcriptionally deregulated genes: overexpressed genes displayed hypomethylation at their transcription start sites while underexpressed genes featured hypomethylation at polycomb binding sites within their promoter CpG island shores. SS patients featured accelerated molecular aging at the levels of both transcription and DNAm. Overall, these findings indicate that NSD1-deposited H3K36 methylation regulates transcription by directing promoter DNA methylation, partially by repressing polycomb repressive complex 2 (PRC2) activity. These findings could explain the phenotypic similarity of SS to OGID disorders that are caused by mutations in PRC2 complex-encoding genes.

Abstract CT142: GD2.Ox40.CD28.z CAR T cell trial in neuroblastoma and osteosarcoma

Background: Chimeric antigen receptor T cells (CARTs) hold promising therapeutic potential for refractory tumors. GD2 is a tumor antigen expressed on neuroblastoma and osteosarcoma. In initial studies, T cells expressing 1st generation GD2-CARTs were shown to be safe and mediated modest antitumor activity in some patients with refractory neuroblastoma.

Methods: We developed a 3rd generation GD2-CART (GD2-CAR.OX40.28.z.ICD9) and conducted a phase I trial (NCT02107963) to determine the feasibility of producing and safety of administering escalating doses of GD2-CARTs in children and young adults with GD2+ solid tumors, including neuroblastoma and osteosarcoma, following cyclophosphamide-based lymphodepletion. Patient samples were evaluated for cytokine profile kinetics, immune phenotype analysis with mass cytometry (CyTOF), transcriptomic evaluation with RNA-sequencing (RNA-seq), epigenetic determination with Assay for Transposase-Accessible Chromatin using sequencing (ATAC-seq), and functional studies with flow cytometry.

Results: 15 patients aged 8-28 years were enrolled on four dose levels, of which 13 patients were infused. No dose-limiting toxicities were observed and administration of up to 1x107 GD2-CART/kg was feasible and safe for children and young adults with neuroblastoma and osteosarcoma. 15.4% (2/13) of patients experienced grade-1 cytokine release syndrome (CRS) and no neurological toxicity was observed. We measured the expansion and persistence of adoptively transferred GD2-CARTs in the peripheral blood. GD2-CARTs expanded in all patients receiving treatment, half of whom had expansion similar to that seen in clinically active CD19 and CD22 CARTs, but the GD2-CARTs had limited persistence. At Day 28 following GD2-CART infusion, 23.1% (3/13) of evaluable patients had progressive disease and 76.9% (10/13) had stable disease (SD). 3/10 SD patients remained stable at 60 days post-GD2-CART, but all patients eventually progressed. Since a major barrier to CART efficacy is inadequate CART expansion, we comprehensively evaluated for phenotypic, transcriptomic, and epigenetic immune profiles of patient apheresis, CART product, and post-treatment peripheral blood samples to identify determinants of CART expansion. GD2-CART expansion is significantly correlated with several T cell markers, and a larger baseline naïve and central memory T cell pool. Unique myeloid populations are associated with CART expansion. ATACseq identifies epigenetic differences in pre-treatment apheresis that may predict good CAR expansion in patients.

Conclusions: GD2-CART therapy following cyclophosphamide conditioning was well tolerated at all four dose levels in pediatric and young adult patients with neuroblastoma and osteosarcoma. Subsequent multi-dimensional analyses suggest key mechanisms underlying CART biology and function and highlight the potential of defining and applying molecular signatures in apheresis and CART product as biomarkers and prognostic indicators of CART expansion, with promise for advancing immunotherapies for solid tumor patients in the future.

Citation Format: Sneha Ramakrishna, Sabina Kaczanowska, Tara Murty, Cristina F. Contreras, Melinda Merchant, John Glod, Norma Gutierrez, Ahmad Alimadadi, David Stroncek, Steven Highfill, Caroline Duault, Priyanka B. Subrahmanyam, Tyson Holmes, Warren Reynolds, Reema Baskar, Antoine Barge, Hayley Lyon, Radim Moravec, Srinika Ranasinghe, Joyce Yu, Roshni Biswas, Samuel Pollack, Stephen Van Nostrand, James Lindsay, Mina Pichavant, Bita Sahaf, Sean C. Bendall, Andrew J. Gentles, Holden Maecker, Catherine C. Hedrick, Crystal Mackall, Rosandra Kaplan. GD2.Ox40.CD28.z CAR T cell trial in neuroblastoma and osteosarcoma [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 CT142.

Lymph node colonization induces tumor-immune tolerance to promote distant metastasis

For many solid malignancies, lymph node (LN) involvement represents a harbinger of distant metastatic disease and, therefore, an important prognostic factor. Beyond its utility as a biomarker, whether and how LN metastasis plays an active role in shaping distant metastasis remains an open question. Here, we develop a syngeneic melanoma mouse model of LN metastasis to investigate how tumors spread to LNs and whether LN colonization influences metastasis to distant tissues. We show that an epigenetically instilled tumor-intrinsic interferon response program confers enhanced LN metastatic potential by enabling the evasion of NK cells and promoting LN colonization. LN metastases resist T cell-mediated cytotoxicity, induce antigen-specific regulatory T cells, and generate tumor-specific immune tolerance that subsequently facilitates distant tumor colonization. These effects extend to human cancers and other murine cancer models, implicating a conserved systemic mechanism by which malignancies spread to distant organs.

Abstract PR05: Lymph node colonization promotes distant tumor metastasis through the induction of tumor-specific immune tolerance

The majority of cancer-associated deaths result from distant organ metastasis, yet the mechanisms that enable this process remain poorly understood. For most solid tumors, colonization of regional or distant lymph nodes (LNs) typically precedes the formation of distant organ metastases, yet it remains unclear whether LN metastasis plays a functional role in disease progression. LNs are major sites of anti-tumor lymphocyte education, including in the context of immunotherapy, yet LN metastasis frequently correlates with further disease progression. Here, we find that LN metastasis represents a critical step in tumor progression through the capacity of such metastases to induce tumor-specific immune tolerance in a manner that promotes further dissemination of tumors to distant organs. Using an in vivo passaging approach of a non-metastatic syngeneic melanoma, we generated 300 unique cell lines exhibiting varying degrees of LN metastatic capacity. We show that the presence of these LN metastases enables distant organ seeding of metastases in a manner that the parental tumor cannot, and this effect is eliminated in mice lacking an adaptive immune response. Furthermore, this promotion of distant seeding by LN metastases is tumor specific. Using flow cytometry and single-cell sequencing to perform organism-wide immune profiling, we identify multiple cellular mediators of tolerance. In particular, we find that LN metastases have the capacity to both resist NK cell cytotoxicity and induce regulatory T cells (Tregs) in vitro. Furthermore, depletion of NK cells in vivo enables non-metastatic tumors to disseminate to LNs, and ablation of Tregs using FoxP3-DTR mice eliminates the occurrence of lymphatic metastases. Adoptive transfer of Tregs from the LNs of mice bearing LN metastasis to naïve mice facilitates metastasis in a manner that Tregs from mice without LN metastases cannot, and we find that these Tregs are induced in an antigen-specific manner. Using genetic mouse models and photoconvertible tracking technologies, we show that Tregs induced within involved LNs preferentially traffic to distant sites compared to other CD4 populations. Through the use of whole exome sequencing, we show that neither the metastatic proclivity nor immunosuppression evolve through the acquisition of driver mutations, loss of neoantigens, loss of MHC class I presentation, or decreases in melanoma antigen expression. Rather, by RNA-seq and ATAC-seq, we show that a conserved interferon signaling axis is upregulated in LN metastases and is rendered stable through epigenetic regulation of chromatin accessibility. Furthermore, using CRISPR/Cas9, we find that these pathways are required for LN metastatic seeding, and validate their conserved significance in additional mouse models of pancreatic ductal adenocarcinoma and head and neck squamous cell carcinoma (HNSCC), along with RNA-seq analysis of malignant populations sorted from HNSCC patients. Together, these findings demonstrate a critical role for LN metastasis in promoting tumor-specific immunosuppression.

This abstract is also being presented as Poster P020.

Citation Format: Nathan E. Reticker-Flynn, Weiruo Zhang, Julia A. Belk, Andrew J. Gentles, Ansuman Satpathy, Sylvia K. Plevritis, Edgar G. Engleman. Lymph node colonization promotes distant tumor metastasis through the induction of tumor-specific immune tolerance [abstract]. In: Abstracts: AACR Virtual Special Conference: Tumor Immunology and Immunotherapy; 2021 Oct 5-6. Philadelphia (PA): AACR; Cancer Immunol Res 2022;10(1 Suppl):Abstract nr PR05.

Clonal architecture predicts clinical outcomes and drug sensitivity in acute myeloid leukemia

The impact of clonal heterogeneity on disease behavior or drug response in acute myeloid leukemia remains poorly understood. Using a cohort of 2,829 patients, we identify features of clonality associated with clinical features and drug sensitivities. High variant allele frequency for 7 mutations (including NRAS and TET2) associate with dismal prognosis; elevated GATA2 variant allele frequency correlates with better outcomes. Clinical features such as white blood cell count and blast percentage correlate with the subclonal abundance of mutations such as TP53 and IDH1. Furthermore, patients with cohesin mutations occurring before NPM1, or transcription factor mutations occurring before splicing factor mutations, show shorter survival. Surprisingly, a branched pattern of clonal evolution is associated with superior clinical outcomes. Finally, several mutations (including NRAS and IDH1) predict drug sensitivity based on their subclonal abundance. Together, these results demonstrate the importance of assessing clonal heterogeneity with implications for prognosis and actionable biomarkers for therapy.

Atlas of clinically distinct cell states and ecosystems across human solid tumors

Determining how cells vary with their local signaling environment and organize into distinct cellular communities is critical for understanding processes as diverse as development, aging, and cancer. Here we introduce EcoTyper, a machine learning framework for large-scale identification and validation of cell states and multicellular communities from bulk, single-cell, and spatially resolved gene expression data. When applied to 12 major cell lineages across 16 types of human carcinoma, EcoTyper identified 69 transcriptionally defined cell states. Most states were specific to neoplastic tissue, ubiquitous across tumor types, and significantly prognostic. By analyzing cell-state co-occurrence patterns, we discovered ten clinically distinct multicellular communities with unexpectedly strong conservation, including three with myeloid and stromal elements linked to adverse survival, one enriched in normal tissue, and two associated with early cancer development. This study elucidates fundamental units of cellular organization in human carcinoma and provides a framework for large-scale profiling of cellular ecosystems in any tissue.

The landscape of tumor cell states and ecosystems in diffuse large B cell lymphoma

Biological heterogeneity in diffuse large B cell lymphoma (DLBCL) is partly driven by cell-of-origin subtypes and associated genomic lesions, but also by diverse cell types and cell states in the tumor microenvironment (TME). However, dissecting these cell states and their clinical relevance at scale remains challenging. Here, we implemented EcoTyper, a machine-learning framework integrating transcriptome deconvolution and single-cell RNA sequencing, to characterize clinically relevant DLBCL cell states and ecosystems. Using this approach, we identified five cell states of malignant B cells that vary in prognostic associations and differentiation status. We also identified striking variation in cell states for 12 other lineages comprising the TME and forming cell state interactions in stereotyped ecosystems. While cell-of-origin subtypes have distinct TME composition, DLBCL ecosystems capture clinical heterogeneity within existing subtypes and extend beyond cell-of-origin and genotypic classes. These results resolve the DLBCL microenvironment at systems-level resolution and identify opportunities for therapeutic targeting (https://ecotyper.stanford.edu/lymphoma).

High-grade serous ovarian tumor cells modulate NK cell function to create an immune-tolerant microenvironment

Tubo-ovarian high-grade serous carcinoma (HGSC) is unresponsive to immune checkpoint blockade despite significant frequencies of exhausted T cells. Here we apply mass cytometry and uncover decidual-like natural killer (dl-NK) cell subpopulations (CD56+CD9+CXCR3+KIR+CD3-CD16-) in newly diagnosed HGSC samples that correlate with both tumor and transitioning epithelial-mesenchymal cell abundance. We show different combinatorial expression patterns of ligands for activating and inhibitory NK receptors within three HGSC tumor compartments: epithelial (E), transitioning epithelial-mesenchymal (EV), and mesenchymal (vimentin expressing [V]), with a more inhibitory ligand phenotype in V cells. In cocultures, NK-92 natural killer cells acquire CD9 from HGSC tumor cells by trogocytosis, resulting in reduced anti-tumor cytokine production and cytotoxicity. Cytotoxicity in these cocultures is restored with a CD9-blocking antibody or CD9 CRISPR knockout, thereby identifying mechanisms of immune suppression in HGSC. CD9 is widely expressed in HGSC tumors and so represents an important new therapeutic target with immediate relevance for NK immunotherapy.

Transient rest restores functionality in exhausted CAR-T cells through epigenetic remodeling

T cell exhaustion limits immune responses against cancer and is a major cause of resistance to chimeric antigen receptor (CAR)-T cell therapeutics. Using murine xenograft models and an in vitro model wherein tonic CAR signaling induces hallmark features of exhaustion, we tested the effect of transient cessation of receptor signaling, or rest, on the development and maintenance of exhaustion. Induction of rest through enforced down-regulation of the CAR protein using a drug-regulatable system or treatment with the multikinase inhibitor dasatinib resulted in the acquisition of a memory-like phenotype, global transcriptional and epigenetic reprogramming, and restored antitumor functionality in exhausted CAR-T cells. This work demonstrates that rest can enhance CAR-T cell efficacy by preventing or reversing exhaustion, and it challenges the notion that exhaustion is an epigenetically fixed state.

Conditional expression of HGAL leads to the development of diffuse large B-cell lymphoma in mice

Diffuse large B-cell lymphomas (DLBCLs) are clinically and genetically heterogeneous tumors. Deregulation of diverse biological processes specific to B cells, such as B-cell receptor (BCR) signaling and motility regulation, contribute to lymphomagenesis. Human germinal center associated lymphoma (HGAL) is a B-cell-specific adaptor protein controlling BCR signaling and B lymphocyte motility. In normal B cells, it is expressed in germinal center (GC) B lymphocytes and promptly downregulated upon further differentiation. The majority of DLBCL tumors, primarily GC B-cell types, but also activated types, express HGAL. To investigate the consequences of constitutive expression of HGAL in vivo, we generated mice that conditionally express human HGAL at different stages of hematopoietic development using 3 restricted Cre-mediated approaches to initiate expression of HGAL in hematopoietic stem cells, pro-B cells, or GC B cells. Following immune stimulation, we observed larger GCs in mice in which HGAL expression was initiated in GC B cells. All 3 mouse strains developed DLBCL at a frequency of 12% to 30% starting at age 13 months, leading to shorter survival. Immunohistochemical studies showed that all analyzed tumors were of the GC B-cell type. Exon sequencing revealed mutations reported in human DLBCL. Our data demonstrate that constitutive enforced expression of HGAL leads to DLBCL development.

Abstract PR007: Lymph node colonization promotes distant tumor metastasis through the induction of systemic tumor-specific immunosuppression

The majority of cancer-associated deaths result from distant organ metastasis rather than the primary tumor, yet the mechanisms that enable this process remain poorly understood. For most solid tumors, colonization of regional or distant lymph nodes (LNs) typically precedes the formation of distant organ metastases, yet it remains unclear whether LN metastasis plays a functional role in disease progression. LNs are education hubs of the adaptive immune system wherein antigens derived from pathogens or malignancies are presented to lymphocytes to elicit an adaptive immune response. Nonetheless, LN metastasis, which is typically attributed to passive drainage of tumor cells through lymphatics, frequently does not lead to the generation of anti-tumor immunity, but instead correlates with further disease progression. Here, we find that LN metastasis represents a critical step in tumor progression through the capacity of such metastases to induce tumor-specific immunosuppression in a manner that promotes further dissemination of tumors to distant organs. Using an in vivo passaging approach of a non-metastatic syngeneic melanoma, we generated 300 unique cell lines exhibiting varying degrees of LN metastatic capacity. We show that the presence of these LN metastases enables distant organ seeding of metastases in a manner that the parental tumor cannot, and this effect is eliminated in mice lacking an adaptive immune response. Furthermore, this promotion of distant seeding by LN metastases is tumor specific. Using flow cytometry and single-cell sequencing to perform organism-wide immune profiling, we identify multiple cellular mediators of tolerance. In particular, we find that LN metastases have the capacity to both resist NK cell cytotoxicity and induce regulatory T cells (Tregs) in vitro. Furthermore, depletion of NK cells in vivo enables non-metastatic tumors to disseminate to LNs, and ablation of Tregs using FoxP3-DTR mice eliminates the occurrence of lymphatic metastases. The immunosuppressive effects of LN metastases can be transferred to tumor-naïve recipients through specific lymphocyte populations in a manner that promotes distant seeding. Through the use of whole exome sequencing, we show that neither the metastatic proclivity nor immunosuppression evolve through the acquisition of driver mutations, loss of neoantigens, loss of MHC class I presentation, or decreases in melanoma antigen expression. Rather, by RNA-seq and ATAC-seq, we show that a conserved interferon signaling axis is upregulated in LN metastases and is rendered stable through epigenetic regulation of chromatin accessibility. Furthermore, using CRISPR/Cas9, we find that these pathways are required for LN metastatic seeding, and validate their conserved significance in additional mouse models of pancreatic ductal adenocarcinoma and head and neck squamous cell carcinoma (HNSCC), along with RNA-seq analysis of malignant populations sorted from HNSCC patients. Together, these findings demonstrate a critical role for LN metastasis in promoting tumor-specific immunosuppression.

Citation Format: Nathan E. Reticker-Flynn, Weiruo Zhang, Serena Chang, Andrew J. Gentles, John B. Sunwoo, Christina S. Kong, Sylvia K. Plevritis, Edgar G. Engleman. Lymph node colonization promotes distant tumor metastasis through the induction of systemic tumor-specific immunosuppression [abstract]. In: Proceedings of the AACR Virtual Special Conference on the Evolving Tumor Microenvironment in Cancer Progression: Mechanisms and Emerging Therapeutic Opportunities; in association with the Tumor Microenvironment (TME) Working Group; 2021 Jan 11-12. Philadelphia (PA): AACR; Cancer Res 2021;81(5 Suppl):Abstract nr PR007.

Abstract 3443: Atlas of clinically-distinct cell states and cellular ecosystems across human solid tumors

Tumors are complex ecosystems consisting of malignant, immune, and stromal elements whose dynamic interactions drive patient survival and response to therapy. A comprehensive understanding of the diversity of cellular states within the tumor microenvironment (TME), and their patterns of co-occurrence, could provide new diagnostic tools for improved disease management and novel targets for therapeutic intervention. To address this challenge, we developed EcoTyper, a novel machine learning framework for large-scale identification of TME cell states and their co-association patterns from bulk, single-cell, and spatially resolved tumor expression data. EcoTyper starts by “purifying” cell type-specific gene expression profiles of epithelial cells, immune, and stromal cell types from bulk tissue transcriptomes using CIBERSORTx (Newman et al., Nat Biotechnol 2019). It then identifies transcriptional states for each cell type, validates them in scRNA-seq data, and uncovers co-occurrence patterns between cell states in order to define tumor cellular ecosystems. Applied to 6,475 tumor and adjacent normal samples from solid tumor types profiled by The Cancer Genome Atlas (TCGA), EcoTyper identified robust transcriptional states across 12 major cell types, including epithelial, fibroblast, endothelial, and 9 immune subsets. These states included both known and novel cellular phenotypes, nearly all of which could be validated in a compendium of scRNA-seq tumor atlases spanning ~140,000 cells. Most cell states were specific to neoplastic tissue, ubiquitous across tumor types, and significantly associated with overall survival, both in TCGA and in 9,062 held-out tumor specimens (Gentles/Newman et al., Nat Medicine 2015). We found that specific cell states co-occur in distinct cellular communities with characteristic patterns of ligand-receptor interactions, genomic features, clinical outcomes, and spatial organization. One such ecosystem defined a normal-like state that was strongly enriched in non-malignant samples. Others delineated novel pro- and anti-tumor inflammatory environments involving specific fibroblast, endothelial, and immune cell transcriptional programs. In summary, large-scale deconvolution of cell type-specific transcriptomes across thousands of solid tumors revealed a comprehensive atlas of TME cell states and cellular ecosystems. Our results provide a high-resolution portrait of cellular heterogeneity in the TME across multiple solid tumor types, with implications for novel diagnostics and immunotherapeutic targets.

References:

1. Newman, A.M., et al., Determining cell type abundance and expression from bulk tissues with digital cytometry. Nat Biotechnol, 2019. 37(7): p. 773-782.

2. Gentles, A.J., et al., The prognostic landscape of genes and infiltrating immune cells across human cancers. Nature medicine, 2015. 21(8): p. 938.

Citation Format: Bogdan A. Luca, Chloé B. Steen, Armon Azizi, Magdalena Matusiak, Joanna Przybyl, Nastaran Neishaboori, Almudena Espín Pérez, Maximilian Diehn, Ash A. Alizadeh, Matt van de Rijn, Andrew J. Gentles, Aaron M. Newman. Atlas of clinically-distinct cell states and cellular ecosystems across human solid tumors [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 3443.

Abstract 3419: Lymph node colonization promotes distant tumor metastasis through the induction of tumor-specific immunosuppression

The majority of cancer-associated deaths result from distant organ metastasis rather than the primary tumor, yet the mechanisms that enable this process remain poorly understood. For most solid tumors, colonization of regional or distant lymph nodes (LNs) typically precedes the formation of distant organ metastases, yet it remains unclear whether LN metastasis plays a functional role in disease progression. LNs are education hubs of the adaptive immune system wherein antigens derived from pathogens or malignancies are presented to lymphocytes to elicit an adaptive immune response. Nonetheless, LN metastasis, which is typically attributed to passive drainage of tumor cells through lymphatics, frequently does not lead to the generation of anti-tumor immunity, but instead correlates with further disease progression. Here, we find that LN metastasis represents a critical step in tumor progression through the capacity of such metastases to induce tumor-specific immunosuppression in a manner that promotes further dissemination of tumors to distant organs. Using an in vivo passaging approach of a non-metastatic syngeneic melanoma, we generated 300 unique cell lines exhibiting varying degrees of LN metastatic capacity. Transcriptional profiling of the lines reveals a conserved enrichment for immune-related programs. We show that the presence of these LN metastases enables distant organ seeding of metastases in a manner that the parental tumor cannot, and this differential seeding is eliminated in mice lacking an adaptive immune response. Furthermore, this promotion of distant seeding by LN metastases is tumor specific. Using mass cytometry to perform organism-wide immune profiling, we identify multiple cellular mediators of tolerance. In particular, we find that LN metastases have the capacity to both resist NK cell cytotoxicity and induce regulatory T cells (Tregs) in vitro. Furthermore, depletion of NK cells in vivo enables non-metastatic tumors to disseminate to LNs, and ablation of Tregs using FoxP3-DTR mice eliminates the occurrence of lymphatic metastases. Through the use of whole exome sequencing, we show that neither the metastatic proclivity nor immunosuppression evolve through the acquisition of driver mutations, loss of neoantigens, loss of MHC class I presentation, or decreases in melanoma antigen expression. Rather, by RNA-seq and ATAC-seq, we show that a conserved interferon signaling axis is upregulated in LN metastases and is rendered stable through epigenetic regulation of chromatin accessibility. Furthermore, using CRISPR/Cas9, we find that these pathways are required for LN metastatic seeding, and validate their conserved significance in additional mouse models of pancreatic ductal adenocarcinoma and head and neck squamous cell carcinoma (HNSCC), along with RNA-seq analysis of malignant populations sorted from HNSCC patients. Together, these findings demonstrate a critical role for LN metastasis in promoting tumor-specific immunosuppression.

Citation Format: Nathan E. Reticker-Flynn, Pamela A. Basto, Weiruo Zhang, Maria M. Martins, Serena Chang, Andrew J. Gentles, John B. Sunwoo, Sylvia K. Plevritis, Edgar G. Engleman. Lymph node colonization promotes distant tumor metastasis through the induction of tumor-specific immunosuppression [abstract]. In: Proceedings of the Annual Meeting of the American Association for Cancer Research 2020; 2020 Apr 27-28 and Jun 22-24. Philadelphia (PA): AACR; Cancer Res 2020;80(16 Suppl):Abstract nr 3419.

Targetable genetic alterations of TCF4 (E2-2) drive immunoglobulin expression in diffuse large B cell lymphoma

The activated B cell (ABC-like) subtype of diffuse large B cell lymphoma (DLBCL) is characterized by chronic activation of signaling initiated by immunoglobulin μ (IgM). By analyzing the DNA copy number profiles of 1000 DLBCL tumors, we identified gains of 18q21.2 as the most frequent genetic alteration in ABC-like DLBCL. Using integrative analysis of matched gene expression profiling data, we found that the TCF4 (E2-2) transcription factor gene was the target of these alterations. Overexpression of TCF4 in ABC-like DLBCL cell lines led to its occupancy on immunoglobulin (IGHM) and MYC gene enhancers and increased expression of these genes at the transcript and protein levels. Inhibition of TCF4 activity with dominant-negative constructs was synthetically lethal to ABC-like DLBCL cell lines harboring TCF4 DNA copy gains, highlighting these gains as an attractive potential therapeutic target. Furthermore, the TCF4 gene was one of the top BRD4-regulated genes in DLBCL cell lines. BET proteolysis-targeting chimera (PROTAC) ARV771 extinguished TCF4, MYC, and IgM expression and killed ABC-like DLBCL cells in vitro. In DLBCL xenograft models, ARV771 treatment reduced tumor growth and prolonged survival. This work highlights a genetic mechanism for promoting immunoglobulin signaling in ABC-like DLBCL and provides a functional rationale for the use of BET inhibitors in this disease.

Maternal Anti-Dengue IgG Fucosylation Predicts Susceptibility to Dengue Disease in Infants

Infant mortality from dengue disease is a devastating global health burden that could be minimized with the ability to identify susceptibility for severe disease prior to infection. Although most primary infant dengue infections are asymptomatic, maternally derived anti-dengue immunoglobulin G (IgGs) present during infection can trigger progression to severe disease through antibody-dependent enhancement mechanisms. Importantly, specific characteristics of maternal IgGs that herald progression to severe infant dengue are unknown. Here, we define ≥10% afucosylation of maternal anti-dengue IgGs as a risk factor for susceptibility of infants to symptomatic dengue infections. Mechanistic experiments show that afucosylation of anti-dengue IgGs promotes FcγRIIIa signaling during infection, in turn enhancing dengue virus replication in FcγRIIIa+ monocytes. These studies identify a post-translational modification of anti-dengue IgGs that correlates with risk for symptomatic infant dengue infections and define a mechanism by which afucosylated antibodies and FcγRIIIa enhance dengue infections.

MYC and Twist1 cooperate to drive metastasis by eliciting crosstalk between cancer and innate immunity

Metastasis is a major cause of cancer mortality. We generated an autochthonous transgenic mouse model whereby conditional expression of MYC and Twist1 enables hepatocellular carcinoma (HCC) to metastasize in >90% of mice. MYC and Twist1 cooperate and their sustained expression is required to elicit a transcriptional program associated with the activation of innate immunity, through secretion of a cytokinome that elicits recruitment and polarization of tumor associated macrophages (TAMs). Systemic treatment with Ccl2 and Il13 induced MYC-HCCs to metastasize; whereas, blockade of Ccl2 and Il13 abrogated MYC/Twist1-HCC metastasis. Further, in 33 human cancers (n = 9502) MYC and TWIST1 predict poor survival (p=4.3×10⁻¹⁰), CCL2/IL13 expression (p<10⁻¹⁰⁹) and TAM infiltration (p<10⁻⁹⁶). Finally, in the plasma of patients with HCC (n = 25) but not cirrhosis (n = 10), CCL2 and IL13 were increased and IL13 predicted invasive tumors. Therefore, MYC and TWIST1 generally appear to cooperate in human cancer to elicit a cytokinome that enables metastasis through crosstalk between cancer and immune microenvironment.

Cancer develops when cells in the body gain mutations that allow them to grow and divide rapidly and uncontrollably. As the disease progresses these cancer cells develop the ability to spread around the body. This process of spreading, called metastasis, is responsible for most cancer-related deaths in humans, but no current treatments target it.

Mutations that increase the levels of two proteins known as MYC and TWIST1 in cells cause many human cancers. In healthy adult cells, normal levels of MYC and TWIST1 act as key regulators that switch thousands of genes on or off. TWIST1 is known to control the movement and spread of cells in the embryo. However, it is not known how MYC and TWIST1 work together to promote the metastasis of cancer cells.

To address this question, Dhanasekaran, Baylot et al. used mice to investigate the roles of MYC and TWIST1 in the metastasis of cancer cells. The experiments showed that these two proteins work together to reprogram mouse cancer cells to release signal molecules known as cytokines. These molecules convert immune cells known as macrophages to a tumor-friendly state that allows cancers cells to spread around the body. Inhibiting two cytokines known as CCL2 and IL13 prevented the cancer cells from moving.

Further experiments analyzed tumor samples from around 10,000 human patients with 33 different cancers. This revealed that patients that had higher levels of MYC and TWIST1 proteins in their tumors also had increased levels of CCL2 and IL13, more activated macrophages and were less likely to recover from their cancer.

The findings of Dhanasekaran, Baylot et al. suggest that MYC and TWIST1 may instigate metastasis in many human cancers, and therapies targeting specific cytokines may prevent these cancers from spreading around the body. Furthermore, screening blood for the levels of cytokines may help to identify the cancer patients who would benefit from such therapies.

Prognostic profiling of the immune cell microenvironment in Ewing´s Sarcoma Family of Tumors

Ewing´s Sarcoma Family of Tumors (ESFT) are clinically aggressive bone and soft tissue tumors in children and young adults. Analysis of the immune tumor microenvironment (TME) provides insight into tumor evolution and novel treatment options. So far, the scarcity of immune cells in ESFT has hindered a comprehensive analysis of rare subtypes. We determined the relative fraction of 22 immune cell types using 197 microarray gene expression datasets of primary ESFT tumor samples by using CIBERSORT, a deconvolution algorithm enumerating infiltrating leucocytes in bulk tumor tissue. The most abundant cells were macrophages (mean 43% of total tumor-infiltrating leukocytes, TILs), predominantly immunosuppressive M2 type macrophages, followed by T cells (mean 23% of TILs). Increased neutrophils, albeit at low number, were associated with a poor overall survival (OS) (p = .038) and increased M2 macrophages predicted a shorter event-free survival (EFS) (p = .033). High frequency of T cells and activated NK cells correlated with prolonged OS (p = .044 and p = .007, respectively). A small patient population (9/32) with combined low infiltrating M2 macrophages, low neutrophils, and high total T cells was identified with favorable outcome. This finding was confirmed in a validation cohort of patients with follow up (11/38). When comparing the immune TME with expression of known stemness genes, hypoxia-inducible factor 1 α (HIF1α) correlated with high abundance of macrophages and neutrophils and decreased T cell levels. The immune TME in ESFTs shows a distinct composition including rare immune cell subsets that in part may be due to expression of HIF1α.

Perfusion MRI-Based Fractional Tumor Burden Differentiates between Tumor and Treatment Effect in Recurrent Glioblastomas and Informs Clinical Decision-Making

BACKGROUND AND PURPOSE

Fractional tumor burden better correlates with histologic tumor volume fraction in treated glioblastoma than other perfusion metrics such as relative CBV. We defined fractional tumor burden classes with low and high blood volume to distinguish tumor from treatment effect and to determine whether fractional tumor burden can inform treatment-related decision-making.

MATERIALS AND METHODS

Forty-seven patients with high-grade gliomas (primarily glioblastoma) with recurrent contrast-enhancing lesions on DSC-MR imaging were retrospectively evaluated after surgical sampling. Histopathologic examination defined treatment effect versus tumor. Normalized relative CBV thresholds of 1.0 and 1.75 were used to define low, intermediate, and high fractional tumor burden classes in each histopathologically defined group. Performance was assessed with an area under the receiver operating characteristic curve. Consensus agreement among physician raters reporting hypothetic changes in treatment-related decisions based on fractional tumor burden was compared with actual real-time treatment decisions.

RESULTS

Mean lower fractional tumor burden, high fractional tumor burden, and relative CBV of the contrast-enhancing volume were significantly different between treatment effect and tumor (P = .002, P < .001, and P < .001), with tumor having significantly higher fractional tumor burden and relative CBV and lower fractional tumor burden. No significance was found with intermediate fractional tumor burden. Performance of the area under the receiver operating characteristic curve was the following: high fractional tumor burden, 0.85; low fractional tumor burden, 0.7; and relative CBV, 0.81. In comparing treatment decisions, there were disagreements in 7% of tumor and 44% of treatment effect cases; in the latter, all disagreements were in cases with scattered atypical cells.

CONCLUSIONS

High fractional tumor burden and low fractional tumor burden define fractions of the contrast-enhancing lesion volume with high and low blood volume, respectively, and can differentiate treatment effect from tumor in recurrent glioblastomas. Fractional tumor burden maps can also help to inform clinical decision-making.

LMO2 Confers Synthetic Lethality to PARP Inhibition in DLBCL

Deficiency in DNA double-strand break (DSB) repair mechanisms has been widely exploited for the treatment of different malignances, including homologous recombination (HR)-deficient breast and ovarian cancers. Here we demonstrate that diffuse large B cell lymphomas (DLBCLs) expressing LMO2 protein are functionally deficient in HR-mediated DSB repair. Mechanistically, LMO2 inhibits BRCA1 recruitment to DSBs by interacting with 53BP1 during repair. Similar to BRCA1-deficient cells, LMO2-positive DLBCLs and T cell acute lymphoblastic leukemia (T-ALL) cells exhibit a high sensitivity to poly(ADP-ribose) polymerase (PARP) inhibitors. Furthermore, chemotherapy and PARP inhibitors synergize to inhibit the growth of LMO2-positive tumors. Together, our results reveal that LMO2 expression predicts HR deficiency and the potential therapeutic use of PARP inhibitors in DLBCL and T-ALL.

Abstract 2703: Lymph node colonization promotes distant tumor metastasis through the induction of systemic immune tolerance

The majority of cancer-associated deaths are the result of distant organ metastasis, an event that is typically preceded by metastasis to regional or distant lymph nodes (LNs). LNs are education hubs of the adaptive immune system wherein antigens derived from pathogens or malignancies are presented to lymphocytes in a manner that facilitates elimination of the threat. Nonetheless, LN metastasis, which is typically attributed to passive drainage of tumor cells through lymphatics, frequently does not lead to the generation of an anti-tumor immune response, but instead correlates with poor prognosis and further disease progression. Here, we find that LN metastasis represents a critical step in tumor progression through the capacity of such metastases to induce systemic immune tolerance in a manner that promotes further dissemination of tumors to distant organs. Through serial in vivo passaging of a syngeneic melanoma in mice, we generate nearly 300 unique cell lines that exhibit an enhanced capacity to metastasize to LNs. Transcriptional profiling of these lines reveals increased expression of immune-related programs. We show that the presence of these LN metastases enables distant organ seeding of metastases in a manner that the parental tumor cannot, and this differential seeding is eliminated in mice that lack an adaptive immune response. To query the effects of the LN metastases on the systemic immune response, we perform organism-wide immune profiling by mass cytometry and identify a number of cellular mediators of tolerance. In particular, we find that LN metastases have the capacity to both resist NK cell cytotoxicity and induce regulatory T cells (Tregs) in vitro. Furthermore, depletion of NK cells in vivo enables non-metastatic tumors to disseminate to LNs, and ablation of Tregs using FoxP3-DTR mice eliminates the occurrence of lymphatic metastases. We further identify an interferon signaling axis that is constitutively activated within the LN metastases in the absence of exogenous interferon signaling. Through the use of ATAC-seq, we find that this program is conferred through epigenetic regulation of chromatin accessibility. Knockout of key interferon-induced genes using CRISPR/Cas9 in the LN-metastatic cells reveals that this program is required for enhanced LN metastatic seeding in vivo, and their overexpression increases LN metastasis of the non-metastatic cells. Using additional mouse models of pancreatic ductal adenocarcinoma and head and neck squamous cell carcinoma (HNSCC), we show that these findings are conserved across multiple malignancies. Additionally, we perform RNA-seq on sorted malignant populations from node-positive and node-negative HNSCC patients and confirm that these differences in transcriptional profiles extend to the human disease. Together, these findings demonstrate a critical role for LN metastasis in promoting tumor immune tolerance.

Citation Format: Nathan E. Reticker-Flynn, Maria M. Martins, Pamela A. Basto, Weiruo Zhang, Alborz Bejnood, Andrew J. Gentles, John B. Sunwoo, Sylvia K. Plevritis, Edgar G. Engleman. Lymph node colonization promotes distant tumor metastasis through the induction of systemic immune tolerance [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2019; 2019 Mar 29-Apr 3; Atlanta, GA. Philadelphia (PA): AACR; Cancer Res 2019;79(13 Suppl):Abstract nr 2703.

Determining cell type abundance and expression from bulk tissues with digital cytometry

Single-cell RNA-sequencing has emerged as a powerful technique for characterizing cellular heterogeneity, but it is currently impractical on large sample cohorts and cannot be applied to fixed specimens collected as part of routine clinical care. We previously developed an approach for digital cytometry, called CIBERSORT, that enables estimation of cell type abundances from bulk tissue transcriptomes. We now introduce CIBERSORTx, a machine learning method that extends this framework to infer cell-type-specific gene expression profiles without physical cell isolation. By minimizing platform-specific variation, CIBERSORTx also allows the use of single-cell RNA-sequencing data for large-scale tissue dissection. We evaluated the utility of CIBERSORTx in multiple tumor types, including melanoma, where single-cell reference profiles were used to dissect bulk clinical specimens, revealing cell-type-specific phenotypic states linked to distinct driver mutations and response to immune checkpoint blockade. We anticipate that digital cytometry will augment single-cell profiling efforts, enabling cost-effective, high-throughput tissue characterization without the need for antibodies, disaggregation or viable cells.

Machine Learning Identifies Stemness Features Associated with Oncogenic Dedifferentiation

Cancer progression involves the gradual loss of a differentiated phenotype and acquisition of progenitor and stem-cell-like features. Here, we provide novel stemness indices for assessing the degree of oncogenic dedifferentiation. We used an innovative one-class logistic regression (OCLR) machine-learning algorithm to extract transcriptomic and epigenetic feature sets derived from non-transformed pluripotent stem cells and their differentiated progeny. Using OCLR, we were able to identify previously undiscovered biological mechanisms associated with the dedifferentiated oncogenic state. Analyses of the tumor microenvironment revealed unanticipated correlation of cancer stemness with immune checkpoint expression and infiltrating immune cells. We found that the dedifferentiated oncogenic phenotype was generally most prominent in metastatic tumors. Application of our stemness indices to single-cell data revealed patterns of intra-tumor molecular heterogeneity. Finally, the indices allowed for the identification of novel targets and possible targeted therapies aimed at tumor differentiation.

Abstract LB-373: Comprehensive analysis of cancer stemness

Cancer progression involves the gradual loss of a differentiated phenotype and acquisition of progenitor and stem cell-like features. Here, we provide new stemness indices for assessing the degree of oncogenic dedifferentiation. We took advantage of an innovative one-class logistic regression machine learning algorithm (OCLR) to extract transcriptomic and epigenetic feature sets derived from non-transformed pluripotent stem cells and their differentiated progenies. Using OCLR, we were able to sort TCGA tumor samples by stemness phenotype and identify previously undiscovered biological mechanisms associated with the dedifferentiated oncogenic state. Analyses of tumor microenvironment revealed the correlation of cancer stemness with immune checkpoint expression and infiltrating immune system cells not previously anticipated. We have shown the de-differentiated oncogenic phenotype increased in the metastatic tumor that further justify their more aggressive phenotype. Application of our stemness indices reveals features of intra-tumor heterogeneity in molecular profiles obtained from the single-cell analyses. Finally, the machine learning-based indices allowed for the identification of chemical compounds and novel targets for the cancer therapies aiming at tumor differentiation. Our findings provide new prognostic signatures that enable cancer biologists and oncologists to quantify the impact of tumor stemness on outcome across cancer types and may help to pave the way for progress in treatment strategies for cancer patients.

Citation Format: Tathiane M. Malta, Artem Sokolov, Andrew J. Gentles, Tomasz Burzykowski, Laila Poisson, John Weinstein, Bożena Kamińska, Joerg Huelsken, Larsson Omberg, Olivier Gevaert, Antonio Colaprico, Patrycja Czerwińska, Sylwia Mazurek, Lopa Mishra, Holger Heyn, Alex Krasnitz, Andrew K. Godwin, Alexander J. Lazar, The Cancer Genome Atlas Research Network, Joshua M. Stuart, Katherine Hoadley, Peter W. Laird, Houtan Noushmehr, Maciej Wiznerowicz. Comprehensive analysis of cancer stemness [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 LB-373.

GFPT2-Expressing Cancer-Associated Fibroblasts Mediate Metabolic Reprogramming in Human Lung Adenocarcinoma

Metabolic reprogramming of the tumor microenvironment is recognized as a cancer hallmark. To identify new molecular processes associated with tumor metabolism, we analyzed the transcriptome of bulk and flow-sorted human primary non-small cell lung cancer (NSCLC) together with ¹⁸FDG-PET scans, which provide a clinical measure of glucose uptake. Tumors with higher glucose uptake were functionally enriched for molecular processes associated with invasion in adenocarcinoma and cell growth in squamous cell carcinoma (SCC). Next, we identified genes correlated to glucose uptake that were predominately overexpressed in a single cell-type comprising the tumor microenvironment. For SCC, most of these genes were expressed by malignant cells, whereas in adenocarcinoma, they were predominately expressed by stromal cells, particularly cancer-associated fibroblasts (CAF). Among these adenocarcinoma genes correlated to glucose uptake, we focused on glutamine-fructose-6-phosphate transaminase 2 (GFPT2), which codes for the glutamine-fructose-6-phosphate aminotransferase 2 (GFAT2), a rate-limiting enzyme of the hexosamine biosynthesis pathway (HBP), which is responsible for glycosylation. GFPT2 was predictive of glucose uptake independent of GLUT1, the primary glucose transporter, and was prognostically significant at both gene and protein level. We confirmed that normal fibroblasts transformed to CAF-like cells, following TGFβ treatment, upregulated HBP genes, including GFPT2, with less change in genes driving glycolysis, pentose phosphate pathway, and TCA cycle. Our work provides new evidence of histology-specific tumor stromal properties associated with glucose uptake in NSCLC and identifies GFPT2 as a critical regulator of tumor metabolic reprogramming in adenocarcinoma.Significance: These findings implicate the hexosamine biosynthesis pathway as a potential new therapeutic target in lung adenocarcinoma. Cancer Res; 78(13); 3445-57. ©2018 AACR.

The Immune Landscape of Cancer

We performed an extensive immunogenomic analysis of more than 10,000 tumors comprising 33 diverse cancer types by utilizing data compiled by TCGA. Across cancer types, we identified six immune subtypes-wound healing, IFN-γ dominant, inflammatory, lymphocyte depleted, immunologically quiet, and TGF-β dominant-characterized by differences in macrophage or lymphocyte signatures, Th1:Th2 cell ratio, extent of intratumoral heterogeneity, aneuploidy, extent of neoantigen load, overall cell proliferation, expression of immunomodulatory genes, and prognosis. Specific driver mutations correlated with lower (CTNNB1, NRAS, or IDH1) or higher (BRAF, TP53, or CASP8) leukocyte levels across all cancers. Multiple control modalities of the intracellular and extracellular networks (transcription, microRNAs, copy number, and epigenetic processes) were involved in tumor-immune cell interactions, both across and within immune subtypes. Our immunogenomics pipeline to characterize these heterogeneous tumors and the resulting data are intended to serve as a resource for future targeted studies to further advance the field.

Subtype assignment of CLL based on B-cell subset associated gene signatures from normal bone marrow - A proof of concept study

Diagnostic and prognostic evaluation of chronic lymphocytic leukemia (CLL) involves blood cell counts, immunophenotyping, IgVH mutation status, and cytogenetic analyses. We generated B-cell associated gene-signatures (BAGS) based on six naturally occurring B-cell subsets within normal bone marrow. Our hypothesis is that by segregating CLL according to BAGS, we can identify subtypes with prognostic implications in support of pathogenetic value of BAGS. Microarray-based gene-expression samples from eight independent CLL cohorts (1,024 untreated patients) were BAGS-stratified into pre-BI, pre-BII, immature, naïve, memory, or plasma cell subtypes; the majority falling within the memory (24.5-45.8%) or naïve (14.5-32.3%) categories. For a subset of CLL patients (n = 296), time to treatment (TTT) was shorter amongst early differentiation subtypes (pre-BI/pre-BII/immature) compared to late subtypes (memory/plasma cell, HR: 0.53 [0.35-0.78]). Particularly, pre-BII subtype patients had the shortest TTT among all subtypes. Correlates derived for BAGS subtype and IgVH mutation (n = 405) revealed an elevated mutation frequency in late vs. early subtypes (71% vs. 45%, P < .001). Predictions for BAGS subtype resistance towards rituximab and cyclophosphamide varied for rituximab, whereas all subtypes were sensitive to cyclophosphamide. This study supports our hypothesis that BAGS-subtyping may be of tangible prognostic and pathogenetic value for CLL patients.

Brd4 regulates the expression of essential autophagy genes and Keap1 in AML cells

We have recently reported that activation of Brd4 is associated with the presence of autophagy in NPMc+ and MLL AML cells. In order to determine the mechanisms underlying this relationship, we have examined the role of Brd4 in regulating the expression of several genes that are central to the process of autophagy. We found that Brd4 binds to the promoters of ATG 3, 7 and CEBPβ, and expression of these genes is markedly reduced by inhibitors of Brd4, as well as by Brd4-shRNA and depletion of CEBPβ. Inhibitors of Brd4 also dramatically suppress the transcription of Keap1, thereby increasing the expression of anti-oxidant genes through the Nrf2 pathway and reducing the cytotoxicity induced by Brd4 inhibitors. Elimination of ATG3 or KEAP1 expression using CRISPR-cas9 mediated genomic editing markedly reduced autophagy. We conclude that Brd4 plays a significant role in autophagy activation through the direct transcriptional regulation of genes essential for it, as well as through the Keap1-Nrf2 axis in NPMc+ and MLL-fusion AML cells.

Module Analysis Captures Pancancer Genetically and Epigenetically Deregulated Cancer Driver Genes for Smoking and Antiviral Response

The availability of increasing volumes of multi-omics profiles across many cancers promises to improve our understanding of the regulatory mechanisms underlying cancer. The main challenge is to integrate these multiple levels of omics profiles and especially to analyze them across many cancers. Here we present AMARETTO, an algorithm that addresses both challenges in three steps. First, AMARETTO identifies potential cancer driver genes through integration of copy number, DNA methylation and gene expression data. Then AMARETTO connects these driver genes with co-expressed target genes that they control, defined as regulatory modules. Thirdly, we connect AMARETTO modules identified from different cancer sites into a pancancer network to identify cancer driver genes. Here we applied AMARETTO in a pancancer study comprising eleven cancer sites and confirmed that AMARETTO captures hallmarks of cancer. We also demonstrated that AMARETTO enables the identification of novel pancancer driver genes. In particular, our analysis led to the identification of pancancer driver genes of smoking-induced cancers and 'antiviral' interferon-modulated innate immune response.

SOFTWARE AVAILABILITY

AMARETTO is available as an R package at https://bitbucket.org/gevaertlab/pancanceramaretto.

Low BUB1 expression is an adverse prognostic marker in gastric adenocarcinoma

Gastric adenocarcinomas are associated with a poor prognosis due to the fact that the tumor has often metastasized by the time of diagnosis and prognostic markers are urgently needed to tailor treatment.

We examined the expression of the mitotic spindle checkpoint protein BUB1 (budding uninhibited by benzimidazoles 1) and Ki-67 protein expression by immunohistochemistry in 218 patients with primary gastric adenocarcinomas.

Tumors with low frequency of BUB1 expression were associated with larger tumor size (pT) (p < 0.001), higher incidence of lymph node metastases (pN) (p = 0.027), distant metastases (pM) (p = 0.006) and higher UICC stage (p < 0.001). Furthermore, BUB1 expression was inversely correlated with residual tumor stage (p = 0.038). Abundant BUB1 protein expression correlated with frequent Ki-67 protein expression (p < 0.001) and low BUB1 expression was associated with shorter survival (p < 0.001). Univariate and multivariate analyses confirmed BUB1 to be an independent prognostic marker in gastric cancer (p = 0.021).

Data normalization considerations for digital tumor dissection

In a recently published article in Genome Biology, Li and colleagues introduced TIMER, a gene expression deconvolution approach for studying tumor-infiltrating leukocytes (TILs) in 23 cancer types profiled by The Cancer Genome Atlas. Methods to characterize TIL biology are increasingly important, and the authors offer several arguments in favor of their strategy. Several of these claims warrant further discussion and highlight the critical importance of data normalization in gene expression deconvolution applications.Please see related Li et al correspondence: www.dx.doi.org/10.1186/s13059-017-1256-5 and Zheng correspondence: www.dx.doi.org/10.1186/s13059-017-1258-3.

Abstract LB-219: Higher levels of mast cells associate with favorable outcomes in non-small cell lung cancer and correlate with lower malignant cell proliferation

The tumor microenvironment (TME) involves complex interactions between malignant and stromal cell types. Much of our knowledge of cancer biology has been derived from studying molecular mechanisms underlying bulk tumors, with a focus on specific malignant pathways that have become dysregulated during tumorigenesis and tumor progression. Although studies have identified interactions between malignant and stromal cells within the TME, few have sought to comprehensively identify such relationships. Here, we performed RNA-seq on bulk and flow-sorted non-small cell primary human lung tumors enriching for malignant cells, endothelial cells, immune cells, and fibroblasts. We derived a map of cell-specific differential gene expression of prognostically associated secreted factors and cell surface markers, and computationally reconstructed pairwise cross-talk between cell types. We found significant novel associations between transcriptional profiles of malignant populations and specific stromal populations, focusing here on mast cells.

We identified presence of infiltrating mast cells to be negatively correlated with malignant cell proliferation. Expression of TPSAB1 (Tryptase Alpha/Beta 1) is largely confined to mast cells, and its high expression in both adenocarcinoma and squamous cell carcinoma is favorably prognostic in both histologies across multiple datasets based on gene expression data. We validated the prognostic relevance of mast cells in NSCLC by immunohistochemical (IHC) staining of a lung tumor tissue microarray (TMA) for MCT (mast cell tryptase). Higher mast cell count was associated with better overall survival across all NSCLC, or when considering either adenocarcinoma or squamous cell carcinoma alone. When mast cell counts were quantified as “High”, “Intermediate”, “Low”, and “Negative” levels without reference to clinical outcomes, these were statistically significantly associated with survival. Negative-, low-, and intermediate-levels of mast cells all conferred worse prognosis than high mast-cell levels. Mast cell levels remained prognostic in multivariate analysis with independent clinical factors including age, stage, and gender. Finally, we directly evaluated the relationship of mast cell levels to tumor proliferation by staining the same samples for the proliferation marker KI67. The percentage of tumor cells staining positive for KI67 was lower in tumors with high vs low numbers of mast cells (p=0.048, Wilcox rank-sum test).

In summary, we identified and validated a specific inverse relationship between levels of infiltrating mast cells and malignant cell proliferation. These results illustrate the utility of transcriptomic profiling of flow-sorted subpopulations from solid tumors in order to identify tumor-microenvironment interactions that may have prognostic and therapeutic relevance.

Citation Format: Andrew J. Gentles, Angela Hui, Weiguo Feng, Ramesh V. Nair, Alice Yu, Majid Shafiq, Erna Forgo, Amanda Khuong, Yue Xu, Chuong D. Hoang, Robert B. West, Matt van de Rijn, Maximilian Diehn, Sylvia K. Plevritis. Higher levels of mast cells associate with favorable outcomes in non-small cell lung cancer and correlate with lower malignant cell proliferation [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr LB-219. doi:10.1158/1538-7445.AM2017-LB-219

Abstract LB-004: Molecular hallmarks of cancer: Stemness

Stemness, defined as the potential for self-renewal and de-differentiation from the cell-of-origin, has been initially attributed to normal stem cells that possess the ability to give rise to all cell types in the adult organism. Cancer stem cells have been identified in most if not all hematological malignancies as well as in solid tumors. CSCs are postulated to be amongst the most resistant cells to various forms of chemotherapy and radiotherapy and contribute to relapse and metastasis and hence inferior clinical outcomes.

Here, we have performed a comprehensive and multi-platform analyses of the stemness features in 33 cancer types totalling more than 10,000 mostly primary, but also some metastatic and recurrent samples. In the first step, we have derived signatures to measure the stemness using molecular profiles of normal cells with various degrees of stemness from publicly available datasets. By multiplatform analyses of the transcriptome, methylome, and chromatin markers using different machine learning computational approaches, we have obtained four independent stemness scores. Initial validation revealed comparable performance of these computational metrics for sorting the TCGA samples. The cancer types with previously documented features of de-differentiation have higher stemness score in contrast to more differentiated cancer types. More detailed analyses of recently published molecular subtypes of gliomas have revealed a strong association of high stemness scores with the worst clinical outcome, further confirming our initial hypotheses.

By adopting machine learning algorithms, we have counted percentage of immune cells of both innate and acquired response infiltrating tumor tissue. Further deconvolution of the TCGA cancers allowed quantification of the tumour microenvironment composed of cancer-associated fibroblasts and the endothelial cells. Integration of gene expression and DNA methylation datasets defines classic hallmarks of pluripotent signatures which shed new light into biological processes regulating and maintaining oncogenic dedifferentiation.

Our ongoing analyses of the TCGA samples sorted by the obtained metrics involve: (1) activity of the selected hallmarks of cancer that have been shown as features of the CSCs; (2) correlation with the tumor pathology grading and clinical outcomes; (3) identification of potential drivers for development and/or repositioning of drugs targeting tumor self-renewal and de-differentiation potential; and (4) development of novel biomarkers for selection of patients that will respond to these novel therapies.

In summary, we defined molecular signatures of cancer stem cells that enabled classification of TCGA cancer types and identification of tumors with stem/progenitor-like phenotypes, associated with poor clinical outcomes. Understanding stemness-like hallmark of cancer will pave the way for novel diagnostics and therapies for cancer patients.

Citation Format: Tathiane Malta, Artem Sokolov, Andrew J. Gentles, Tomasz Burzykowski, Olivier Gevaert, The Stemness Analysis Working Group TCGA PanCancerNetwork, Peter W. Laird, Houtan Noushmehr, Maciej Wiznerowicz. Molecular hallmarks of cancer: Stemness [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr LB-004. doi:10.1158/1538-7445.AM2017-LB-004

Identification of an atypical etiological head and neck squamous carcinoma subtype featuring the CpG island methylator phenotype

Head and neck squamous cell carcinoma (HNSCC) is broadly classified into HNSCC associated with human papilloma virus (HPV) infection, and HPV negative HNSCC, which is typically smoking-related. A subset of HPV negative HNSCCs occur in patients without smoking history, however, and these etiologically 'atypical' HNSCCs disproportionately occur in the oral cavity, and in female patients, suggesting a distinct etiology. To investigate the determinants of clinical and molecular heterogeneity, we performed unsupervised clustering to classify 528 HNSCC patients from The Cancer Genome Atlas (TCGA) into putative intrinsic subtypes based on their profiles of epigenetically (DNA methylation) deregulated genes. HNSCCs clustered into five subtypes, including one HPV positive subtype, two smoking-related subtypes, and two atypical subtypes. One atypical subtype was particularly genomically stable, but featured widespread gene silencing associated with the 'CpG island methylator phenotype' (CIMP). Further distinguishing features of this 'CIMP-Atypical' subtype include an antiviral gene expression profile associated with pro-inflammatory M1 macrophages and CD8+ T cell infiltration, CASP8 mutations, and a well-differentiated state corresponding to normal SOX2 copy number and SOX2OT hypermethylation. We developed a gene expression classifier for the CIMP-Atypical subtype that could classify atypical disease features in two independent patient cohorts, demonstrating the reproducibility of this subtype. Taken together, these findings provide unprecedented evidence that atypical HNSCC is molecularly distinct, and postulates the CIMP-Atypical subtype as a distinct clinical entity that may be caused by chronic inflammation.

Human AML-iPSCs Reacquire Leukemic Properties after Differentiation and Model Clonal Variation of Disease

Understanding the relative contributions of genetic and epigenetic abnormalities to acute myeloid leukemia (AML) should assist integrated design of targeted therapies. In this study, we generated induced pluripotent stem cells (iPSCs) from AML patient samples harboring MLL rearrangements and found that they retained leukemic mutations but reset leukemic DNA methylation/gene expression patterns. AML-iPSCs lacked leukemic potential, but when differentiated into hematopoietic cells, they reacquired the ability to give rise to leukemia in vivo and reestablished leukemic DNA methylation/gene expression patterns, including an aberrant MLL signature. Epigenetic reprogramming was therefore not sufficient to eliminate leukemic behavior. This approach also allowed us to study the properties of distinct AML subclones, including differential drug susceptibilities of KRAS mutant and wild-type cells, and predict relapse based on increased cytarabine resistance of a KRAS wild-type subclone. Overall, our findings illustrate the value of AML-iPSCs for investigating the mechanistic basis and clonal properties of human AML.

Role of KEAP1/NRF2 and TP53 Mutations in Lung Squamous Cell Carcinoma Development and Radiation Resistance

Lung squamous cell carcinoma (LSCC) pathogenesis remains incompletely understood, and biomarkers predicting treatment response remain lacking. Here, we describe novel murine LSCC models driven by loss of Trp53 and Keap1, both of which are frequently mutated in human LSCCs. Homozygous inactivation of Keap1 or Trp53 promoted airway basal stem cell (ABSC) self-renewal, suggesting that mutations in these genes lead to expansion of mutant stem cell clones. Deletion of Trp53 and Keap1 in ABSCs, but not more differentiated tracheal cells, produced tumors recapitulating histologic and molecular features of human LSCCs, indicating that they represent the likely cell of origin in this model. Deletion of Keap1 promoted tumor aggressiveness, metastasis, and resistance to oxidative stress and radiotherapy (RT). KEAP1/NRF2 mutation status predicted risk of local recurrence after RT in patients with non-small lung cancer (NSCLC) and could be noninvasively identified in circulating tumor DNA. Thus, KEAP1/NRF2 mutations could serve as predictive biomarkers for personalization of therapeutic strategies for NSCLCs.

SIGNIFICANCE

We developed an LSCC mouse model involving Trp53 and Keap1, which are frequently mutated in human LSCCs. In this model, ABSCs are the cell of origin of these tumors. KEAP1/NRF2 mutations increase radioresistance and predict local tumor recurrence in radiotherapy patients. Our findings are of potential clinical relevance and could lead to personalized treatment strategies for tumors with KEAP1/NRF2 mutations. Cancer Discov; 7(1); 86-101. ©2016 AACR.This article is highlighted in the In This Issue feature, p. 1.