Malignant pleural mesotheliomas (MPMs) often show CDKN2A and NF2 inactivation, but other highly recurrent mutations have not been described. To identify additional driver genes, we used an integrated genomic analysis of 53 MPM tumor samples to guide a focused sequencing effort that uncovered somatic inactivating mutations in BAP1 in 23% of MPMs. The BAP1 nuclear deubiquitinase is known to target histones (together with ASXL1 as a Polycomb repressor subunit) and the HCF1 transcriptional co-factor, and we show that BAP1 knockdown in MPM cell lines affects E2F and Polycomb target genes. These findings implicate transcriptional deregulation in the pathogenesis of MPM.

Molecularly targeted therapies aim to obstruct cell autonomous programs required for tumor growth. We show that mitogen-activated protein kinase (MAPK) and cyclin-dependent kinase 4/6 inhibitors act in combination to suppress the proliferation of KRAS-mutant lung cancer cells while simultaneously provoking a natural killer (NK) cell surveillance program leading to tumor cell death. The drug combination, but neither agent alone, promotes retinoblastoma (RB) protein-mediated cellular senescence and activation of the immunomodulatory senescence-associated secretory phenotype (SASP). SASP components tumor necrosis factor-α and intercellular adhesion molecule−1 are required for NK cell surveillance of drug-treated tumor cells, which contributes to tumor regressions and prolonged survival in a KRAS-mutant lung cancer mouse model. Therefore, molecularly targeted agents capable of inducing senescence can produce tumor control through non−cell autonomous mechanisms involving NK cell surveillance.

Crucial transitions in cancer-including tumor initiation, local expansion, metastasis, and therapeutic resistance-involve complex interactions between cells within the dynamic tumor ecosystem. Transformative single-cell genomics technologies and spatial multiplex in situ methods now provide an opportunity to interrogate this complexity at unprecedented resolution. The Human Tumor Atlas Network (HTAN), part of the National Cancer Institute (NCI) Cancer Moonshot Initiative, will establish a clinical, experimental, computational, and organizational framework to generate informative and accessible three-dimensional atlases of cancer transitions for a diverse set of tumor types. This effort complements both ongoing efforts to map healthy organs and previous largescale cancer genomics approaches focused on bulk sequencing at a single point in time. Generating single-cell, multiparametric, longitudinal atlases and integrating them with clinical outcomes should help identify novel predictive biomarkers and features as well as therapeutically relevant cell types, cell states, and cellular interactions across transitions. The resulting tumor atlases should have a profound impact on our understanding of cancer biology and have the potential to improve cancer detection, prevention, and therapeutic discovery for better precision-medicine treatments of cancer patients and those at risk for cancer.Cancer forms and progresses through a series of critical transitions-from pre-malignant to malignant states, from locally contained to metastatic disease, and from treatment-responsive to treatment-resistant tumors (Figure 1). Although specifics differ across tumor types and patients, all transitions involve complex dynamic interactions between diverse pre-malignant, malignant, and non-malignant cells (e.g., stroma cells and immune cells), often organized in specific patterns within the tumor

PD-1 blockade unleashes CD8 T cells1, including those specific for mutation-associated neoantigens (MANA), but factors in the tumour microenvironment can inhibit these T cell responses. Single-cell transcriptomics have revealed global T cell dysfunction programs in tumour-infiltrating lymphocytes (TIL). However, the majority of TIL do not recognize tumour antigens2, and little is known about transcriptional programs of MANA-specific TIL. Here, we identify MANA-specific T cell clones using the MANA functional expansion of specific T cells assay3 in neoadjuvant anti-PD-1-treated non-small cell lung cancers (NSCLC). We use their T cell receptors as a ‘barcode’ to track and analyse their transcriptional programs in the tumour microenvironment using coupled single-cell RNA sequencing and T cell receptor sequencing. We find both MANA- and virus-specific clones in TIL, regardless of response, and MANA-, influenza- and Epstein–Barr virus-specific TIL each have unique transcriptional programs. Despite exposure to cognate antigen, MANA-specific TIL express an incompletely activated cytolytic program. MANA-specific CD8 T cells have hallmark transcriptional programs of tissue-resident memory (TRM) cells, but low levels of interleukin-7 receptor (IL-7R) and are functionally less responsive to interleukin-7 (IL-7) compared with influenza-specific TRM cells. Compared with those from responding tumours, MANA-specific clones from non-responding tumours express T cell receptors with markedly lower ligand-dependent signalling, are largely confined to HOBIThigh TRM subsets, and coordinately upregulate checkpoints, killer inhibitory receptors and inhibitors of T cell activation. These findings provide important insights for overcoming resistance to PD-1 blockade.