The ectonucleotidase CD39 identifies tumor-reactive CD8+ T cells predictive of immune checkpoint blockade efficacy in human lung cancer.

Improved identification of anti-tumor T cells is needed to advance cancer immunotherapies. CD39 expression is a promising surrogate of tumor-reactive CD8+ T cells. Here, we comprehensively profiled CD39 expression in human lung cancer. CD39 expression enriched for CD8+ T cells with features of exhaustion, tumor reactivity, and clonal expansion. Flow cytometry of 440 lung cancer biospecimens revealed weak association between CD39+ CD8+ T cells and tumoral features, such as programmed death-ligand 1 (PD-L1), tumor mutation burden, and driver mutations. Immune checkpoint blockade (ICB), but not cytotoxic chemotherapy, increased intratumoral CD39+ CD8+ T cells. Higher baseline frequency of CD39+ CD8+ T cells conferred improved clinical outcomes from ICB therapy. Furthermore, a gene signature of CD39+ CD8+ T cells predicted benefit from ICB, but not chemotherapy, in a phase III clinical trial of non-small cell lung cancer. These findings highlight CD39 as a proxy of tumor-reactive CD8+ T cells in human lung cancer.

Phase and context shape the function of composite oncogenic mutations.

Cancers develop as a result of driver mutations1,2 that lead to clonal outgrowth and the evolution of disease3,4. The discovery and functional characterization of individual driver mutations are central aims of cancer research, and have elucidated myriad phenotypes5 and therapeutic vulnerabilities6. However, the serial genetic evolution of mutant cancer genes7,8 and the allelic context in which they arise is poorly understood in both common and rare cancer genes and tumour types. Here we find that nearly one in four human tumours contains a composite mutation of a cancer-associated gene, defined as two or more nonsynonymous somatic mutations in the same gene and tumour. Composite mutations are enriched in specific genes, have an elevated rate of use of less-common hotspot mutations acquired in a chronology driven in part by oncogenic fitness, and arise in an allelic configuration that reflects context-specific selective pressures. cis-acting composite mutations are hypermorphic in some genes in which dosage effects predominate (such as TERT), whereas they lead to selection of function in other genes (such as TP53). Collectively, composite mutations are driver alterations that arise from context- and allele-specific selective pressures that are dependent in part on gene and mutation function, and which lead to complex-often neomorphic-functions of biological and therapeutic importance.

Author Correction: Tumour lineage shapes BRCA-mediated phenotypes.

An Amendment to this paper has been published and can be accessed via a link at the top of the paper.

Tumour lineage shapes BRCA-mediated phenotypes.

Mutations in BRCA1 and BRCA2 predispose individuals to certain cancers1-3, and disease-specific screening and preventative strategies have reduced cancer mortality in affected patients4,5. These classical tumour-suppressor genes have tumorigenic effects associated with somatic biallelic inactivation, although haploinsufficiency may also promote the formation and progression of tumours6,7. Moreover, BRCA1/2-mutant tumours are often deficient in the repair of double-stranded DNA breaks by homologous recombination8-13, and consequently exhibit increased therapeutic sensitivity to platinum-containing therapy and inhibitors of poly-(ADP-ribose)-polymerase (PARP)14,15. However, the phenotypic and therapeutic relevance of mutations in BRCA1 or BRCA2 remains poorly defined in most cancer types. Here we show that in the 2.7% and 1.8% of patients with advanced-stage cancer and germline pathogenic or somatic loss-of-function alterations in BRCA1/2, respectively, selective pressure for biallelic inactivation, zygosity-dependent phenotype penetrance, and sensitivity to PARP inhibition were observed only in tumour types associated with increased heritable cancer risk in BRCA1/2 carriers (BRCA-associated cancer types). Conversely, among patients with non-BRCA-associated cancer types, most carriers of these BRCA1/2 mutation types had evidence for tumour pathogenesis that was independent of mutant BRCA1/2. Overall, mutant BRCA is an indispensable founding event for some tumours, but in a considerable proportion of other cancers, it appears to be biologically neutral-a difference predominantly conditioned by tumour lineage-with implications for disease pathogenesis, screening, design of clinical trials and therapeutic decision-making.

Genome-wide loss of heterozygosity predicts aggressive, treatment-refractory behavior in pituitary neuroendocrine tumors.

Pituitary neuroendocrine tumors (PitNETs) exhibiting aggressive, treatment-refractory behavior are the rare subset that progress after surgery, conventional medical therapies, and an initial course of radiation and are characterized by unrelenting growth and/or metastatic dissemination. Two groups of patients with PitNETs were sequenced: a prospective group of patients (n = 66) who consented to sequencing prior to surgery and a retrospective group (n = 26) comprised of aggressive/higher risk PitNETs. A higher mutational burden and fraction of loss of heterozygosity (LOH) was found in the aggressive, treatment-refractory PitNETs compared to the benign tumors (p = 1.3 × 10-10 and p = 8.5 × 10-9, respectively). Within the corticotroph lineage, a characteristic pattern of recurrent chromosomal LOH in 12 specific chromosomes was associated with treatment-refractoriness (occurring in 11 of 14 treatment-refractory versus 1 of 14 benign corticotroph PitNETs, p = 1.7 × 10-4). Across the cohort, a higher fraction of LOH was identified in tumors with TP53 mutations (p = 3.3 × 10-8). A machine learning approach identified loss of heterozygosity as the most predictive variable for aggressive, treatment-refractory behavior, outperforming the most common gene-level alteration, TP53, with an accuracy of 0.88 (95% CI: 0.70-0.96). Aggressive, treatment-refractory PitNETs are characterized by significant aneuploidy due to widespread chromosomal LOH, most prominently in the corticotroph tumors. This LOH predicts treatment-refractoriness with high accuracy and represents a novel biomarker for this poorly defined PitNET category.

Alterations in driver genes are predictive of survival in patients with resected pancreatic ductal adenocarcinoma.

BACKGROUND: KRAS, TP53, CDKN2A, and SMAD4 are established driver genes in pancreatic ductal adenocarcinoma (PDAC). This study was aimed at determining whether the mutational status of driver genes and those involved in DNA repair pathways are associated with clinical outcomes for individuals who undergo resection. METHODS: Eligible individuals were those who underwent resection of PDAC and consented to targeted sequencing of their primary tumor via Memorial Sloan Kettering-Integrated Mutation Profiling of Actionable Cancer Targets (MSK-IMPACT). Genomic alterations were determined on the basis of MSK-IMPACT results from formalin-fixed, paraffin-embedded samples. Associations between genomic alterations and clinical outcomes were assessed. RESULTS: Targeted sequencing was performed on 283 primary tumors resected between 2004 and 2017. The median follow-up was 23 months among survivors. Alterations in KRAS and TP53 were associated with worse overall survival (OS) in comparison to wild type (median for KRAS, 38.8 months [95% CI, 33.0-45.5 months] vs 91.0 months [95% CI, 34.8 months to not available (NA)]; P = .043; median for TP53, 37.4 months [95% CI, 32.1-42.8 months] vs 65.0 months [95% CI, 33.0 months to NA]; P = .035). KRAS G12D mutations were associated with worse OS (median, 31.6 months [95% CI, 25.3-45.5 months] vs 39.2 months [95% CI, 37.4-75.2 months]; P = .012). TP53 truncating mutations (median, 39.6 months [95% CI, 32.4-75.2 months] vs 33.9 months [95% CI, 24.0-39.0 months]; P = .020) and those associated with loss of heterozygosity (median, 26.6 months [95% CI, 21.6-44.2 months] vs 39.2 months [95% CI, 34.5-49.1 months]; P = .048) had decreased OS. TP53 alterations were independently associated with OS in a multivariate analysis (hazard ratio, 1.54; 95% CI, 1.01-2.33; P = .042). Individuals with germline alterations in homologous recombination deficiency (HRD) genes had improved OS in comparison with those without them (median, not reached vs 37.0 months; 95% CI, 33.0-49.8 months; P = .035). CONCLUSIONS: In patients with resected PDAC, genomic alterations in KRAS and TP53 are associated with worse outcomes, whereas alterations in HRD genes are associated with a favorable prognosis. Further studies are needed to better define these alterations as biomarkers in resected PDAC.

Environmental perturbations lead to extensive directional shifts in RNA processing.

Environmental perturbations have large effects on both organismal and cellular traits, including gene expression, but the extent to which the environment affects RNA processing remains largely uncharacterized. Recent studies have identified a large number of genetic variants associated with variation in RNA processing that also have an important role in complex traits; yet we do not know in which contexts the different underlying isoforms are used. Here, we comprehensively characterized changes in RNA processing events across 89 environments in five human cell types and identified 15,300 event shifts (FDR = 15%) comprised of eight event types in over 4,000 genes. Many of these changes occur consistently in the same direction across conditions, indicative of global regulation by trans factors. Accordingly, we demonstrate that environmental modulation of splicing factor binding predicts shifts in intron retention, and that binding of transcription factors predicts shifts in alternative first exon (AFE) usage in response to specific treatments. We validated the mechanism hypothesized for AFE in two independent datasets. Using ATAC-seq, we found altered binding of 64 factors in response to selenium at sites of AFE shift, including ELF2 and other factors in the ETS family. We also performed AFE QTL mapping in 373 individuals and found an enrichment for SNPs predicted to disrupt binding of the ELF2 factor. Together, these results demonstrate that RNA processing is dramatically changed in response to environmental perturbations through specific mechanisms regulated by trans factors.

High-throughput allele-specific expression across 250 environmental conditions.

Gene-by-environment (GxE) interactions determine common disease risk factors and biomedically relevant complex traits. However, quantifying how the environment modulates genetic effects on human quantitative phenotypes presents unique challenges. Environmental covariates are complex and difficult to measure and control at the organismal level, as found in GWAS and epidemiological studies. An alternative approach focuses on the cellular environment using in vitro treatments as a proxy for the organismal environment. These cellular environments simplify the organism-level environmental exposures to provide a tractable influence on subcellular phenotypes, such as gene expression. Expression quantitative trait loci (eQTL) mapping studies identified GxE interactions in response to drug treatment and pathogen exposure. However, eQTL mapping approaches are infeasible for large-scale analysis of multiple cellular environments. Recently, allele-specific expression (ASE) analysis emerged as a powerful tool to identify GxE interactions in gene expression patterns by exploiting naturally occurring environmental exposures. Here we characterized genetic effects on the transcriptional response to 50 treatments in five cell types. We discovered 1455 genes with ASE (FDR < 10%) and 215 genes with GxE interactions. We demonstrated a major role for GxE interactions in complex traits. Genes with a transcriptional response to environmental perturbations showed sevenfold higher odds of being found in GWAS. Additionally, 105 genes that indicated GxE interactions (49%) were identified by GWAS as associated with complex traits. Examples include GIPR-caffeine interaction and obesity and include LAMP3-selenium interaction and Parkinson disease. Our results demonstrate that comprehensive catalogs of GxE interactions are indispensable to thoroughly annotate genes and bridge epidemiological and genome-wide association studies.

RNA-DNA differences in human mitochondria restore ancestral form of 16S ribosomal RNA.

RNA transcripts are generally identical to the underlying DNA sequences. Nevertheless, RNA-DNA differences (RDDs) were found in the nuclear human genome and in plants and animals but not in human mitochondria. Here, by deep sequencing of human mitochondrial DNA (mtDNA) and RNA, we identified three RDD sites at mtDNA positions 295 (C-to-U), 13710 (A-to-U, A-to-G), and 2617 (A-to-U, A-to-G). Position 2617, within the 16S rRNA, harbored the most prevalent RDDs (>30% A-to-U and ∼15% A-to-G of the reads in all tested samples). The 2617 RDDs appeared already at the precursor polycistrone mitochondrial transcript. By using traditional Sanger sequencing, we identified the A-to-U RDD in six different cell lines and representative primates (Gorilla gorilla, Pongo pigmaeus, and Macaca mulatta), suggesting conservation of the mechanism generating such RDD. Phylogenetic analysis of more than 1700 vertebrate mtDNA sequences supported a thymine as the primate ancestral allele at position 2617, suggesting that the 2617 RDD recapitulates the ancestral 16S rRNA. Modeling U or G (the RDDs) at position 2617 stabilized the large ribosomal subunit structure in contrast to destabilization by an A (the pre-RDDs). Hence, these mitochondrial RDDs are likely functional.

Interactions between glucocorticoid treatment and cis-regulatory polymorphisms contribute to cellular response phenotypes.

Glucocorticoids (GCs) mediate physiological responses to environmental stress and are commonly used as pharmaceuticals. GCs act primarily through the GC receptor (GR, a transcription factor). Despite their clear biomedical importance, little is known about the genetic architecture of variation in GC response. Here we provide an initial assessment of variability in the cellular response to GC treatment by profiling gene expression and protein secretion in 114 EBV-transformed B lymphocytes of African and European ancestry. We found that genetic variation affects the response of nearby genes and exhibits distinctive patterns of genotype-treatment interactions, with genotypic effects evident in either only GC-treated or only control-treated conditions. Using a novel statistical framework, we identified interactions that influence the expression of 26 genes known to play central roles in GC-related pathways (e.g. NQO1, AIRE, and SGK1) and that influence the secretion of IL6.

Convergent evolution of BRCA2 reversion mutations under therapeutic pressure by PARP inhibition and platinum chemotherapy.

Reversion mutations that restore wild-type function of the BRCA gene have been described as a key mechanism of resistance to Poly(ADP-ribose) polymerase (PARP) inhibitor therapy in BRCA-associated cancers. Here, we report a case of a patient with metastatic castration-resistant prostate cancer (mCRPC) with a germline BRCA2 mutation who developed acquired resistance to PARP inhibition. Extensive genomic interrogation of cell-free DNA (cfDNA) and tissue at baseline, post-progression, and postmortem revealed ten unique BRCA2 reversion mutations across ten sites. While several of the reversion mutations were private to a specific site, nine out of ten tumors contained at least one mutation, suggesting a powerful clonal selection for reversion mutations in the presence of therapeutic pressure by PARP inhibition. Variable cfDNA shed was seen across tumor sites, emphasizing a potential shortcoming of cfDNA monitoring for PARPi resistance. This report provides a genomic portrait of the temporal and spatial heterogeneity of prostate cancer under the selective pressure of a PARP inhibition and exposes limitations in the current strategies for detection of reversion mutations.

Therapy-Induced Senescence Contributes to the Efficacy of Abemaciclib in Patients with Dedifferentiated Liposarcoma.

PURPOSE: We conducted research on CDK4/6 inhibitors (CDK4/6i) simultaneously in the preclinical and clinical spaces to gain a deeper understanding of how senescence influences tumor growth in humans. PATIENTS AND METHODS: We coordinated a first-in-kind phase II clinical trial of the CDK4/6i abemaciclib for patients with progressive dedifferentiated liposarcoma (DDLS) with cellular studies interrogating the molecular basis of geroconversion. RESULTS: Thirty patients with progressing DDLS enrolled and were treated with 200 mg of abemaciclib twice daily. The median progression-free survival was 33 weeks at the time of the data lock, with 23 of 30 progression-free at 12 weeks (76.7%, two-sided 95% CI, 57.7%-90.1%). No new safety signals were identified. Concurrent preclinical work in liposarcoma cell lines identified ANGPTL4 as a necessary late regulator of geroconversion, the pathway from reversible cell-cycle exit to a stably arrested inflammation-provoking senescent cell. Using this insight, we were able to identify patients in which abemaciclib induced tumor cell senescence. Senescence correlated with increased leukocyte infiltration, primarily CD4-positive cells, within a month of therapy. However, those individuals with both senescence and increased TILs were also more likely to acquire resistance later in therapy. These suggest that combining senolytics with abemaciclib in a subset of patients may improve the duration of response. CONCLUSIONS: Abemaciclib was well tolerated and showed promising activity in DDLS. The discovery of ANGPTL4 as a late regulator of geroconversion helped to define how CDK4/6i-induced cellular senescence modulates the immune tumor microenvironment and contributes to both positive and negative clinical outcomes. See related commentary by Weiss et al., p. 649.

Clinical Benefit From Immunotherapy in Patients With SCLC Is Associated With Tumor Capacity for Antigen Presentation.

INTRODUCTION: A small percentage of patients with SCLC experience durable responses to immune checkpoint blockade (ICB). Defining determinants of immune response may nominate strategies to broaden the efficacy of immunotherapy in patients with SCLC. Prior studies have been limited by small numbers or concomitant chemotherapy administration. METHODS: CheckMate 032, a multicenter, open-label, phase 1/2 trial evaluating nivolumab alone or with ipilimumab was the largest study of ICB alone in patients with SCLC. We performed comprehensive RNA sequencing of 286 pretreatment SCLC tumor samples, assessing outcome on the basis of defined SCLC subtypes (SCLC-A, -N, -P, and -Y), and expression signatures associated with durable benefit, defined as progression-free survival more than or equal to 6 months. Potential biomarkers were further explored by immunohistochemistry. RESULTS: None of the subtypes were associated with survival. Antigen presentation machinery signature (p = 0.000032) and presence of more than or equal to 1% infiltrating CD8+ T cells by immunohistochemistry (hazard ratio = 0.51, 95% confidence interval: 0.27-0.95) both correlated with survival in patients treated with nivolumab. Pathway enrichment analysis revealed the association between durable benefit from immunotherapy and antigen processing and presentation. Analysis of epigenetic determinants of antigen presentation identified LSD1 gene expression as a correlate of worse survival outcomes for patients treated with either nivolumab or the combination of nivolumab and ipilimumab. CONCLUSIONS: Tumor antigen processing and presentation is a key correlate of ICB efficacy in patients with SCLC. As antigen presentation machinery is frequently epigenetically suppressed in SCLC, this study defines a targetable mechanism by which we might improve clinical benefit of ICB for patients with SCLC.

A Phase II Study of Epacadostat and Pembrolizumab in Patients with Advanced Sarcoma.

PURPOSE: Epacadostat, an indole 2,3 dioxygenase 1 (IDO1) inhibitor, proposed to shift the tumor microenvironment toward an immune-stimulated state, showed early promise in melanoma but has not been studied in sarcoma. This study combined epacadostat with pembrolizumab, which has modest activity in select sarcoma subtypes. PATIENTS AND METHODS: This phase II study enrolled patients with advanced sarcoma into five cohorts including (i) undifferentiated pleomorphic sarcoma (UPS)/myxofibrosarcoma, (ii) liposarcoma (LPS), (iii) leiomyosarcoma (LMS), (iv) vascular sarcoma, including angiosarcoma and epithelioid hemangioendothelioma (EHE), and (v) other subtypes. Patients received epacadostat 100 mg twice daily plus pembrolizumab at 200 mg/dose every 3 weeks. The primary endpoint was best objective response rate (ORR), defined as complete response (CR) and partial response (PR), at 24 weeks by RECIST v.1.1. RESULTS: Thirty patients were enrolled [60% male; median age 54 years (range, 24-78)]. The best ORR at 24 weeks was 3.3% [PR, n = 1 (leiomyosarcoma); two-sided 95% CI, 0.1%-17.2%]. The median PFS was 7.6 weeks (two-sided 95% CI, 6.9-26.7). Treatment was well tolerated. Grade 3 treatment-related adverse events occurred in 23% (n = 7) of patients. In paired pre- and post-treatment tumor samples, no association was found between treatment and PD-L1 or IDO1 tumor expression or IDO-pathway-related gene expression by RNA sequencing. No significant changes in serum tryptophan or kynurenine levels were observed after baseline. CONCLUSIONS: Combination epacadostat and pembrolizumab was well tolerated and showed limited antitumor activity in sarcoma. Correlative analyses suggested that inadequate IDO1 inhibition was achieved.

Multi-institutional analysis of aneuploidy and outcomes to chemoradiation and durvalumab in stage III non-small cell lung cancer.

There is a need to identify predictive biomarkers to guide treatment strategies in stage III non-small cell lung cancer (NSCLCs). In this multi-institutional cohort of 197 patients with stage III NSCLC treated with concurrent chemoradiation (cCRT) and durvalumab consolidation, we identify that low tumor aneuploidy is independently associated with prolonged progression-free survival (HR 0.63; p=0.03) and overall survival (HR 0.50; p=0.03). Tumors with high aneuploidy had a significantly greater incidence of distant metastasis and shorter median distant-metastasis free survival (p=0.04 and p=0.048, respectively), but aneuploidy level did not associate with local-regional outcomes. Multiplexed immunofluorescence analysis in a cohort of NSCLC found increased intratumoral CD8-positive, PD-1-positive cells, double-positive PD-1 CD8 cells, and FOXP3-positive T-cell in low aneuploid tumors. Additionally, in a cohort of 101 patients treated with cCRT alone, tumor aneuploidy did not associate with disease outcomes. These data support the need for upfront treatment intensification strategies in stage III NSCLC patients with high aneuploid tumors and suggest that tumor aneuploidy is a promising predictive biomarker.

Impact of Aneuploidy and Chromosome 9p Loss on Tumor Immune Microenvironment and Immune Checkpoint Inhibitor Efficacy in NSCLC.

INTRODUCTION: Although gene-level copy number alterations have been studied as a potential biomarker of immunotherapy efficacy in NSCLC, the impact of aneuploidy burden and chromosomal arm-level events on immune checkpoint inhibitor (ICI) efficacy in NSCLC is uncertain. METHODS: Patients who received programmed cell death protein 1 or programmed death-ligand 1 (PD-L1) inhibitor at two academic centers were included. Across all 22 chromosomes analyzed, an arm was considered altered if at least 70% of its territory was either gained or deleted. Among nonsquamous NSCLCs which underwent targeted next-generation sequencing, we retrospectively quantified aneuploidy using the adjusted fraction of chromosomal arm alterations (FAA), defined as the number of altered chromosome arms divided by the number of chromosome arms assessed, adjusted for tumor purity. RESULTS: Among 2293 nonsquamous NSCLCs identified, the median FAA increased with more advanced cancer stage and decreased with higher PD-L1 tumor proportion score (TPS) levels (median FAA in TPS < 1%: 0.09, TPS 1%-49%: 0.08, TPS ≥ 50%: 0.05, p < 0.0001). There was a very weak correlation between FAA and tumor mutational burden when taken as continuous variables (R: 0.07, p = 0.0005). A total of 765 advanced nonsquamous NSCLCs with available FAA values were treated with ICIs. With decreasing FAA tertiles, there was a progressive improvement in objective response rate (ORR 15.1% in upper tertile versus 23.2% in middle tertile versus 28.4% in lowest tertile, p = 0.001), median progression-free survival (mPFS 2.5 versus 3.3 versus 4.1 mo, p < 0.0001), and median overall survival (mOS 12.5 versus 13.9 versus 16.4 mo, p = 0.006), respectively. In the arm-level enrichment analysis, chromosome 9p loss (OR = 0.22, Q = 0.0002) and chromosome 1q gain (OR = 0.43, Q = 0.002) were significantly enriched in ICI nonresponders after false discovery rate adjustment. Compared with NSCLCs without chromosome 9p loss (n = 452), those with 9p loss (n = 154) had a lower ORR (28.1% versus 7.8%, p < 0.0001), a shorter mPFS (4.1 versus 2.3 mo, p < 0.0001), and a shorter mOS (18.0 versus 9.6 mo, p < 0.0001) to immunotherapy. In addition, among NSCLCs with high PD-L1 expression (TPS ≥ 50%), chromosome 9p loss was associated with lower ORR (43% versus 6%, p < 0.0001), shorter mPFS (6.4 versus 2.6 mo, p = 0.0006), and shorter mOS (30.2 versus 14.3 mo, p = 0.0008) to immunotherapy compared with NSCLCs without 9p loss. In multivariable analysis, adjusting for key variables including FAA, chromosome 9p loss, but not 1q gain, retained a significant impact on ORR (hazard ratio [HR] = 0.25, p < 0.001), mPFS (HR = 1.49, p = 0.001), and mOS (HR = 1.47, p = 0.003). Multiplexed immunofluorescence and computational deconvolution of RNA sequencing data revealed that tumors with either high FAA levels or chromosome 9p loss had significantly fewer tumor-associated cytotoxic immune cells. CONCLUSIONS: Nonsquamous NSCLCs with high aneuploidy and chromosome 9p loss have a distinct tumor immune microenvironment and less favorable outcomes to ICIs.

Neoplasia risk in patients with Lynch syndrome treated with immune checkpoint blockade.

Metastatic and localized mismatch repair-deficient (dMMR) tumors are exquisitely sensitive to immune checkpoint blockade (ICB). The ability of ICB to prevent dMMR malignant or pre-malignant neoplasia development in patients with Lynch syndrome (LS) is unknown. Of 172 cancer-affected patients with LS who had received ≥1 ICB cycles, 21 (12%) developed subsequent malignancies after ICB exposure, 91% (29/32) of which were dMMR, with median time to development of 21 months (interquartile range, 6-38). Twenty-four of 61 (39%) ICB-treated patients who subsequently underwent surveillance colonoscopy had premalignant polyps. Within matched pre-ICB and post-ICB follow-up periods, the overall rate of tumor development was unchanged; however, on subgroup analysis, a decreased incidence of post-ICB visceral tumors was observed. These data suggest that ICB treatment of LS-associated tumors does not eliminate risk of new neoplasia development, and LS-specific surveillance strategies should continue. These data have implications for immunopreventative strategies and provide insight into the immunobiology of dMMR tumors.

Long-term Follow-up and Patterns of Response, Progression, and Hyperprogression in Patients after PD-1 Blockade in Advanced Sarcoma.

PURPOSE: Programmed cell death protein 1 (PD-1) blockade can mediate objective responses in advanced sarcomas, but their durability has not been established and it is unclear if hyperprogressive disease (HPD) occurs in sarcomas treated with PD-1 inhibitors. EXPERIMENTAL DESIGN: We pooled patients who were treated prospectively with nivolumab or pembrolizumab as monotherapy or with bempegaldesleukin, epacadostat, ipilimumab, or talimogene laherparepvec. We did a new independent assessment for HPD and analyzed clinical, pathologic, and genomic data from baseline tumor biopsies. Our primary endpoint was the incidence of HPD; secondary endpoints were clinical or genomic correlates of response or HPD. RESULTS: We treated 134 patients with advanced sarcoma from 2015 to 2019. Twenty-one patients (16%) had a complete or partial response (CR/PR), and 30% of responses were durable for over 2 years. Forty-eight (36%) patients had stable disease (SD), 45 (34%) had progressive disease without HPD (PD), and 15 (11%) had HPD. Five patients (4%) were not evaluable for HPD. The sarcoma subtypes, sites of metastasis, clinical course, and genomic alterations in patients with PD and HPD were similar, except HPD tumors were smaller at baseline. CONCLUSIONS: In patients with advanced sarcoma, PD-1 blockade can mediate durable responses. HPD occurs in sarcoma at an incidence that is similar to what has been reported in other solid tumors, but patients with HPD were clinically and biologically similar to those who had PD. Further research is required to establish whether HPD is a biologically distinct phenomenon and whether a theoretical risk of HPD should influence patient management.

Clinical, genomic, and transcriptomic correlates of response to immune checkpoint blockade-based therapy in a cohort of patients with angiosarcoma treated at a single center.

BACKGROUND: Angiosarcoma is a histologically and molecularly heterogeneous vascular neoplasm with aggressive clinical behavior. Emerging data suggests that immune checkpoint blockade (ICB) is efficacious against some angiosarcomas, particularly cutaneous angiosarcoma of the head and neck (CHN). METHODS: Patients with histologically confirmed angiosarcoma treated with ICB-based therapy at a comprehensive cancer center were retrospectively identified. Clinical characteristics and the results of targeted exome sequencing, transcriptome sequencing, and immunohistochemistry analyses were examined for correlation with clinical benefit. Durable clinical benefit was defined as a progression-free survival (PFS) of ≥16 weeks. RESULTS: For the 35 patients included in the analyses, median PFS and median overall survival (OS) from the time of first ICB-based treatment were 11.9 (95% CI 7.4 to 31.9) and 42.5 (95% CI 19.6 to 114.2) weeks, respectively. Thirteen patients (37%) had PFS ≥16 weeks. Clinical factors associated with longer PFS and longer OS in multivariate analyses were ICB plus other therapy regimens, CHN disease, and white race. Three of 10 patients with CHN angiosarcoma evaluable for tumor mutational burden (TMB) had a TMB ≥10. Five of six patients with CHN angiosarcoma evaluable for mutational signature analysis had a dominant mutational signature associated with ultraviolet (UV) light. No individual gene or genomic pathway was significantly associated with PFS or OS; neither were TMB or UV signature status. Analyses of whole transcriptomes from nine patient tumor samples found upregulation of angiogenesis, inflammatory response, and KRAS signaling pathways, among others, in patients with PFS ≥16 weeks, as well as higher levels of cytotoxic T cells, dendritic cells, and natural killer cells. Patients with PFS <16 weeks had higher numbers of cancer-associated fibroblasts. Immunohistochemistry findings for 12 patients with baseline samples available suggest that neither PD-L1 expression nor presence of tumor-infiltrating lymphocytes at baseline appears necessary for a response to ICB-based therapy. CONCLUSIONS: ICB-based therapy benefits only a subset of angiosarcoma patients. Patients with CHN angiosarcoma are more likely to have PFS ≥16 weeks, a dominant UV mutational signature, and higher TMB than angiosarcomas arising from other primary sites. However, clinical benefit was seen in other angiosarcomas also and was not restricted to tumors with a high TMB, a dominant UV signature, PD-L1 expression, or presence of tumor infiltrating lymphocytes at baseline.

Pilot study of bempegaldesleukin in combination with nivolumab in patients with metastatic sarcoma.

PD-1 blockade (nivolumab) efficacy remains modest for metastatic sarcoma. In this paper, we present an open-label, non-randomized, non-comparative pilot study of bempegaldesleukin, a CD122-preferential interleukin-2 pathway agonist, with nivolumab in refractory sarcoma at Memorial Sloan Kettering/MD Anderson Cancer Centers (NCT03282344). We report on the primary outcome of objective response rate (ORR) and secondary endpoints of toxicity, clinical benefit, progression-free survival, overall survival, and durations of response/treatment. In 84 patients in 9 histotype cohorts, all patients experienced ≥1 adverse event and treatment-related adverse event; 1 death was possibly treatment-related. ORR was highest in angiosarcoma (3/8) and undifferentiated pleomorphic sarcoma (2/10), meeting predefined endpoints. Results of our exploratory investigation of predictive biomarkers show: CD8 + T cell infiltrates and PD-1 expression correlate with improved ORR; upregulation of immune-related pathways correlate with improved efficacy; Hedgehog pathway expression correlate with resistance. Exploration of this combination in selected sarcomas, and of Hedgehog signaling as a predictive biomarker, warrants further study in larger cohorts.