An atlas of epithelial cell states and plasticity in lung adenocarcinoma.

Understanding the cellular processes that underlie early lung adenocarcinoma (LUAD) development is needed to devise intervention strategies1. Here we studied 246,102 single epithelial cells from 16 early-stage LUADs and 47 matched normal lung samples. Epithelial cells comprised diverse normal and cancer cell states, and diversity among cancer cells was strongly linked to LUAD-specific oncogenic drivers. KRAS mutant cancer cells showed distinct transcriptional features, reduced differentiation and low levels of aneuploidy. Non-malignant areas surrounding human LUAD samples were enriched with alveolar intermediate cells that displayed elevated KRT8 expression (termed KRT8+ alveolar intermediate cells (KACs) here), reduced differentiation, increased plasticity and driver KRAS mutations. Expression profiles of KACs were enriched in lung precancer cells and in LUAD cells and signified poor survival. In mice exposed to tobacco carcinogen, KACs emerged before lung tumours and persisted for months after cessation of carcinogen exposure. Moreover, they acquired Kras mutations and conveyed sensitivity to targeted KRAS inhibition in KAC-enriched organoids derived from alveolar type 2 (AT2) cells. Last, lineage-labelling of AT2 cells or KRT8+ cells following carcinogen exposure showed that KACs are possible intermediates in AT2-to-tumour cell transformation. This study provides new insights into epithelial cell states at the root of LUAD development, and such states could harbour potential targets for prevention or intervention.

Evolution of immune and stromal cell states and ecotypes during gastric adenocarcinoma progression.

Understanding tumor microenvironment (TME) reprogramming in gastric adenocarcinoma (GAC) progression may uncover novel therapeutic targets. Here, we performed single-cell profiling of precancerous lesions, localized and metastatic GACs, identifying alterations in TME cell states and compositions as GAC progresses. Abundant IgA+ plasma cells exist in the premalignant microenvironment, whereas immunosuppressive myeloid and stromal subsets dominate late-stage GACs. We identified six TME ecotypes (EC1-6). EC1 is exclusive to blood, while EC4, EC5, and EC2 are highly enriched in uninvolved tissues, premalignant lesions, and metastases, respectively. EC3 and EC6, two distinct ecotypes in primary GACs, associate with histopathological and genomic characteristics, and survival outcomes. Extensive stromal remodeling occurs in GAC progression. High SDC2 expression in cancer-associated fibroblasts (CAFs) is linked to aggressive phenotypes and poor survival, and SDC2 overexpression in CAFs contributes to tumor growth. Our study provides a high-resolution GAC TME atlas and underscores potential targets for further investigation.

Proteogenomic landscape of gastric adenocarcinoma peritoneal metastases.

Advanced gastric adenocarcinoma (GAC) often leads to peritoneal carcinomatosis (PC) and is associated with very poor outcome. Here we report the comprehensive proteogenomic study of ascites derived cells from a prospective GAC cohort (n = 26 patients with peritoneal carcinomatosis, PC). A total of 16,449 proteins were detected from whole cell extracts (TCEs). Unsupervised hierarchical clustering resulted in three distinct groups that reflected extent of enrichment in tumor cells. Integrated analysis revealed enriched biological pathways and notably, some druggable targets (cancer-testis antigens, kinases, and receptors) that could be exploited to develop effective therapies and/or tumor stratifications. Systematic comparison of expression levels of proteins and mRNAs revealed special expression patterns of key therapeutics target notably high mRNA and low protein expression of HAVCR2 (TIM-3), and low mRNA but high protein expression of cancer-testis antigens CTAGE1 and CTNNA2. These results inform strategies to target GAC vulnerabilities.

Single cell clonotypic and transcriptional evolution of multiple myeloma precursor disease.

Multiple myeloma remains an incurable disease, and the cellular and molecular evolution from precursor conditions, including monoclonal gammopathy of undetermined significance and smoldering multiple myeloma, is incompletely understood. Here, we combine single-cell RNA and B cell receptor sequencing from fifty-two patients with myeloma precursors in comparison with myeloma and normal donors. Our comprehensive analysis reveals early genomic drivers of malignant transformation, distinct transcriptional features, and divergent clonal expansion in hyperdiploid versus non-hyperdiploid samples. Additionally, we observe intra-patient heterogeneity with potential therapeutic implications and identify distinct patterns of evolution from myeloma precursor disease to myeloma. We also demonstrate distinctive characteristics of the microenvironment associated with specific genomic changes in myeloma cells. These findings add to our knowledge about myeloma precursor disease progression, providing valuable insights into patient risk stratification, biomarker discovery, and possible clinical applications.

Pan-cancer T cell atlas links a cellular stress response state to immunotherapy resistance.

Tumor-infiltrating T cells offer a promising avenue for cancer treatment, yet their states remain to be fully characterized. Here we present a single-cell atlas of T cells from 308,048 transcriptomes across 16 cancer types, uncovering previously undescribed T cell states and heterogeneous subpopulations of follicular helper, regulatory and proliferative T cells. We identified a unique stress response state, TSTR, characterized by heat shock gene expression. TSTR cells are detectable in situ in the tumor microenvironment across various cancer types, mostly within lymphocyte aggregates or potential tertiary lymphoid structures in tumor beds or surrounding tumor edges. T cell states/compositions correlated with genomic, pathological and clinical features in 375 patients from 23 cohorts, including 171 patients who received immune checkpoint blockade therapy. We also found significantly upregulated heat shock gene expression in intratumoral CD4/CD8+ cells following immune checkpoint blockade treatment, particularly in nonresponsive tumors, suggesting a potential role of TSTR cells in immunotherapy resistance. Our well-annotated T cell reference maps, web portal and automatic alignment/annotation tool could provide valuable resources for T cell therapy optimization and biomarker discovery.

The Single-Cell Immunogenomic Landscape of B and Plasma Cells in Early-Stage Lung Adenocarcinoma.

Tumor-infiltrating B and plasma cells (TIB) are prevalent in lung adenocarcinoma (LUAD); however, they are poorly characterized. We performed paired single-cell RNA and B-cell receptor (BCR) sequencing of 16 early-stage LUADs and 47 matching multiregion normal tissues. By integrative analysis of ∼50,000 TIBs, we define 12 TIB subsets in the LUAD and adjacent normal ecosystems and demonstrate extensive remodeling of TIBs in LUADs. Memory B cells and plasma cells (PC) were highly enriched in tumor tissues with more differentiated states and increased frequencies of somatic hypermutation. Smokers exhibited markedly elevated PCs and PCs with distinct differentiation trajectories. BCR clonotype diversity increased but clonality decreased in LUADs, smokers, and with increasing pathologic stage. TIBs were mostly localized within CXCL13+ lymphoid aggregates, and immune cell sources of CXCL13 production evolved with LUAD progression and included elevated fractions of CD4 regulatory T cells. This study provides a spatial landscape of TIBs in early-stage LUAD. SIGNIFICANCE: While TIBs are highly enriched in LUADs, they are poorly characterized. This study provides a much-needed understanding of the transcriptional, clonotypic states and phenotypes of TIBs, unraveling their potential roles in the immunopathology of early-stage LUADs and constituting a road map for the development of TIB-targeted immunotherapies for the treatment of this morbid malignancy. This article is highlighted in the In This Issue feature, p. 2483.

TIGIT is the central player in T-cell suppression associated with CAR T-cell relapse in mantle cell lymphoma.

BACKGROUND: Chimeric antigen receptor (CAR) T-cell therapy using brexucabtagene autoleucel (BA) induces remission in many patients with mantle cell lymphoma (MCL), and BA is the only CAR T-cell therapy approved by the FDA for MCL. However, development of relapses to BA is recognized with poor patient outcomes. Multiple CAR T-cell therapies have been approved for other lymphomas and the resistance mechanisms have been investigated. However, the mechanisms underlying BA relapse in MCL have not been investigated and whether any previously reported resistance mechanisms apply to BA-relapsed patients with MCL is unknown. METHODS: To interrogate BA resistance mechanisms in MCL, we performed single-cell RNA sequencing on 39 longitudinally collected samples from 15 BA-treated patients, and multiplex cytokine profiling on 80 serial samples from 20 patients. RESULTS: We demonstrate that after BA relapse, the proportion of T cells, especially cytotoxic T cells (CTLs), decreased among non-tumor cells, while the proportion of myeloid cells correspondingly increased. TIGIT, LAG3, and CD96 were the predominant checkpoint molecules expressed on exhausted T cells and CTLs; only TIGIT was significantly increased after relapse. CTLs expanded during remission, and then contracted during relapse with upregulated TIGIT expression. Tumor cells also acquired TIGIT expression after relapse, leading to the enhanced interaction of tumor cell TIGIT with monocyte CD155/PVR. In myeloid cells, post-relapse HLA-II expression was reduced relative to pretreatment and during remission. Myeloid-derived suppressor cells (MDSCs) were enriched after relapse with elevated expression of activation markers, including CLU (clusterin) and VCAN (versican). Extracellular chemokines (CCL4, CXCL9, CXCL13), soluble checkpoint inhibitors (sPD-L1, sTIM3, s4-1BB), and soluble receptors (sIL-2R, sTNFRII) were decreased during remission but elevated after relapse. CONCLUSIONS: Our data demonstrate that multiple tumor-intrinsic and -extrinsic factors are associated with T-cell suppression and BA relapse. Among these, TIGIT appears to be the central player given its elevated expression after BA relapse in not only CTLs but also MCL cells. The acquisition of TIGIT expression on tumor cells is MCL-specific and has not been reported in other CAR T-treated diseases. Together, our data suggest that co-targeting TIGIT may prevent CAR T relapses and thus promote long-term progression-free survival in MCL patients.

Follicular Lymphoma Microenvironment Characteristics Associated with Tumor Cell Mutations and MHC Class II Expression.

Follicular lymphoma (FL) is a B-cell malignancy with a complex tumor microenvironment that is rich in nonmalignant immune cells. We applied single-cell RNA sequencing to characterize the diverse tumor and immune cell populations of FL and identified major phenotypic subsets of FL T cells, including a cytotoxic CD4 T-cell population. We characterized four major FL subtypes with differential representation or relative depletion of distinct T-cell subsets. By integrating exome sequencing, we observed that somatic mutations are associated with, but not definitive for, reduced MHC expression on FL cells. In turn, expression of MHCII genes by FL cells was associated with significant differences in the proportions and targetable immunophenotypic characteristics of T cells. This provides a classification framework of the FL microenvironment in association with FL genotypes and MHC expression, and informs different potential immunotherapeutic strategies based upon tumor cell MHCII expression. SIGNIFICANCE: We have characterized the FL-infiltrating T cells, identified cytotoxic CD4 T cells as an important component that is associated with tumor cell-intrinsic characteristics, and identified sets of targetable immune checkpoints on T cells that differed from FLs with normal versus low MHC expression. See related commentary by Melnick, p. 374. This article is highlighted in the In This Issue feature, p. 369.

Androgen receptor blockade promotes response to BRAF/MEK-targeted therapy.

Treatment with therapy targeting BRAF and MEK (BRAF/MEK) has revolutionized care in melanoma and other cancers; however, therapeutic resistance is common and innovative treatment strategies are needed1,2. Here we studied a group of patients with melanoma who were treated with neoadjuvant BRAF/MEK-targeted therapy ( NCT02231775 , n = 51) and observed significantly higher rates of major pathological response (MPR; ≤10% viable tumour at resection) and improved recurrence-free survival (RFS) in female versus male patients (MPR, 66% versus 14%, P = 0.001; RFS, 64% versus 32% at 2 years, P = 0.021). The findings were validated in several additional cohorts2-4 of patients with unresectable metastatic melanoma who were treated with BRAF- and/or MEK-targeted therapy (n = 664 patients in total), demonstrating improved progression-free survival and overall survival in female versus male patients in several of these studies. Studies in preclinical models demonstrated significantly impaired anti-tumour activity in male versus female mice after BRAF/MEK-targeted therapy (P = 0.006), with significantly higher expression of the androgen receptor in tumours of male and female BRAF/MEK-treated mice versus the control (P = 0.0006 and P = 0.0025). Pharmacological inhibition of androgen receptor signalling improved responses to BRAF/MEK-targeted therapy in male and female mice (P = 0.018 and P = 0.003), whereas induction of androgen receptor signalling (through testosterone administration) was associated with a significantly impaired response to BRAF/MEK-targeted therapy in male and female patients (P = 0.021 and P < 0.0001). Together, these results have important implications for therapy.

MTAP deficiency creates an exploitable target for antifolate therapy in 9p21-loss cancers.

Methylthioadenosine phosphorylase, an essential enzyme for the adenine salvage pathway, is often deficient (MTAPdef) in tumors with 9p21 loss and hypothetically renders tumors susceptible to synthetic lethality by antifolates targeting de novo purine synthesis. Here we report our single arm phase II trial (NCT02693717) that assesses pemetrexed in MTAPdef urothelial carcinoma (UC) with the primary endpoint of overall response rate (ORR). Three of 7 enrolled MTAPdef patients show response to pemetrexed (ORR 43%). Furthermore, a historic cohort shows 4 of 4 MTAPdef patients respond to pemetrexed as compared to 1 of 10 MTAP-proficient patients. In vitro and in vivo preclinical data using UC cell lines demonstrate increased sensitivity to pemetrexed by inducing DNA damage, and distorting nucleotide pools. In addition, MTAP-knockdown increases sensitivity to pemetrexed. Furthermore, in a lung adenocarcinoma retrospective cohort (N = 72) from the published BATTLE2 clinical trial (NCT01248247), MTAPdef associates with an improved response rate to pemetrexed. Our data demonstrate a synthetic lethal interaction between MTAPdef and de novo purine inhibition, which represents a promising therapeutic strategy for larger prospective trials.

Inhibiting Type I Arginine Methyltransferase Activity Promotes T Cell-Mediated Antitumor Immune Responses.

Protein arginine methyltransferases (PRMT) are a widely expressed class of enzymes responsible for catalyzing arginine methylation on numerous protein substrates. Among them, type I PRMTs are responsible for generating asymmetric dimethylarginine. By controlling multiple basic cellular processes, such as DNA damage responses, transcriptional regulation, and mRNA splicing, type I PRMTs contribute to cancer initiation and progression. A type I PRMT inhibitor, GSK3368715, has been developed and has entered clinical trials for solid and hematologic malignancies. Although type I PRMTs have been reported to play roles in modulating immune cell function, the immunologic role of tumor-intrinsic pathways controlled by type I PRMTs remains uncharacterized. Here, our The Cancer Genome Atlas dataset analysis revealed that expression of type I PRMTs associated with poor clinical response and decreased immune infiltration in patients with melanoma. In cancer cell lines, inhibition of type I PRMTs induced an IFN gene signature, amplified responses to IFN and innate immune signaling, and decreased expression of the immunosuppressive cytokine VEGF. In immunocompetent mouse tumor models, including a model of T-cell exclusion that represents a common mechanism of anti-programmed cell death protein 1 (PD-1) resistance in humans, type I PRMT inhibition increased T-cell infiltration, produced durable responses dependent on CD8+ T cells, and enhanced efficacy of anti-PD-1 therapy. These data indicate that type I PRMT inhibition exhibits immunomodulatory properties and synergizes with immune checkpoint blockade (ICB) to induce durable antitumor responses in a T cell-dependent manner, suggesting that type I PRMT inhibition can potentiate an antitumor immunity in refractory settings.

Identification of MicroRNA-mRNA Networks in Melanoma and Their Association with PD-1 Checkpoint Blockade Outcomes.

Metastatic melanoma is a deadly malignancy with poor outcomes historically. Immuno-oncology (IO) agents, targeting immune checkpoint molecules such as cytotoxic T-lymphocyte associated protein-4 (CTLA-4) and programmed cell death-1 (PD-1), have revolutionized melanoma treatment and outcomes, achieving significant response rates and remarkable long-term survival. Despite these vast improvements, roughly half of melanoma patients do not achieve long-term clinical benefit from IO therapies and there is an urgent need to understand and mitigate mechanisms of resistance. MicroRNAs are key post-transcriptional regulators of gene expression that regulate many aspects of cancer biology, including immune evasion. We used network analysis to define two core microRNA-mRNA networks in melanoma tissues and cell lines corresponding to 'MITF-low' and 'Keratin' transcriptomic subsets of melanoma. We then evaluated expression of these core microRNAs in pre-PD-1-inhibitor-treated melanoma patients and observed that higher expression of miR-100-5p and miR-125b-5p were associated with significantly improved overall survival. These findings suggest that miR-100-5p and 125b-5p are potential markers of response to PD-1 inhibitors, and further evaluation of these microRNA-mRNA interactions may yield further insight into melanoma resistance to PD-1 inhibitors.

9p21 loss confers a cold tumor immune microenvironment and primary resistance to immune checkpoint therapy.

Immune checkpoint therapy (ICT) provides substantial clinical benefits to cancer patients, but a large proportion of cancers do not respond to ICT. To date, the genomic underpinnings of primary resistance to ICT remain elusive. Here, we performed immunogenomic analysis of data from TCGA and clinical trials of anti-PD-1/PD-L1 therapy, with a particular focus on homozygous deletion of 9p21.3 (9p21 loss), one of the most frequent genomic defects occurring in ~13% of all cancers. We demonstrate that 9p21 loss confers "cold" tumor-immune phenotypes, characterized by reduced abundance of tumor-infiltrating leukocytes (TILs), particularly, T/B/NK cells, altered spatial TILs patterns, diminished immune cell trafficking/activation, decreased rate of PD-L1 positivity, along with activation of immunosuppressive signaling. Notably, patients with 9p21 loss exhibited significantly lower response rates to ICT and worse outcomes, which were corroborated in eight ICT trials of >1,000 patients. Further, 9p21 loss synergizes with PD-L1/TMB for patient stratification. A "response score" was derived by incorporating 9p21 loss, PD-L1 expression and TMB levels in pre-treatment tumors, which outperforms PD-L1, TMB, and their combination in identifying patients with high likelihood of achieving sustained response from otherwise non-responders. Moreover, we describe potential druggable targets in 9p21-loss tumors, which could be exploited to design rational therapeutic interventions.

Resolving the Spatial and Cellular Architecture of Lung Adenocarcinoma by Multiregion Single-Cell Sequencing.

Little is known of the geospatial architecture of individual cell populations in lung adenocarcinoma (LUAD) evolution. Here, we perform single-cell RNA sequencing of 186,916 cells from five early-stage LUADs and 14 multiregion normal lung tissues of defined spatial proximities from the tumors. We show that cellular lineages, states, and transcriptomic features geospatially evolve across normal regions to LUADs. LUADs also exhibit pronounced intratumor cell heterogeneity within single sites and transcriptional lineage-plasticity programs. T regulatory cell phenotypes are increased in normal tissues with proximity to LUAD, in contrast to diminished signatures and fractions of cytotoxic CD8+ T cells, antigen-presenting macrophages, and inflammatory dendritic cells. We further find that the LUAD ligand-receptor interactome harbors increased expression of epithelial CD24, which mediates protumor phenotypes. These data provide a spatial atlas of LUAD evolution, and a resource for identification of targets for its treatment. SIGNIFICANCE: The geospatial ecosystem of the peripheral lung and early-stage LUAD is not known. Our multiregion single-cell sequencing analyses unravel cell populations, states, and phenotypes in the spatial and ecologic evolution of LUAD from the lung that comprise high-potential targets for early interception.This article is highlighted in the In This Issue feature, p. 2355.

Longitudinal single-cell profiling reveals molecular heterogeneity and tumor-immune evolution in refractory mantle cell lymphoma.

The mechanisms driving therapeutic resistance and poor outcomes of mantle cell lymphoma (MCL) are incompletely understood. We characterize the cellular and molecular heterogeneity within and across patients and delineate the dynamic evolution of tumor and immune cell compartments at single cell resolution in longitudinal specimens from ibrutinib-sensitive patients and non-responders. Temporal activation of multiple cancer hallmark pathways and acquisition of 17q are observed in a refractory MCL. Multi-platform validation is performed at genomic and cellular levels in PDX models and larger patient cohorts. We demonstrate that due to 17q gain, BIRC5/survivin expression is upregulated in resistant MCL tumor cells and targeting BIRC5 results in marked tumor inhibition in preclinical models. In addition, we discover notable differences in the tumor microenvironment including progressive dampening of CD8+ T cells and aberrant cell-to-cell communication networks in refractory MCLs. This study reveals diverse and dynamic tumor and immune programs underlying therapy resistance in MCL.

Single-Cell Expression Landscape of SARS-CoV-2 Receptor ACE2 and Host Proteases in Normal and Malignant Lung Tissues from Pulmonary Adenocarcinoma Patients.

The novel coronavirus SARS-CoV-2 is the causative agent of the COVID-19 pandemic. Severely symptomatic COVID-19 is associated with lung inflammation, pneumonia, and respiratory failure, thereby raising concerns of elevated risk of COVID-19-associated mortality among lung cancer patients. Angiotensin-converting enzyme 2 (ACE2) is the major receptor for SARS-CoV-2 entry into lung cells. The single-cell expression landscape of ACE2 and other SARS-CoV-2-related genes in pulmonary tissues of lung cancer patients remains unknown. We sought to delineate single-cell expression profiles of ACE2 and other SARS-CoV-2-related genes in pulmonary tissues of lung adenocarcinoma (LUAD) patients. We examined the expression levels and cellular distribution of ACE2 and SARS-CoV-2-priming proteases TMPRSS2 and TMPRSS4 in 5 LUADs and 14 matched normal tissues by single-cell RNA-sequencing (scRNA-seq) analysis. scRNA-seq of 186,916 cells revealed epithelial-specific expression of ACE2, TMPRSS2, and TMPRSS4. Analysis of 70,030 LUAD- and normal-derived epithelial cells showed that ACE2 levels were highest in normal alveolar type 2 (AT2) cells and that TMPRSS2 was expressed in 65% of normal AT2 cells. Conversely, the expression of TMPRSS4 was highest and most frequently detected (75%) in lung cells with malignant features. ACE2-positive cells co-expressed genes implicated in lung pathobiology, including COPD-associated HHIP, and the scavengers CD36 and DMBT1. Notably, the viral scavenger DMBT1 was significantly positively correlated with ACE2 expression in AT2 cells. We describe normal and tumor lung epithelial populations that express SARS-CoV-2 receptor and proteases, as well as major host defense genes, thus comprising potential treatment targets for COVID-19 particularly among lung cancer patients.