Multi-modal molecular programs regulate melanoma cell state.

Melanoma cells display distinct intrinsic phenotypic states. Here, we seek to characterize the molecular regulation of these states using multi-omic analyses of whole exome, transcriptome, microRNA, long non-coding RNA and DNA methylation data together with reverse-phase protein array data on a panel of 68 highly annotated early passage melanoma cell lines. We demonstrate that clearly defined cancer cell intrinsic transcriptomic programs are maintained in melanoma cells ex vivo and remain highly conserved within melanoma tumors, are associated with distinct immune features within tumors, and differentially correlate with checkpoint inhibitor and adoptive T cell therapy efficacy. Through integrative analyses we demonstrate highly complex multi-omic regulation of melanoma cell intrinsic programs that provide key insights into the molecular maintenance of phenotypic states. These findings have implications for cancer biology and the identification of new therapeutic strategies. Further, these deeply characterized cell lines will serve as an invaluable resource for future research in the field.

Experimental models of undifferentiated pleomorphic sarcoma and malignant peripheral nerve sheath tumor.

Undifferentiated pleomorphic sarcoma (UPS) and malignant peripheral nerve sheath tumor (MPNST) are aggressive soft tissue sarcomas that do not respond well to current treatment modalities. The limited availability of UPS and MPNST cell lines makes it challenging to identify potential therapeutic targets in a laboratory setting. Understanding the urgent need for improved treatments for these tumors and the limited cellular models available, we generated additional cell lines to study these rare cancers. Patient-derived tumors were used to establish 4 new UPS models, including one radiation-associated UPS-UPS271.1, UPS511, UPS0103, and RIS620, one unclassified spindle cell sarcoma-USC060.1, and 3 new models of MPNST-MPNST007, MPNST3813E, and MPNST4970. This study examined the utility of the new cell lines as sarcoma models by assessing their tumorigenic potential and mutation status for known sarcoma-related genes. All the cell lines formed colonies and migrated in vitro. The in vivo tumorigenic potential of the cell lines and corresponding xenografts was determined by subcutaneous injection or xenograft re-passaging into immunocompromised mice. USC060.1 and UPS511 cells formed tumors in mice upon subcutaneous injection. UPS0103 and RIS620 tumor implants formed tumors in vivo, as did MPNST007 and MPNST3813E tumor implants. Targeted sequencing analysis of a panel of genes frequently mutated in sarcomas identified TP53, RB1, and ATRX mutations in a subset of the cell lines. These new cellular models provide the scientific community with powerful tools for detailed studies of tumorigenesis and for investigating novel therapies for UPS and MPNST.

The histologic phenotype of lung cancers is associated with transcriptomic features rather than genomic characteristics.

Histology plays an essential role in therapeutic decision-making for lung cancer patients. However, the molecular determinants of lung cancer histology are largely unknown. We conduct whole-exome sequencing and microarray profiling on 19 micro-dissected tumor regions of different histologic subtypes from 9 patients with lung cancers of mixed histology. A median of 68.9% of point mutations and 83% of copy number aberrations are shared between different histologic components within the same tumors. Furthermore, different histologic components within the tumors demonstrate similar subclonal architecture. On the other hand, transcriptomic profiling reveals shared pathways between the same histologic subtypes from different patients, which is supported by the analyses of the transcriptomic data from 141 cell lines and 343 lung cancers of different histologic subtypes. These data derived from mixed histologic subtypes in the setting of identical genetic background and exposure history support that the histologic fate of lung cancer cells is associated with transcriptomic features rather than the genomic profiles in most tumors.

Somatic Mutations in Circulating Cell-Free DNA and Risk for Hepatocellular Carcinoma in Hispanics.

Hispanics are disproportionally affected by liver fibrosis and hepatocellular carcinoma (HCC). Advanced liver fibrosis is a major risk factor for HCC development. We aimed at identifying somatic mutations in plasma cell-free DNA (cfDNA) of Hispanics with HCC and Hispanics with advanced liver fibrosis but no HCC. Targeted sequencing of over 262 cancer-associated genes identified nonsynonymous mutations in 22 of the 27 HCC patients. Mutations were detected in known HCC-associated genes (e.g., CTNNB1, TP53, NFE2L2, and ARID1A). No difference in cfDNA concentrations was observed between patients with mutations and those without detectable mutations. HCC patients with higher cfDNA concentrations or higher number of mutations had a shorter overall survival (p < 0.001 and p = 0.045). Nonsynonymous mutations were also identified in 17 of the 51 subjects with advanced liver fibrosis. KMT2C was the most commonly mutated gene. Nine genes were mutated in both subjects with advanced fibrosis and HCC patients. Again, no significant difference in cfDNA concentrations was observed between subjects with mutations and those without detectable mutations. Furthermore, higher cfDNA concentrations and higher number of mutations correlated with a death outcome in subjects with advanced fibrosis. In conclusion, cfDNA features are promising non-invasive markers for HCC risk prediction and overall survival.

Enhancer reprogramming in PRC2-deficient malignant peripheral nerve sheath tumors induces a targetable de-differentiated state.

Malignant peripheral nerve sheath tumors (MPNSTs) are soft tissue sarcomas that frequently harbor genetic alterations in polycomb repressor complex 2 (PRC2) components-SUZ12 and EED. Here, we show that PRC2 loss confers a dedifferentiated early neural-crest phenotype which is exclusive to PRC2-mutant MPNSTs and not a feature of neurofibromas. Neural crest phenotype in PRC2 mutant MPNSTs was validated via cross-species comparative analysis using spontaneous and transgenic MPNST models. Systematic chromatin state profiling of the MPNST cells showed extensive epigenomic reprogramming or chromatin states associated with PRC2 loss and identified gains of active enhancer states/super-enhancers on early neural crest regulators in PRC2-mutant conditions around genomic loci that harbored repressed/poised states in PRC2-WT MPNST cells. Consistently, inverse correlation between H3K27me3 loss and H3K27Ac gain was noted in MPNSTs. Epigenetic editing experiments established functional roles for enhancer gains on DLX5-a key regulator of neural crest phenotype. Consistently, blockade of enhancer activity by bromodomain inhibitors specifically suppressed this neural crest phenotype and tumor burden in PRC2-mutant PDXs. Together, these findings reveal accumulation of dedifferentiated neural crest like state in PRC2-mutant MPNSTs that can be targeted by enhancer blockade.

Multi-omic molecular profiling reveals potentially targetable abnormalities shared across multiple histologies of brain metastasis.

The deadly complication of brain metastasis (BM) is largely confined to a relatively narrow cross-section of systemic malignancies, suggesting a fundamental role for biological mechanisms shared across commonly brain metastatic tumor types. To identify and characterize such mechanisms, we performed genomic, transcriptional, and proteomic profiling using whole-exome sequencing, mRNA-seq, and reverse-phase protein array analysis in a cohort of the lung, breast, and renal cell carcinomas consisting of BM and patient-matched primary or extracranial metastatic tissues. While no specific genomic alterations were associated with BM, correlations with impaired cellular immunity, upregulated oxidative phosphorylation (OXPHOS), and canonical oncogenic signaling pathways including phosphoinositide 3-kinase (PI3K) signaling, were apparent across multiple tumor histologies. Multiplexed immunofluorescence analysis confirmed significant T cell depletion in BM, indicative of a fundamentally altered immune microenvironment. Moreover, functional studies using in vitro and in vivo modeling demonstrated heightened oxidative metabolism in BM along with sensitivity to OXPHOS inhibition in murine BM models and brain metastatic derivatives relative to isogenic parentals. These findings demonstrate that pathophysiological rewiring of oncogenic signaling, cellular metabolism, and immune microenvironment broadly characterizes BM. Further clarification of this biology will likely reveal promising targets for therapeutic development against BM arising from a broad variety of systemic cancers.

A Cryptic BCR-PDGFRB Fusion Resulting in a Chronic Myeloid Neoplasm With Monocytosis and Eosinophilia: A Novel Finding With Treatment Implications.

RNA-seq was used to identify the partner gene and confirm the presence of a BCR-PDGFRB fusion. Identification of this fusion product resulted in successful treatment and long-term remission of this myeloid neoplasm. Based on our results, we suggest that despite current WHO recommendations, screening for PDGFRB rearrangement in cases of leukocytosis with eosinophilia and no other etiologic explanation is necessary, even if the karyotype is normal.

Telomere dysfunction activates YAP1 to drive tissue inflammation.

Germline telomere maintenance defects are associated with an increased incidence of inflammatory diseases in humans, yet whether and how telomere dysfunction causes inflammation are not known. Here, we show that telomere dysfunction drives pATM/c-ABL-mediated activation of the YAP1 transcription factor, up-regulating the major pro-inflammatory factor, pro-IL-18. The colonic microbiome stimulates cytosolic receptors activating caspase-1 which cleaves pro-IL-18 into mature IL-18, leading to recruitment of interferon (IFN)-γ-secreting T cells and intestinal inflammation. Correspondingly, patients with germline telomere maintenance defects exhibit DNA damage (γH2AX) signaling together with elevated YAP1 and IL-18 expression. In mice with telomere dysfunction, telomerase reactivation in the intestinal epithelium or pharmacological inhibition of ATM, YAP1, or caspase-1 as well as antibiotic treatment, dramatically reduces IL-18 and intestinal inflammation. Thus, telomere dysfunction-induced activation of the ATM-YAP1-pro-IL-18 pathway in epithelium is a key instigator of tissue inflammation.

Blocking immunosuppressive neutrophils deters pY696-EZH2-driven brain metastases.

The functions of immune cells in brain metastases are unclear because the brain has traditionally been considered "immune privileged." However, we found that a subgroup of immunosuppressive neutrophils is recruited into the brain, enabling brain metastasis development. In brain metastatic cells, enhancer of zeste homolog 2 (EZH2) is highly expressed and phosphorylated at tyrosine-696 (pY696)-EZH2 by nuclear-localized Src tyrosine kinase. Phosphorylation of EZH2 at Y696 changes its binding preference from histone H3 to RNA polymerase II, which consequently switches EZH2's function from a methyltransferase to a transcription factor that increases c-JUN expression. c-Jun up-regulates protumorigenic inflammatory cytokines, including granulocyte colony-stimulating factor (G-CSF), which recruits Arg1+- and PD-L1+ immunosuppressive neutrophils into the brain to drive metastasis outgrowth. G-CSF-blocking antibodies or immune checkpoint blockade therapies combined with Src inhibitors impeded brain metastasis in multiple mouse models. These findings indicate that pY696-EZH2 can function as a methyltransferase-independent transcription factor to facilitate the brain infiltration of immunosuppressive neutrophils, which could be clinically targeted for brain metastasis treatment.

Spatially resolved analyses link genomic and immune diversity and reveal unfavorable neutrophil activation in melanoma.

Complex tumor microenvironmental (TME) features influence the outcome of cancer immunotherapy (IO). Here we perform immunogenomic analyses on 67 intratumor sub-regions of a PD-1 inhibitor-resistant melanoma tumor and 2 additional metastases arising over 8 years, to characterize TME interactions. We identify spatially distinct evolution of copy number alterations influencing local immune composition. Sub-regions with chromosome 7 gain display a relative lack of leukocyte infiltrate but evidence of neutrophil activation, recapitulated in The Cancer Genome Atlas (TCGA) samples, and associated with lack of response to IO across three clinical cohorts. Whether neutrophil activation represents cause or consequence of local tumor necrosis requires further study. Analyses of T-cell clonotypes reveal the presence of recurrent priming events manifesting in a dominant T-cell clonotype over many years. Our findings highlight the links between marked levels of genomic and immune heterogeneity within the physical space of a tumor, with implications for biomarker evaluation and immunotherapy response.

Molecular Analysis of Clinically Defined Subsets of High-Grade Serous Ovarian Cancer.

The diversity and heterogeneity within high-grade serous ovarian cancer (HGSC), which is the most lethal gynecologic malignancy, is not well understood. Here, we perform comprehensive multi-platform omics analyses, including integrated analysis, and immune monitoring on primary and metastatic sites from highly clinically annotated HGSC samples based on a laparoscopic triage algorithm from patients who underwent complete gross resection (R0) or received neoadjuvant chemotherapy (NACT) with excellent or poor response. We identify significant distinct molecular abnormalities and cellular changes and immune cell repertoire alterations between the groups, including a higher rate of NF1 copy number loss, and reduced chromothripsis-like patterns, higher levels of strong-binding neoantigens, and a higher number of infiltrated T cells in the R0 versus the NACT groups.

Comprehensive T cell repertoire characterization of non-small cell lung cancer.

Immunotherapy targeting T cells is increasingly utilized to treat solid tumors including non-small cell lung cancer (NSCLC). This requires a better understanding of the T cells in the lungs of patients with NSCLC. Here, we report T cell repertoire analysis in a cohort of 236 early-stage NSCLC patients. T cell repertoire attributes are associated with clinicopathologic features, mutational and immune landscape. A considerable proportion of the most prevalent T cells in tumors are also prevalent in the uninvolved tumor-adjacent lungs and appear specific to shared background mutations or viral infections. Patients with higher T cell repertoire homology between the tumor and uninvolved tumor-adjacent lung, suggesting a less tumor-focused T cell response, exhibit inferior survival. These findings indicate that a concise understanding of antigens and T cells in NSCLC is needed to improve therapeutic efficacy and reduce toxicity with immunotherapy, particularly adoptive T cell therapy.

Multiplex profiling of peritoneal metastases from gastric adenocarcinoma identified novel targets and molecular subtypes that predict treatment response.

OBJECTIVE: Peritoneal carcinomatosis (PC) occurs frequently in patients with gastric adenocarcinoma (GAC) and confers a poor prognosis. Multiplex profiling of primary GACs has been insightful but the underpinnings of PC's development/progression remain largely unknown. We characterised exome/transcriptome/immune landscapes of PC cells from patients with GAC aiming to identify novel therapeutic targets. DESIGN: We performed whole-exome sequencing (WES) and whole transcriptome sequencing (RNA-seq) on 44 PC specimens (43 patients with PC) including an integrative analysis of WES, RNA-seq, immune profile, clinical and pathological phenotypes to dissect the molecular pathogenesis, identifying actionable targets and/or biomarkers and comparison with TCGA primary GACs. RESULTS: We identified distinct alterations in PC versus primary GACs, such as more frequent CDH1 and TAF1 mutations, 6q loss and chr19 gain. Alterations associated with aggressive PC phenotypes emerged with increased mutations in TP53, CDH1, TAF1 and KMT2C, higher level of 'clock-like' mutational signature, increase in whole-genome doublings, chromosomal instability (particularly, copy number losses), reprogrammed microenvironment, enriched cell cycle pathways, MYC activation and impaired immune response. Integrated analysis identified two main molecular subtypes: 'mesenchymal-like' and 'epithelial-like' with discriminating response to chemotherapy (31% vs 71%). Patients with the less responsive 'mesenchymal-like' subtype had high expression of immune checkpoint T-Cell Immunoglobulin And Mucin Domain-Containing Protein 3 (TIM-3), its ligand galectin-9, V-domain Ig suppressor of T cell activation (VISTA) and transforming growth factor-ß as potential therapeutic immune targets. CONCLUSIONS: We have uncovered the unique mutational landscape, copy number alteration and gene expression profile of PC cells and defined PC molecular subtypes, which correlated with PC therapy resistance/response. Novel targets and immune checkpoint proteins have been identified with a potential to be translated into clinics.

A pilot study of pembrolizumab in smoldering myeloma: report of the clinical, immune, and genomic analysis.

Multiple myeloma is, in most patients, an incurable cancer. Its precursors can be identified with routine tests setting the stage for early intervention to prevent active myeloma. We investigated the efficacy and safety of pembrolizumab, an antiprogrammed cell death 1 antibody, in smoldering myeloma patients with intermediate/high risk of progression to symptomatic myeloma. Thirteen patients were treated with a median number of 8 cycles. One patient achieved a stringent complete response with bone marrow next-generation sequencing negativity at 10-4 that is ongoing at 27 months (8%); 11 had stable disease (85%), and 1 progressed (8%). Three patients discontinued therapy due to immune-related adverse events: 2 with transaminitis and 1 due to tubulointerstitial nephritis. Immune profiling of bone marrow samples at baseline showed markers associated with a preexisting immune response in the responder compared with nonresponders and features of increased T-cell exhaustion in nonresponders. Consistent with this, transcriptome sequencing of bone marrow samples at baseline revealed an increased interferon-γ signature in the responder compared with the nonresponders. In summary, our results suggest that smoldering myeloma may be immunogenic in a subset of patients, and therapies that enhance antitumor T-cell responses may be effective in preventing its progression. This trial was registered at www.clinicaltrials.gov as #NCT02603887.

PRDM16s transforms megakaryocyte-erythroid progenitors into myeloid leukemia-initiating cells.

Oncogenic mutations confer on cells the ability to propagate indefinitely, but whether oncogenes alter the cell fate of these cells is unknown. Here, we show that the transcriptional regulator PRDM16s causes oncogenic fate conversion by transforming cells fated to form platelets and erythrocytes into myeloid leukemia stem cells (LSCs). Prdm16s expression in megakaryocyte-erythroid progenitors (MEPs), which normally lack the potential to generate granulomonocytic cells, caused AML by converting MEPs into LSCs. Prdm16s blocked megakaryocytic/erythroid potential by interacting with super enhancers and activating myeloid master regulators, including PU.1. A CRISPR dropout screen confirmed that PU.1 is required for Prdm16s-induced leukemia. Ablating PU.1 attenuated leukemogenesis and reinstated the megakaryocytic/erythroid potential of leukemic MEPs in mouse models and human AML with PRDM16 rearrangement. Thus, oncogenic PRDM16 s expression gives MEPs an LSC fate by activating myeloid gene regulatory networks.

Multi-region exome sequencing reveals genomic evolution from preneoplasia to lung adenocarcinoma.

There has been a dramatic increase in the detection of lung nodules, many of which are preneoplasia atypical adenomatous hyperplasia (AAH), adenocarcinoma in situ (AIS), minimally invasive adenocarcinoma (MIA) or invasive adenocarcinoma (ADC). The molecular landscape and the evolutionary trajectory of lung preneoplasia have not been well defined. Here, we perform multi-region exome sequencing of 116 resected lung nodules including AAH (n = 22), AIS (n = 27), MIA (n = 54) and synchronous ADC (n = 13). Comparing AAH to AIS, MIA and ADC, we observe progressive genomic evolution at the single nucleotide level and demarcated evolution at the chromosomal level supporting the early lung carcinogenesis model from AAH to AIS, MIA and ADC. Subclonal analyses reveal a higher proportion of clonal mutations in AIS/MIA/ADC than AAH suggesting neoplastic transformation of lung preneoplasia is predominantly associated with a selective sweep of unfit subclones. Analysis of multifocal pulmonary nodules from the same patients reveal evidence of convergent evolution.

Metabolic reprogramming toward oxidative phosphorylation identifies a therapeutic target for mantle cell lymphoma.

Metabolic reprogramming is linked to cancer cell growth and proliferation, metastasis, and therapeutic resistance in a multitude of cancers. Targeting dysregulated metabolic pathways to overcome resistance, an urgent clinical need in all relapsed/refractory cancers, remains difficult. Through genomic analyses of clinical specimens, we show that metabolic reprogramming toward oxidative phosphorylation (OXPHOS) and glutaminolysis is associated with therapeutic resistance to the Bruton's tyrosine kinase inhibitor ibrutinib in mantle cell lymphoma (MCL), a B cell lymphoma subtype with poor clinical outcomes. Inhibition of OXPHOS with a clinically applicable small molecule, IACS-010759, which targets complex I of the mitochondrial electron transport chain, results in marked growth inhibition in vitro and in vivo in ibrutinib-resistant patient-derived cancer models. This work suggests that targeting metabolic pathways to subvert therapeutic resistance is a clinically viable approach to treat highly refractory malignancies.

Combined Analysis of Antigen Presentation and T-cell Recognition Reveals Restricted Immune Responses in Melanoma.

The quest for tumor-associated antigens (TAA) and neoantigens is a major focus of cancer immunotherapy. Here, we combine a neoantigen prediction pipeline and human leukocyte antigen (HLA) peptidomics to identify TAAs and neoantigens in 16 tumors derived from seven patients with melanoma and characterize their interactions with their tumor-infiltrating lymphocytes (TIL). Our investigation of the antigenic and T-cell landscapes encompassing the TAA and neoantigen signatures, their immune reactivity, and their corresponding T-cell identities provides the first comprehensive analysis of cancer cell T-cell cosignatures, allowing us to discover remarkable antigenic and TIL similarities between metastases from the same patient. Furthermore, we reveal that two neoantigen-specific clonotypes killed 90% of autologous melanoma cells, both in vitro and in vivo, showing that a limited set of neoantigen-specific T cells may play a central role in melanoma tumor rejection. Our findings indicate that combining HLA peptidomics with neoantigen predictions allows robust identification of targetable neoantigens, which could successfully guide personalized cancer immunotherapies.Significance: As neoantigen targeting is becoming more established as a powerful therapeutic approach, investigating these molecules has taken center stage. Here, we show that a limited set of neoantigen-specific T cells mediates tumor rejection, suggesting that identifying just a few antigens and their corresponding T-cell clones could guide personalized immunotherapy. Cancer Discov; 8(11); 1366-75. ©2018 AACR. This article is highlighted in the In This Issue feature, p. 1333.

Variants with a low allele frequency detected in genomic DNA affect the accuracy of mutation detection in cell-free DNA by next-generation sequencing.

BACKGROUND: Next-generation sequencing of cell-free DNA (cfDNA) has been shown to be a useful noninvasive test for detecting mutations in solid tumors. METHODS: Targeted gene sequencing was performed with a panel of 263 cancer-related genes for cfDNA and genomic DNA of peripheral blood mononuclear cells (PBMCs) obtained from presurgical specimens of 6 lung cancer patients, and mutation calls in these samples were compared with those of primary tumors and corresponding patient-derived xenografts (PDXs). RESULTS: Approximately 67% of the mutations detected in the tumor samples (primary tumors and/or PDXs) were also detected in genomic DNA from PBMCs as background mutations. These background mutations consisted of germline polymorphisms and a group of mutations with low allele frequencies, mostly <10%. These variants with a low allele frequency were repeatedly detected in all types of samples from the same patients and at similarly low allele frequency levels in PBMCs from different patients; this indicated that their detection might be derived from common causes, such as homologous sequences in the human genome. Allele frequencies of mutations detected in both primary tumors and cfDNA showed 2 patterns: 1) low allele frequencies (approximately 1%-10%) in cfDNA but high allele frequencies (usually >10% or >3-fold increase) in primary tumors and further enrichment in PDXs and 2) similar allele frequencies across samples. CONCLUSIONS: Because only a small fraction of total cfDNA might be derived from tumor cells, only mutations with the first allele frequency pattern may be regarded as tumor-specific mutations in cfDNA. Effective filtering of background mutations will be required to improve the accuracy of mutation calls in cfDNA. Cancer 2018;124:1061-9. © 2017 American Cancer Society.