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.

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.

Abstract 2883: Gut microbiome dysbiosis promotes immune suppression and lung cancer development

Mounting evidence supports synergistic roles for the gut microbiome in cancer progression. Yet, the interplay between the gut microbiome and immune responses in cancer is still poorly understood. We recently showed that gut microbiome changes are closely associated with development of Kras-mutant lung adenocarcinoma (KM-LUAD) in a human-relevant, tobacco-associated mouse model (Gprc5a-/-; G). Knockout of the antimicrobial protein Lcn2 in these mice (Gprc5a-/-/Lcn2-/-; GL) further reduced microbial diversity while enhancing inflammation and tumor development. We thus hypothesized that microbial dysbiosis in the gut, such as that incurred by loss of Lcn2, may exacerbate LUAD development. Here, we investigated the effects of gut microbiome modulation on LUAD pathogenesis using fecal microbiota transfer (FMT) in both syngeneic and tobacco carcinogenesis models. Syngeneic G mice (transplant of G LUAD cells) that received FMT from GL donors (G < GL) exhibited significantly increased tumor growth relative to littermates with FMT from G mice (G < G). These effects were recapitulated in an independent syngeneic model (KrasG12D LKR13 cells in wild type mice). Tobacco carcinogen-exposed G < GL mice also exhibited increased lung tumor development compared with similarly exposed G < G littermates. 16S rDNA-Seq analysis of fecal pellets revealed significant differences in gut beta diversity between syngeneic G < G and G < GL mice. G < GL mice additionally displayed elevated relative abundance of tumor-promoting Alistipes, while Ruminoccocus and Akkermansia, taxa associated with favorable response to immunotherapy, were reduced. We next performed single-cell RNA-sequencing to comprehensively probe the tumor immune microenvironment (TIME) and the immune milieu near the gut of tumors and mesenteric lymph nodes (MLNs), respectively. The TIME in G < GL mice displayed an overall enhanced immunosuppressive phenotype evidenced by prominently increased fractions of T regulatory and Cd4+ Izumo1r+ exhausted T cells and, conversely, reduced levels of activated Isg15+ Cd8a+ T cells. MLNs from G < GL mice showed markedly increased fractions of memory B cells expressing the immunosuppressor Bank1 and reduced levels of follicular B cells and Cd8a+ Clec9a+ class 1 dendritic cells (cDC1). Flow cytometry further showed enhanced immunosuppression in G < GL relative to G < G mice, including increased fractions of myeloid-derived suppressor cells in the TIME of the former group. Our findings show that gut microbiome dysbiosis fosters lung cancer development by promoting immunosuppression, perhaps via a local and systemic gut microbiota-immune system crosstalk. Modulating the gut microbiome may be a promising strategy for interception or early treatment of lung cancer.

Citation Format: Zahraa Rahal, Fuduan Peng, Yuejiang Liu, Matthew C. Ross, Ansam Sinjab, Ke Liang, Jiping Feng, Chidera O. Chukwuocha, Manvi Sharma, Elizabeth Tang, Camille Abaya, Joseph Petrosino, Junya Fujimoto, Seyed Javad Moghaddam, Linghua Wang, Kristi L. Hoffman, Humam Kadara. Gut microbiome dysbiosis promotes immune suppression and lung cancer development [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 2883.

Abstract 113: An atlas of epithelial cell states and plasticity in lung adenocarcinoma

Understanding cellular processes underlying early lung adenocarcinoma (LUAD) development is needed to devise intervention strategies. While most if not all single-cell RNA sequencing (scRNA-seq) studies of lung cancer provided details on immune and stromal states, little insight is drawn to epithelial cells given their paucity (~4%) when performing unbiased scRNA-seq analysis without prior enrichment. Here, we performed in-depth scRNA-seq of enriched (by sorting for EPCAM+) epithelial cell subsets from 16 early-stage LUADs and 47 matching normal lung (NL) tissues. We also studied tissues from the same LUADs and adjacent NL by whole exome sequencing and a subset by high-resolution spatial protein and transcriptome analysis. We also performed scRNA-seq analysis of murine lungs from a human-relevant model of LUAD development following exposure to tobacco carcinogen, including strains with an alveolar type II (AT2) cell-specific lineage reporter. After extensive quality control, we retained 246,102 high quality human epithelial cells which comprised diverse normal alveolar and airway lineages as well as cancer cell populations. Diversity among cancer cells was strongly linked to LUAD oncogenic drivers. KRAS-mutant cancer cells were unique in their transcriptional features, strikingly reduced differentiation, low levels of copy number changes, and increased variability amongst the cells themselves. The local epithelial niche of LUADs, relative to that of NL, was enriched with intermediary cells in lung alveolar differentiation. A subset of these cells displayed elevated KRT8 expression (KRT8+ alveolar cells; KACs), increased plasticity and frequency of KRASG12D mutations, and its gene expression profiles were enriched in lung precancer and LUAD and signified poor survival. Notably, KACs harboring KRAS mutations were only found in the ecosystem of KRAS-mutant LUADs. Murine KACs were evident in lungs of tobacco carcinogen-exposed mice that develop KRAS-mutant LUADs but not in the saline-treated control group. While murine KACs emerged prior to tumor onset, they persisted for months after carcinogen cessation, and like their human counterparts, acquired driver Kras mutations, were poorly differentiated, and harbored KRAS-specific transcriptional programs. Spatial transcriptomics analysis showed that KAC and KRAS signatures were elevated in both murine and human tumors as well as in KACs that were in the local spatial vicinity of the LUADs. Organoids derived from lungs of tumor-bearing reporter mice were markedly enriched with KACs and were conspicuously sensitive to targeted inhibition of KRAS-G12D. This study provides new insights into the landscape of normal epithelial and malignant cells in LUAD, the role of alveolar intermediate subsets in development of the malignancy, particularly that driven by mutant KRAS, and, thus, potential targets for early interception.

Citation Format: Guangchun Han, Ansam Sinjab, Warapen Treekitkarnmongkol, Zahraa Rahal, Yuejiang Liu, Alejandra G. Serrano, Jiping Feng, Ke Liang, Khaja Khan, Wei Lu, Sharia Hernandez, Xuanye Cao, Enyu Dai, Yunhe Liu, Guangsheng Pei, Jian Hu, Lorena I. Gomez Bolanos, Edwin R. Parra, Tina Cascone, Boris Sepesi, Seyed Javad Moghaddam, Paul Scheet, Marcelo V. Negrao, John V. Heymach, Mingyao Li, Jichao Chen, Steven M. Dubinett, Junya Fujimoto, Luisa M. Solis, Ignacio I. Wistuba, Christopher S. Stevenson, Avrum Spira, Linghua Wang, Humam Kadara. An atlas of epithelial cell states and plasticity in lung adenocarcinoma [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 113.

Abstract 5104: Targeting IL-1B synergizes with PD-1 blockade for enhanced T and B cell immune responses and inhibition of early lung cancer development

Early immunotherapy with inhibitors of immune checkpoints such as PD-1 has revolutionized lung adenocarcinoma (LUAD) treatment. Still, many patients do not respond or relapse following PD-1 blockade. Tumor-promoting inflammation, such as that mediated by the pleiotropic cytokine interleukin 1 beta (IL-1B), fosters immunosuppression in the tumor immune microenvironment (TIME). Our group and others showed that IL-1B blockade inhibits lung cancer development. We thus hypothesized that addition of IL-1B blockade to anti-PD-1 treatment may enhance outcomes against LUAD. Using a human-relevant, tobacco-associated, mouse model of LUAD development, we compared the effects of combined PD-1 and IL-1B blockade relative to treatment with single-agents (anti-IL-1B or anti-PD-1) and control antibody on early lung tumor development and the TIME. Drugs (anti-PD-1 + anti-IL-1B, anti-PD-1, anti-IL-1B, control IgG) were administered at end of exposure to the tobacco-specific carcinogen NNK, to evaluate effects on formation of early lesions (preventive), or at 3 months post-NNK (prophylactic) to interrogate LUAD development (8 groups). Comprehensive interrogation of the lung ecosystem and the TIME was performed using deep single-cell RNA-sequencing (scRNA-seq) analysis in a subset of the mice (n = 3 to 4) from each of the 8 groups (n = 31 total; 143,897 cells after stringent quality control). Mice treated with combined PD-1 and IL-1B blockade displayed reduced development of lung tumors when compared to animals treated with anti-IL-1B, anti-PD-1, or control antibodies. Fractions of cytotoxic Cd8+ T cells were conspicuously higher and those of tumor cells and exhausted Cd8+ T cells evidently lower in lungs of mice treated with combined PD-1 and IL-1B blockade relative to monotherapy- or control antibody-treated animals. Igha+ plasma cells were strikingly highest in lungs of mice treated with combined PD-1 and IL-1B blockade and nearly absent in monotherapy- and control-treated groups. Lungs of mice treated with combined PD-1 and IL-1B blockade showed higher fractions of Cd80+/Cd86+ memory B cells and, consistently, T follicular helper T cells, while exhibiting reduced fractions of naïve B and Cd24a+/Tgfb1+ B cells suggestive of enhanced activation of B cell responses by the combinatorial treatment. These effects were, overall, present, or much more pronounced, in animals that were prophylactically treated. Flow cytometry analysis of lung tissues and immune profiling of bronchioalveolar lavage fluid overall confirmed augmented immune cell responses by combined PD-1 and IL-1B blockade. Our findings show that blocking IL-1B synergizes with anti-PD-1 in regression of early tumor cells and reversal of immunosuppression. Combined blockade of PD-1 and IL-1B may be a promising strategy for early treatment of lung cancer that warrants further clinical studies.

Citation Format: Warapen Treekitkarnmongkol, Guangchun Han, Zahraa Rahal, Jiping Feng, Ansam Sinjab, Tina Cascone, Christopher S. Stevenson, Cheryl Sweeney, Matt Edwards, Avrum Spira, Junya Fujimoto, Seyed Javad Moghaddam, Linghua Wang, Humam Kadara. Targeting IL-1B synergizes with PD-1 blockade for enhanced T and B cell immune responses and inhibition of early lung cancer development. [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 5104.

Neoadjuvant chemotherapy plus nivolumab with or without ipilimumab in operable non-small cell lung cancer: the phase 2 platform NEOSTAR trial

Neoadjuvant ipilimumab + nivolumab (Ipi+Nivo) and nivolumab + chemotherapy (Nivo+CT) induce greater pathologic response rates than CT alone in patients with operable non-small cell lung cancer (NSCLC). The impact of adding ipilimumab to neoadjuvant Nivo+CT is unknown. Here we report the results and correlates of two arms of the phase 2 platform NEOSTAR trial testing neoadjuvant Nivo+CT and Ipi+Nivo+CT with major pathologic response (MPR) as the primary endpoint. MPR rates were 32.1% (7/22, 80% confidence interval (CI) 18.7-43.1%) in the Nivo+CT arm and 50% (11/22, 80% CI 34.6-61.1%) in the Ipi+Nivo+CT arm; the primary endpoint was met in both arms. In patients without known tumor EGFR/ALK alterations, MPR rates were 41.2% (7/17) and 62.5% (10/16) in the Nivo+CT and Ipi+Nivo+CT groups, respectively. No new safety signals were observed in either arm. Single-cell sequencing and multi-platform immune profiling (exploratory endpoints) underscored immune cell populations and phenotypes, including effector memory CD8+ T, B and myeloid cells and markers of tertiary lymphoid structures, that were preferentially increased in the Ipi+Nivo+CT cohort. Baseline fecal microbiota in patients with MPR were enriched with beneficial taxa, such as Akkermansia, and displayed reduced abundance of pro-inflammatory and pathogenic microbes. Neoadjuvant Ipi+Nivo+CT enhances pathologic responses and warrants further study in operable NSCLC. (ClinicalTrials.gov registration: NCT03158129 .).

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.

Game of clones: Battles in the field of carcinogenesis

Recent advances in bulk sequencing approaches as well as genomic decoding at the single-cell level have revealed surprisingly high somatic mutational burdens in normal tissues, as well as increased our understanding of the landscape of "field cancerization", that is, molecular and immune alterations in mutagen-exposed normal-appearing tissues that recapitulated those present in tumors. Charting the somatic mutational landscapes in normal tissues can have strong implications on our understanding of how tumors arise from mutagenized epithelium. Making sense of those mutations to understand the progression along the pathologic continuum of normal epithelia, preneoplasias, up to malignant tissues will help pave way for identification of ideal targets that can guide new strategies for preventing or eliminating cancers at their earliest stages of development. In this review, we will provide a brief history of field cancerization and its implications on understanding early stages of cancer pathogenesis and deviation from the pathologically "normal" state. The review will provide an overview of how mutations accumulating in normal tissues can lead to a patchwork of mutated cell clones that compete while maintaining an overall state of functional homeostasis. The review also explores the role of clonal competition in directing the fate of normal tissues and summarizes multiple mechanisms elicited in this phenomenon and which have been linked to cancer development. Finally, we highlight the importance of understanding mutations in normal tissues, as well as clonal competition dynamics (in both the epithelium and the microenvironment) and their significance in exploring new approaches to combatting cancer.

Cell-by-Cell: Unlocking Lung Cancer Pathogenesis

For lung cancers, cellular trajectories and fates are strongly pruned by cell intrinsic and extrinsic factors. Over the past couple of decades, the combination of comprehensive molecular and genomic approaches, as well as the use of relevant pre-clinical models, enhanced micro-dissection techniques, profiling of rare preneoplastic lesions and surrounding tissues, as well as multi-region tumor sequencing, have all provided in-depth insights into the early biology and evolution of lung cancers. The advent of single-cell sequencing technologies has revolutionized our ability to interrogate these same models, tissues, and cohorts at an unprecedented resolution. Single-cell tracking of lung cancer pathogenesis is now transforming our understanding of the roles and consequences of epithelial-microenvironmental cues and crosstalk during disease evolution. By focusing on non-small lung cancers, specifically lung adenocarcinoma subtype, this review aims to summarize our knowledge base of tumor cells-of-origin and tumor-immune dynamics that have been primarily fueled by single-cell analysis of lung adenocarcinoma specimens at various stages of disease pathogenesis and of relevant animal models. The review will provide an overview of how recent reports are rewriting the mechanistic details of lineage plasticity and intra-tumor heterogeneity at a magnified scale thanks to single-cell studies of early- to late-stage lung adenocarcinomas. Future advances in single-cell technologies, coupled with analysis of minute amounts of rare clinical tissues and novel animal models, are anticipated to help transform our understanding of how diverse micro-events elicit macro-scale consequences, and thus to significantly advance how basic genomic and molecular knowledge of lung cancer evolution can be translated into successful targets for early detection and prevention of this lethal disease.

Chronic Exposure to Waterpipe Smoke Elicits Immunomodulatory and Carcinogenic Effects in the Lung

Effects of waterpipe smoking on lung pathobiology and carcinogenesis remain sparse despite the worldwide emergence of this tobacco vector. To address this gap, we investigated the effects of chronic waterpipe smoke (WPS) exposure on lung pathobiology, host immunity, and tumorigenesis using an experimental animal model that is prone to tobacco carcinogens and an exploratory observational analysis of human waterpipe smokers and nonsmokers. Mice exhibited elevated incidence of lung tumors following heavy WPS exposure (5 days/week for 20 weeks) compared to littermates with light WPS (once/week for 20 weeks) or control air. Lungs of mice exposed to heavy WPS showed augmented CD8+ and CD4+ T cell counts along with elevated protumor immune phenotypes including increased IL17A in T/B cells, PD-L1 on tumor and immune cells, and the proinflammatory cytokine IL1β in myeloid cells. RNA-sequencing (RNA-seq) analysis showed reduced antitumor immune gene signatures in animals exposed to heavy WPS relative to control air. We also performed RNA-seq analysis of airway epithelia from bronchial brushings of cancer-free waterpipe smokers and nonsmokers undergoing diagnostic bronchoscopy. Transcriptomes of normal airway cells in waterpipe smokers, relative to waterpipe nonsmokers, harbored gene programs that were associated with poor clinical outcomes in patients with lung adenocarcinoma, alluding to a WPS-associated molecular injury, like that established in response to cigarette smoking. Our findings support the notion that WPS exhibits carcinogenic effects and constitutes a possible risk factor for lung cancer as well as warrant future studies that can guide evidence-based policies for mitigating waterpipe smoking.

PREVENTION RELEVANCE

Potential carcinogenic effects of waterpipe smoking are very poorly understood despite its emergence as a socially acceptable form of smoking. Our work highlights carcinogenic effects of waterpipe smoking in the lung and, thus, accentuate the need for inclusion of individuals with exclusive waterpipe smoking in prevention and smoking cessation studies.

Abstract 1586: An emerging role for inflammation-associated alveolar intermediate cells in early phenotypic development of KRAS-mutant lung adenocarcinoma

Lung adenocarcinomas (LUAD) represent the most common lung cancer subtype and frequently harbor somatic mutations in the KRAS oncogene (KM-LUADs). While enhanced screening has improved early diagnosis of KM-LUAD, patient prognosis remains moderate to poor. Decoding the earliest events driving KM-LUADs can inform of ideal targets for its interception. Previous work showed that tobacco carcinogen (NNK) exposure leads to a pervasive field of injury comprised of molecular (e.g., KRAS mutations) and inflammatory changes that are shared between LUADs and their adjacent normal-appearing ecosystem. We and others have also shown that early immune and inflammatory alterations are implicated in the progression of normal lung (NL) epithelia and premalignant lesions (PMLs) to KM-LUAD. Yet, we still do not know the identities of specific epithelial subsets or how they promote a field of injury and inspire KM-LUAD pathogenesis. Here, we performed single-cell RNA-sequencing (scRNA-seq) of lungs from a human-relevant mouse model that develops PMLs and somatic KM-LUADs following NNK exposure. Analysis of 203,991 cells including 19,513 epithelial subsets after NNK cessation and at the onset of KM-LUADs revealed a unique population of alveolar cells that closely associated with tumor inception. These cells were highly evident in NNK- but not in control saline-exposed animals. Trajectory analysis showed that tumor clones developed through these transitionary cells, henceforth referred to as alveolar intermediate cells (AICs). Notably, AICs persisted for months after NNK cessation and acquired the same driver Kras mutations found in the resultant LUADs, thus supporting a role for AICs as KM-LUAD progenitors. Intriguingly, AICs harbored elevated expression of key components of p53 signaling (Trp53, Cdkn2a) and pro-inflammatory responses (IL-1β receptor Il1r1, NF-κB), and augmented cell-cell communication with Il1b+ macrophages which were enriched in LUAD-bearing lungs. Indeed, targeting IL-1β attenuated KM-LUAD development and increased anti-tumor immunity. In parallel, murine AIC expression profiles were significantly enriched in transcriptomes of human PMLs and LUADs. We thus probed our in-house and expanding scRNA-seq cohort of enriched (by sorting) epithelial subsets from human LUADs and NL, including 191,491 alveolar cells. AICs were not only evident in human lung tissues, but their fractions were also significantly increased in LUADs relative to NL. In conclusion, we identified a unique alveolar cell state that typified KM-LUAD progenitors, associated with inflammatory cues, and progressed along the pathologic continuum of damaged epithelium to KM-LUADs. Ongoing studies are evaluating whether AICs, in concert with tumor-initiating inflammation, trigger a field of injury that may underlie early phenotypic initiation and development of KM-LUAD.

Citation Format: Ansam Sinjab, Guangchun Han, Warapen Treekitkarnmongkol, Dapeng Hao, Enyu Dai, Luisa M. Solis, Seyed Javad Moghaddam, Junya Fujimoto, Jichao Chen, Matthew Edwards, Christopher S. Stevenson, Avrum E. Spira, Linghua Wang, Humam Kadara. An emerging role for inflammation-associated alveolar intermediate cells in early phenotypic development of KRAS-mutant lung adenocarcinoma [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 1586.

Abstract 4077: Immune suppressive phenotypes and pre-existing bystander T cell repertoire in early-stage lung adenocarcinoma revealed by integrative single-cell profiling

Lung adenocarcinomas (LUADs) constitute the most frequently diagnosed lung cancer subtype and are typified by poor clinical outcome and inferior response to therapy. Strategies to treat LUAD at its earliest stages heavily rely on improved understanding of its pathogenesis. Accumulating evidence shows that progression of cancers, including LUADs, is tightly linked to a dynamic interplay between cancer cells and the tumor immune microenvironment (TIME), which can block activation of tumor-infiltrating T cells (TILs), thereby hindering an effective anti-tumor immune response. Yet, how the TIME is shaped and the mechanisms mediating immune evasion in early-stage LUAD remain poorly understood. To fill this void, we dissected the phenotypic and clonotypic states as well as interaction networks of TILs by performing deep cellular and molecular profiling of 225,641 T cell transcriptomes and their T cell receptor (TCR) repertoire in 16 early-stage LUADs and matched normal lung (NL) samples using simultaneous single-cell RNA and TCR sequencing (scRNA/TCR-seq). Clustering analysis revealed a highly heterogeneous T cell landscape with distinct cell states identified including naïve, effector, regulatory, helper and exhausted cell subsets, as well as cellular phenotypes indicating stress and interferon response. The landscape of TILs in LUAD tissues was far more heterogeneous when compared to NLs, and the compositions and phenotypes of TILs differed significantly across patients and in association with clinicopathological variables. Relative to NLs, LUADs were specifically enriched with regulatory and stress response subsets, whereas fractions of cytotoxic T cell subsets (e.g., CD4+ CTLs, CD8+ tissue-resident memory/TRMs, and CD8+ natural killer T/NKT cells) were decreased. scTCR-seq analysis revealed significantly increased TCR clonotype diversity and reduced clonality in LUADs compared to NLs. We then tracked phenotype transitions of TILs by integrative analysis of TCR clonotypes and transcriptional states. Intriguingly, expanded TCR clonotypes were mostly TILs with effector phenotypes that were largely shared between LUADs and NLs, suggesting that their clonal expansion may not be tumor-specific. Of note, expanded clones were also enriched with immunosuppressive phenotypes (e.g., regulatory, exhausted CD8, and stress response T cells). In contrast, we found that contracted TCR clonotypes were characterized by memory/effector phenotypes. Together, our results reveal and characterize a diverse TIL landscape in early-stage LUAD encompassing a pre-existing and expanded T cell compartment that is likely tumor-unspecific, as well as switching in T cell phenotypes, and that together comprise viable targets for early immunotherapeutic interception of this malignancy.

Citation Format: Guangchun Han, Ansam Sinjab, Dapeng Hao, Lorena Gomez Bolanos, Enyu Dai, Luisa Maren Solis, Edwin Parra, Stephen Swisher, Tina Cascone, Boris Sepesi, Junya Fujimoto, Steven Dubinett, Ignacio Wistuba, Christopher Stevenson, Avrum Spira, Humam Kadara, Linghua Wang. Immune suppressive phenotypes and pre-existing bystander T cell repertoire in early-stage lung adenocarcinoma revealed by integrative single-cell profiling [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 4077.

Abstract 2126: Single-cell sequencing of early-stage lung adenocarcinomas reveals prominent intratumoral heterogeneity and epithelial plasticity programs

Decoding the complex molecular and cellular processes during lung adenocarcinoma (LUAD) development is needed to devise early intervention strategies. To comprehensively capture LUAD neoplastic heterogeneity and cellular plasticity, we performed single-cell RNA-sequencing (scRNA-seq) of 257,481 enriched epithelial cells (EPCAM+ sorting) from 16 early-stage LUADs, each with 3 matched normal lung (NL) samples at defined spatial proximities to the tumor (n=47). 29,076 LUAD-derived cells clustered by patient and harbored distinct gene expression features (e.g., oxidative stress response), signifying interpatient LUAD heterogeneity. We also identified, using whole exome sequencing (WES) of matching lung and germline control samples, recurrent oncogenic driver alterations (e.g., EGFR, TP53, KRAS). Transcriptomic features of malignant cells were shared between LUADs (e.g., loss of lineage-specific gene expression) or private such as those associated with driver mutation status (e.g., KRAS). Indeed, clusters of malignant cells were overall segregated based on driver mutations (e.g., KRAS, EGFR). Malignant cells from KRAS-mutant LUADs (KM-LUADs) had increased activation of NF-kB, estrogen and hypoxia signaling, comprising a unique gene module (GM) that correlated with a less differentiated state. We also found hallmark pathways (cholesterol metabolism, DNA replication, cell fate decision) specific to EGFR-mutant LUADs (EM-LUADs). Notably, cells from one EM-LUAD and its 3 multiregion NL tissues clustered closely and had activated pro-tumor lymphoid signatures (CD4 naïve, Treg). Mutation burden increased with tumor proximity and intriguingly, EGFR exon20 mutation was evident in the tumor (VAF = 0.29) and its most proximal NL (VAF = 0.05), signifying a mutational field effect. Copy number variations (CNVs) derived from WES of all samples were overall consistent with those inferred from scRNA-seq data. Relative to EM-LUADs, malignant cells from KM-LUADs displayed lower CNV burdens. Interpatient CNV heterogeneity was prominent even among LUADs harboring the same oncogenic drivers. Notably, intratumor heterogeneity (ITH) was high among epithelial cells within single regions from the same LUAD. Among LUADs, malignant cell clades with KRAS mutations and lower CNV scores displayed less differentiated states. To investigate biological pathways driving ITH, we derived 6 GMs with tumor-relevant functional features, including a transcription/translation regulation GM that consistently correlated with reduced differentiation. Our analysis of a large number of lung epithelial cells from LUAD patients reveals in-depth insights into LUAD taxonomy which can help identify epithelial heterotypes, unravel the continuum of early differentiation events and expand our understanding of early LUAD pathogenesis.

Citation Format: Guangchun Han, Ansam Sinjab, Warapen Treekitkarnmongkol, Dapeng Hao, Enyu Dai, Luisa M. Solis, Edwin R. Parra, Stephen Swisher, Tina Cascone, Boris Sepesi, Jichao Chen, Steven Dubinett, Junya Fujimoto, Ignacio I. Wistuba, Christopher S. Stevenson, Avrum E. Spira, Linghua Wang, Humam Kadara. Single-cell sequencing of early-stage lung adenocarcinomas reveals prominent intratumoral heterogeneity and epithelial plasticity programs [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 2126.

Abstract 6148: Immunogenomic landscape of tumor-infiltrating B and plasma cells in early-stage lung adenocarcinoma

Lung adenocarcinoma (LUAD) is the leading cause of cancer-related deaths in lifetime smokers. Relative to recent advances in immunotherapy of advanced-stage LUAD, there are very limited strategies for early immune-based treatment or interception of the malignancy in its primitive stages. This is largely due to a poor understanding of the roles and functional phenotypes of distinct immune cell subsets and how they evolve early in LUAD pathogenesis. For instance, while T cells have been a central focus of cancer immunopathology and immunotherapy, the roles of tumor-infiltrating B and plasma cells (TIBs) in the activity of the adaptive immune response along the pathogenic course of solid tumors such as LUAD are extremely poorly understood. To fill these voids, we conducted pan-cancer single-cell RNA sequencing (scRNA-seq) analysis of TIBs using public and in-house cohorts of >15 cancers and matched normal samples. We found that fractions of TIBs, including plasma cells (PCs), were evidently high in the tumor microenvironment (TME) of LUADs, particularly in smokers. We then performed multi-region paired scRNA-seq and single-cell B cell receptor sequencing (scBCR-seq) of tumors and three matched normal lung (NL) tissue samples with varying spatial proximity from each of the tumors of 16 early-stage LUAD patients. Fractions of TIBs including PCs and memory B cells were immensely increased in the TME of early-stage LUADs compared to uninvolved NL, and conversely, the abundance of naïve B cells was markedly decreased. TIB fractions were progressively increased along the pathologic continuum of NL, premalignant lesions (PMLs), up to invasive LUAD. Consistently, expression of the B cell chemotactic CXCL13 - CXCR5 axis in T cells and TIBs, respectively, was increased in both PMLs and LUADs relative to NL. Simultaneous scBCR-seq revealed markedly reduced clonality of BCR repertoires in LUAD compared to NL. Multi-region NL tissues showed progressively increasing BCR clonotype diversity and immunoglobin somatic hypermutation (SHM) with increasing proximity to tumors. BCR clonality was strikingly lower in smoker LUADs relative to non-smoker tumors as well as progressively attenuated with increasing pathologic stage. TIBs in the LUAD TME were mostly composed of terminally differentiated IgA+ or IgG+ PCs and memory B cells with an immunosuppressive phenotype. To understand how TIBs shape the LUAD TME, we also profiled interactions of TIBs with other TME cell components and identified TIB subsets showing strong co-occurrence patterns with immunosuppressive T cell subsets. By comprehensively defining transcriptional profiles, SHM and antibody repertories, as well as cellular interactions of TIBs at single-cell resolution, our results map out the geospatial landscape of TIBs in early-stage LUAD and provide a valuable resource to leverage targets for innovative immunomodulatory strategies.

Citation Format: Dapeng Hao, Guangchun Han, Ansam Sinjab, Lorena Gomez Bolanos, Rossana Lazcano Segura, Enyu Dai, Luisa Maren Solis Soto, Edwin Parra, Jennifer Wargo, Stephen Swisher, Tina Cascone, Boris Sepesi, Junya Fujimoto, Steven Dubinett, Ignacio Wistuba, Christopher Stevenson, Avrum Spira, Humam Kadara, Linghua Wang. Immunogenomic landscape of tumor-infiltrating B and plasma cells in early-stage lung adenocarcinoma [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 6148.

Abstract 223: Chronic exposure to waterpipe smoke elicits immunomodulatory and carcinogenic effects in the lung

Effects of waterpipe smoking on lung pathobiology and carcinogenesis remain sparse despite the worldwide emergence of this tobacco vector as a socially acceptable form of smoking, particularly among the youth. To address this gap, we investigated the effects of chronic waterpipe smoke (WPS) exposure on lung pathobiology, host immunity, and tumorigenesis using both an experimental animal model and an exploratory observational analysis of human waterpipe smokers and non-smokers. Mice exposed to increasing doses of WPS (once or five times per week for 20 weeks), through an exclusively devised exposure system for this study, were more prone to develop lung tumors compared to control-air exposed littermates. This effect was accompanied by various pro-tumor immune phenotypes, including increased IL-17A+ levels in T/B cells as well as elevated expression of the immune checkpoint PD-L1 and the pro-inflammatory cytokine IL-1β in myeloid cells. While flow cytometry analysis revealed increased CD4 and CD8 T cell infiltration in WPS-exposed mice, gene signatures of cytotoxic and expanded immune response were, conversely, decreased in animals exposed to heavy WPS relative to control air. Interestingly, mice heavily exposed to WPS exhibited increased expression of Tmprss4, Cd55, and Ace2, cell receptors and mediators of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) entry and, thus, COVID-19 pathogenesis. We also perform RNA-sequencing analysis of bronchial airway epithelial brushings of cancer-free waterpipe smokers and non-smokers undergoing diagnostic bronchoscopy. Transcriptomes of normal airway cells in waterpipe smokers, relative to waterpipe non-smokers, harbored gene programs that were associated with poor clinical outcomes in lung adenocarcinoma (LUAD) patients, alluding to a WPS-associated molecular injury, like that established in response to cigarette smoking. Overall, our analyses demonstrate immunomodulatory and carcinogenic effects of WPS on the murine and human lung. Our study also shows that WPS exposure leads to a field of injury that may be associated with increased risk for lung cancer. This study is unique for interrogating carcinogenic effects of chronic exposure to WPS with the longest, to our knowledge, follow-up time in vivo. Our findings accentuate the need for additional studies that can guide evidence-based policies to counteract shortfalls in public health control of waterpipe smoking.

Citation Format: Zahraa Rahal, Maya Hassane, Nareg Karaoghlanian, Jiexin Zhang, Ansam Sinjab, J. Jack Lee, Maria Gabriela Raso, Junya Fujimoto, Hassan Chami, Alan Shihadeh, Humam Kadara. Chronic exposure to waterpipe smoke elicits immunomodulatory and carcinogenic effects in the lung [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 223.

Abstract PR006: An emerging role for inflammation-associated alveolar intermediate cells in early phenotypic development of KRAS-mutant lung adenocarcinoma

Lung adenocarcinomas (LUADs) represent the most common lung cancer subtype and frequently harbor somatic mutations in the KRAS oncogene (KM-LUADs). While enhanced screening has improved early diagnosis of KM-LUAD, patient prognosis remains moderate to poor. Decoding the earliest events driving KM-LUADs can inform of ideal targets for its interception. Previous work showed that tobacco carcinogen (NNK) exposure leads to a pervasive field of injury comprised of molecular (e.g., KRAS mutations) and inflammatory changes that are shared between LUADs and their adjacent normal-appearing ecosystem. We and others have also shown that early immune and inflammatory alterations are implicated in the progression of normal lung (NL) epithelia and premalignant lesions (PMLs) to KM-LUAD. Yet, we still do not know the identities of specific epithelial subsets or how they promote a field of injury and inspire KM-LUAD pathogenesis. Here, we performed single-cell RNA-sequencing (scRNA-seq) of lungs from a human-relevant mouse model that develops PMLs and somatic KM-LUADs following NNK exposure. Analysis of 203,991 cells including 19,513 epithelial subsets after NNK cessation and at the onset of KM-LUADs revealed a unique population of alveolar cells that closely associated with tumor inception. These cells were highly evident in NNK- but not in control saline-exposed animals. Trajectory analysis showed that tumor clones developed through these transitionary cells, henceforth referred to as alveolar intermediate cells (AICs). Notably, AICs persisted for months after NNK cessation and acquired the same driver Kras mutations found in the resultant LUADs, thus supporting a role for AICs as KM-LUAD progenitors. Intriguingly, AICs harbored elevated expression of key components of p53 signaling (Trp53, Cdkn2a) and pro-inflammatory responses (IL-1β receptor Il1r1, NF-κB), and augmented cell-cell communication with Il1b+ macrophages which were enriched in LUAD-bearing lungs. Indeed, targeting IL-1β attenuated KM-LUAD development and increased anti-tumor immunity. In parallel, murine AIC expression profiles were significantly enriched in transcriptomes of human PMLs and LUADs. We thus probed our in-house and expanding scRNA-seq cohort of enriched (by sorting) epithelial subsets from human LUADs and NL, including 191,491 alveolar cells. AICs were not only evident in human lung tissues, but their fractions were also significantly increased in LUADs relative to NL. In conclusion, we identified a unique alveolar cell state that typified KM-LUAD progenitors, associated with inflammatory cues, and progressed along the pathologic continuum of damaged epithelium to KM-LUADs. Ongoing studies are evaluating whether AICs, in concert with tumor-initiating inflammation, trigger a field of injury that may underlie early phenotypic initiation and development of KM-LUAD.

Citation Format: Ansam Sinjab, Guangchun Han, Warapen Treekitkarnmongkol, Dapeng Hao, Enyu Dai, Luisa M Solis, Seyed Javad Moghaddam, Junya Fujimoto, Jichao Chen, Matthew Edwards, Christopher S. Stevenson, Avrum E. Spira, Linghua Wang, Humam Kadara. An emerging role for inflammation-associated alveolar intermediate cells in early phenotypic development of KRAS-mutant lung adenocarcinoma [abstract]. In: Proceedings of the AACR Special Conference on the Evolutionary Dynamics in Carcinogenesis and Response to Therapy; 2022 Mar 14-17. Philadelphia (PA): AACR; Cancer Res 2022;82(10 Suppl):Abstract nr PR006.

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.

Early Diagnosis and Screening for Lung Cancer

Cancer interception refers to actively blocking the cancer development process by preventing progression of premalignancy to invasive disease. The rate-limiting steps for effective lung cancer interception are the incomplete understanding of the earliest molecular events associated with lung carcinogenesis, the lack of preclinical models of pulmonary premalignancy, and the challenge of developing highly sensitive and specific methods for early detection. Recent advances in cancer interception are facilitated by developments in next-generation sequencing, computational methodologies, as well as the renewed emphasis in precision medicine and immuno-oncology. This review summarizes the current state of knowledge in the areas of molecular abnormalities in lung cancer continuum, preclinical human models of lung cancer pathogenesis, and the advances in early lung cancer diagnostics.

Abstract 130: Resolving the spatial and cellular architecture of lung adenocarcinoma by multi-region single-cell sequencing

Lung adenocarcinoma (LUAD) is the most commonly diagnosed histological subtype of lung cancer. While earlier work has underscored genomic and immune alterations in LUAD, the roles of individual cell populations in early-stage human LUAD evolution in space remain unknown. Here, we provide a detailed cellular atlas of early-stage LUAD and its spatial ecosystem along the peripheral lung. We performed single-cell RNA sequencing of 186,916 cells including enriched epithelial fractions from five early-stage LUADs with fourteen multi-region normal lung tissues of defined spatial proximities from the primary LUADs. We show that major epithelial and immune cellular lineages, states, and transcriptomic features geospatially and progressively evolve across normal regions and with increasing LUAD proximity. Analysis of 70,030 lung epithelial cells unraveled diverse lineage trajectories, transcriptional lineage plasticity programs underlying KRAS-mutant cells, and intratumoral heterogeneity within single sites. T regulatory cell programs including multiple immune checkpoints increased in tissues with closer proximity to LUADs, in sharp contrast to signatures of CD8+ cytotoxic T cells, antigen presentation by macrophages, and inflammatory dendritic cells. We found that some spatial signatures (e.g. a B cell signature score) were increased along the pathologic spectrum of normal lung, preneoplastic lesions, and matched invasive LUADs. LUAD cell-cell communication networks were enriched with ligand-receptor interactions involving CD24, LGALS9 and TIM3 immune checkpoints, including crosstalk between CD24 antigen in LUAD epithelial cells and SIGLEC10 in myeloid subsets. CD24 was markedly increased in preneoplasias relative to normal lung and further in LUAD, and its expression was highly positively correlated with immunosuppressive phenotypes. These data provide an atlas of cellular states and phenotypes underlying early-stage LUAD evolution in space, and a scalable resource for identification of targets for early treatment.

Citation Format: Ansam Sinjab, Guangchun Han, Warapen Treekitkarnmongkol, Kieko Hara, Patrick Brennan, Minghao Dang, Dapeng Hao, Ruiping Wang, Enyu Dai, Hitoshi Dejima, Jiexin Zhang, Elena Bogatenkova, Beatriz Sanchez-Espiridion, Kyle Chang, Danielle R. Little, Samer Bazzi, Linh Tran, Kostyantyn Krysan, Carmen Behrens, Dzifa Duose, Edwin R. Parra, Maria Gabriela Raso, Luisa M. Solis, Junya Fukuoka, Jianjun Zhang, Boris Sepesi, Tina Cascone, Lauren A. Byers, Don L. Gibbons, Jichao Chen, Seyed Javad Moghaddam, Edwin J. Ostrin, Daniel G. Rosen, John V. Heymach, Paul Scheet, Steven Dubinett, Ignacio I. Wistuba, Junya Fujimoto, Christopher S. Stevenson, Avrum E. Spira, Linghua Wang, Humam Kadara. Resolving the spatial and cellular architecture of lung adenocarcinoma by multi-region single-cell sequencing [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 130.

Abstract 702: Single-cell expression landscape of SARS-CoV-2 receptor ACE2 and host proteases in human lung adenocarcinoma

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-asociated mortality among lung cancer patients. Angiotensin converting enzyme 2 (ACE2) is the major receptor for SARS-CoV-2 entry into lung cells. Yet, the single-cell expression landscape of ACE2 and other SARS-CoV-2-related genes in pulmonary tissues of lung cancer patients remains unknown. To fill these voids, we leveraged our ongoing efforts in single-cell transcriptomic analysis of 186,916 cells including a large number of epithelial cells (n=70,030) derived from 5 LUADs and 14 matching uninvolved normal lung tissues, to delineate expression levels and cellular distribution of ACE2 and SARS-CoV-2 priming proteases TMPRSS2 and TMPRSS4. Single-cell RNA sequencing 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, expression of TMPRSS4 was highest and most frequently detected (75%) in malignant lung cells. 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. In conclusion, we describe normal and tumor lung epithelial populations that express SARS-CoV-2 receptor and proteases, as well as major host defense genes, and that thus comprise potential treatment targets for COVID-19 particularly among lung cancer patients.

Citation Format: Guangchun Han, Ansam Sinjab, Kieko Hara, Warapen Treekitkarnmongkol, Patrick Brennan, Kyle Chang, Elena Bokatenkova, Beatriz Sanchez-Espiridion, Carmen Behrens, Luisa M. Solis, Boning Gao, Luc Girard, Jianjun Zhang, Boris Sepesi, Tina Cascone, Lauren A. Byers, Don L. Gibbons, Jichao Chen, Seyed Javad Moghaddam, Edwin J. Ostrin, Paul Scheet, Junya Fujimoto, Jerry Shay, John V. Heymach, John D. Minna, Steven Dubinett, Ignacio I. Wistuba, Christopher S. Stevenson, Avrum E. Spira, Linghua Wang, Humam Kadara. Single-cell expression landscape of SARS-CoV-2 receptor ACE2 and host proteases in human lung adenocarcinoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 702.

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

Simple Summary

The coronavirus disease 2019 (COVID-19) pandemic continues to spread rapidly on a global scale. When presenting with severe respiratory complications, COVID-19 results in markedly high death rates, particularly among patients with comorbidities such as cancer. Motivated by the ongoing global health crisis, we leveraged a growing in-house cohort of pulmonary tissues from lung cancer patients to analyze, at high resolution, the expression of host proteins implicated in the entryway of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) into lung epithelial cells. Our results identify key pathways in lung pathobiology and inflammation that offer the potential to identify novel markers and therapeutic targets that can be repurposed for clinical management of COVID-19, particularly among lung cancer patients, a population that represents over half a million individuals in the United States alone.

Abstract

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.

Augmented Lipocalin-2 Is Associated with Chronic Obstructive Pulmonary Disease and Counteracts Lung Adenocarcinoma Development

Rationale: Early pathogenesis of lung adenocarcinoma (LUAD) remains largely unknown. We found that, relative to wild-type littermates, the innate immunomodulator Lcn2 (lipocalin-2) was increased in normal airways from mice with knockout of the airway lineage gene Gprc5a (Gprc5a^-/-) and that are prone to developing inflammation and LUAD. Yet, the role of LCN2 in lung inflammation and LUAD is poorly understood.Objectives: Delineate the role of Lcn2 induction in LUAD pathogenesis.Methods: Normal airway brushings, uninvolved lung tissues, and tumors from Gprc5a^-/- mice before and after tobacco carcinogen exposure were analyzed by RNA sequencing. LCN2 mRNA was analyzed in public and in-house data sets of LUAD, lung squamous cancer (LUSC), chronic obstructive pulmonary disease (COPD), and LUAD/LUSC with COPD. LCN2 protein was immunohistochemically analyzed in a tissue microarray of 510 tumors. Temporal lung tumor development, gene expression programs, and host immune responses were compared between Gprc5a^-/- and Gprc5a^-/-/Lcn2^-/- littermates.Measurements and Main Results: Lcn2 was progressively elevated during LUAD development and positively correlated with proinflammatory cytokines and inflammation gene sets. LCN2 was distinctively elevated in human LUADs, but not in LUSCs, relative to normal lungs and was associated with COPD among smokers and patients with LUAD. Relative to Gprc5a^-/- mice, Gprc5a^-/-/Lcn2^-/- littermates exhibited significantly increased lung tumor development concomitant with reduced T-cell abundance (CD4⁺) and richness, attenuated antitumor immune gene programs, and increased immune cell expression of protumor inflammatory cytokines.Conclusions: Augmented LCN2 expression is a molecular feature of COPD-associated LUAD and counteracts LUAD development in vivo by maintaining antitumor immunity.

Field Carcinogenesis in Cancer Evolution: What the Cell Is Going On?

Field carcinogenesis describes the prevalence of tumor-related alterations in normal appearing tissues. Here, we summarize recent efforts in profiling molecular field dynamics for resolving early events in cancer evolution. We also highlight gaps in our knowledge of the molecular and cellular heterogeneity of field carcinogenesis and propose directions to tackle these voids using single-cell-based approaches and unique tissue sampling models. By interrogating both the mutagenized epithelium and its microenvironment, we surmise that single-cell-guided studies will help chart the spatiotemporal molecular and cellular "atlas" of field carcinogenesis, will further delineate preneoplastic initiation and progression, and will help identify cancer prevention and early intervention targets.

Abstract 1518: A single-cell transcriptomic atlas of lung adenocarcinoma and adjacent normal-appearing tissue

Lung adenocarcinoma (LUAD) is the most frequently diagnosed histological subtype of lung cancer and accounts for most smoking-related cancer deaths, warranting strategies for early intervention. Earlier work revealed genome-wide aberrations in LUADs and the adjacent premalignant field, known as “field carcinogenesis”, that are pertinent to LUAD pathogenesis. Yet, we still poorly understand the cellular and molecular architecture of LUAD and its nearby “field”. To fill this void, we performed, using the 10X Genomics system and NovaSeq 6000 platform, single-cell RNA sequencing (scRNA-seq) analysis of 4 early-stage resected LUADs as well as 11 matched normal lung tissues with differing spatial proximity from the tumors. We first analyzed 15,739 cells from a LUAD and matched tumor-adjacent and -distant normal tissues from patient 1. The resulting fraction of EPCAM+ cells from this analysis was approximately 4%. For deeper single-cell resolution of LUAD, we thus performed scRNA-seq of separately sorted/enriched epithelial cells (EPCAM+; n = 50,883) and non-epithelial cells (EPCAM-; n= 91,093) from patients 2, 3 and 4, each with an early-stage LUAD and three matched spatially-distributed normal-appearing tissues. Overall, we achieved on average approximately 150,000 reads and 1,956 genes detected per cell. Divergent populations of malignant and non-malignant cells and multiple epithelial and immune subsets clustered in an overall spatially modulated pattern according to proximity to the tumor. Hierarchical clustering revealed multiple distinct populations of airway lineage cells, including alveolar type 1 (AT1), alveolar type 2 (AT2), AT1/AT2 bipotential cells, club, goblet, basal and ciliated cells, with cells from the resected LUAD present in some but not all clusters. By analysis of epithelial-enriched fractions, we were also able to identify and interrogate rare cell subpopulations including CFTR-expressing ionocytes. By inference of copy number variation along with analysis of key oncogene expression patterns, we also identified epithelial cells that represent potential tumor cells-of-origin. Spatial reprogramming of field carcinogenesis encompassed immune cell populations including tumor-exclusive FOXP3+ regulatory T cells with marked over-expression of the major immune checkpoints CTLA4 and TIGIT. In contrast, for at least some of the cases, tumors were nearly devoid of natural killer cells, and abundance of CD8+ T cells dampened with increasing proximity to the tumor. Our single-cell surveys offer insights into novel cues in LUAD pathogenesis. Efforts are underway to interrogate the single-cell epithelial and immune landscape of additional LUADs and matched nearby normal-appearing lung, including those from non-smokers, to better understand the evolution of the disease and, thus, identify a low-hanging fruit of targets for early management of this fatal malignancy.

Citation Format: Ansam Sinjab, Guangchun Han, Warapen Treekitkarnmongkol, Patrick Brennan, Kieko Hara, Kyle Chang, Elena Bogatenkova, Beatriz Sanchez-Espiridion, Carmen Behrens, Jianjun Zhang, Boris Sepesi, Tina Cascone, Don L. Gibbons, Jichao Chen, George Eapen, Edwin J. Ostrin, Junya Fujimoto, Avrum E. Spira, Paul A. Scheet, Ignacio I. Wistuba, Linghua Wang, Humam Kadara. A single-cell transcriptomic atlas of lung adenocarcinoma and adjacent normal-appearing tissue [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 1518.

Anti-Tumor Effects of Biomimetic Sulfated Glycosaminoglycans on Lung Adenocarcinoma Cells in 2D and 3D In Vitro Models

Lung cancer development relies on cell proliferation and migration, which in turn requires interaction with extracellular matrix (ECM) components such as glycosaminoglycans (GAGs). The mechanisms through which GAGs regulate cancer cell functions are not fully understood but they are, in part, mediated by controlled interactions with cytokines and growth factors (GFs). In order to mechanistically understand the effect of the degree of sulfation (DS) of GAGs on lung adenocarcinoma (LUAD) cells, we synthesized sulfated alginate (AlgSulf) as sulfated GAG mimics with DS = 0.0, 0.8, 2.0, and 2.7. Human (H1792) and mouse (MDA-F471) LUAD cell lines were treated with AlgSulf of various DSs at two concentrations 10 and 100 µg/mL and their anti-tumor properties were assessed using 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT), trypan blue exclusion, and wound healing assays for 2D models and sphere formation assay for the 3D model. The proliferation and number of live MDA-F471 cells at the concentration of 100 µg/mL decreased significantly with the increase in the DS of biomimetic GAGs. In addition, the increase in the DS of biomimetic GAGs decreased cell migration (p < 0.001 for DS = 2.0 and 2.7 compared to control) and decreased the diameter and number of spheres formed (p < 0.001). The increased DS of biomimetic GAGs attenuated the expression of cancer stem cell (CSC)/progenitor markers in the 3D cultures. In conclusion, GAG-mimetic AlgSulf with increased DS exhibit enhanced anti-proliferative and migratory properties while also reducing growth of KRAS-mutant LUAD spheres in vitro. We suggest that these anti-tumor effects by GAG-mimetic AlgSulf are possibly due to differential binding to GFs and consequential decreased cell stemness. AlgSulf may be suitable for applications in cancer therapy after further in vivo validation.

Insights Into Lung Cancer Immune-Based Biology, Prevention, and Treatment

Lung cancer is the number one cause of cancer-related deaths. The malignancy is characterized by dismal prognosis and poor clinical outcome mostly due to advanced-stage at diagnosis, thereby inflicting a heavy burden on public health worldwide. Recent breakthroughs in immunotherapy have greatly benefited a subset of lung cancer patients, and more importantly, they are undauntedly bringing forth a paradigm shift in the drugs approved for cancer treatment, by introducing "tumor-type agnostic therapies". Yet, and to fulfill immunotherapy's potential of personalized cancer treatment, demarcating the immune and genomic landscape of cancers at their earliest possible stages will be crucial to identify ideal targets for early treatment and to predict how a particular patient will fare with immunotherapy. Recent genomic surveys of premalignant lung cancer have shed light on early alterations in the evolution of lung cancer. More recently, the advent of immunogenomic technologies has provided prodigious opportunities to study the multidimensional landscape of lung tumors as well as their microenvironment at the molecular, genomic, and cellular resolution. In this review, we will summarize the current state of immune-based therapies for cancer, with a focus on lung malignancy, and highlight learning outcomes from clinical and preclinical studies investigating the naïve immune biology of lung cancer. The review also collates immunogenomic-based evidence from seminal reports which collectively warrant future investigations of premalignancy, the tumor-adjacent normal-appearing lung tissue, pulmonary inflammatory conditions such as chronic obstructive pulmonary disease, as well as systemic microbiome imbalance. Such future directions enable novel insights into the evolution of lung cancers and, thus, can provide a low-hanging fruit of targets for early immune-based treatment of this fatal malignancy.

Genome-Wide and Phenotypic Evaluation of Stem Cell Progenitors Derived From Gprc5a-Deficient Murine Lung Adenocarcinoma With Somatic Kras Mutations

Lung adenocarcinomas (LUADs) with somatic mutations in the KRAS oncogene comprise the most common molecular subtype of lung cancer in smokers and present with overall dismal prognosis and resistance to most therapies. Our group recently demonstrated that tobacco carcinogen-exposed mice with knockout of the airway lineage G-protein coupled receptor, Gprc5a, develop LUADs with somatic mutations in Kras. Earlier work has suggested that cancer stem cells (CSCs) play crucial roles in clonal evolution of tumors and in therapy resistance. To date, our understanding of CSCs in LUADs with somatic Kras mutations remains lagging. Here we derived CSCs (as spheres in 3D cultures) with self-renewal properties from a murine Kras-mutant LUAD cell line we previously established from a tobacco carcinogen-exposed Gprc5a^−/− mouse. Using syngeneic Gprc5a^−/− models, we found that these CSCs, compared to their parental isoforms, exhibited increased tumorigenic potential in vivo, particularly in female animals. Using whole-transcriptome sequencing coupled with pathways analysis and confirmatory PCR, we identified gene features (n = 2,600) differentially expressed in the CSCs compared to parental cells and that were enriched with functional modules associated with an augmented malignant phenotype including stemness, tumor-promoting inflammation and anti-oxidant responses. Further, based on in silico predicted activation of GSK3β in CSCs, we found that tideglusib, an irreversible inhibitor of the kinase, exhibited marked anti-growth effects in the cultured CSCs. Our study underscores molecular cues in the pathogenesis of Kras-mutant LUAD and presents new models to study the evolution, and thus high-potential targets, of this aggressive malignancy.

Genome-Wide Gene Expression Changes in the Normal-Appearing Airway during the Evolution of Smoking-Associated Lung Adenocarcinoma

Smoking perpetuates in cytologically normal airways a molecular "field of injury" that is pertinent to lung cancer and early detection. The evolution of airway field changes prior to lung oncogenesis is poorly understood largely due to the long latency of lung cancer in smokers. Here, we studied airway expression changes prior to lung cancer onset in mice with knockout of the Gprc5a gene (Gprc5a^-/-) and tobacco carcinogen (NNK) exposure and that develop the most common type of lung cancer, lung adenocarcinoma, within 6 months following exposure. Airway epithelial brushings were collected from Gprc5a^-/- mice before exposure and at multiple times post-NNK until time of lung adenocarcinoma development and then analyzed by RNA sequencing. Temporal airway profiles were identified by linear models and analyzed by comparative genomics in normal airways of human smokers with and without lung cancer. We identified significantly altered profiles (n = 926) in the NNK-exposed mouse normal airways relative to baseline epithelia, a subset of which were concordantly modulated with smoking status in the human airway. Among airway profiles that were significantly modulated following NNK, we found that expression changes (n = 22) occurring as early as 2 months following exposure were significantly associated with lung cancer status when examined in airways of human smokers. Furthermore, a subset of a recently reported human bronchial gene classifier (Percepta; n = 56) was enriched in the temporal mouse airway profiles. We underscore evolutionarily conserved profiles in the normal-appearing airway that develop prior to lung oncogenesis and that comprise viable markers for early lung cancer detection in suspect smokers. Cancer Prev Res; 11(4); 237-48. ©2018 AACR.

Smoking and Lung Cancer: A Geo-Regional Perspective

Lung cancer is the leading cause of cancer-related deaths worldwide. Non-small cell lung cancer (NSCLC) represents the most frequently diagnosed subtype of this morbid malignancy. NSCLC is causally linked to tobacco consumption with more than 500 million smokers worldwide at high risk for this fatal malignancy. We are currently lagging in our knowledge of the early molecular (e.g., genomic) effects of smoking in NSCLC pathogenesis that would constitute ideal markers for early detection. This limitation is further amplified when considering the variable etiologic factors in NSCLC pathogenesis among different regions around the globe. In this review, we present our current knowledge of genomic alterations arising during early stages of smoking-induced lung cancer initiation and progression, including discussing the premalignant airway field of injury induced by smoking. The review also underscores the wider spectra and higher age-adjusted rates of tobacco (e.g., water-pipe smoke) consumption, along with elevated environmental carcinogenic exposures and relatively poorer socioeconomic status, in low-middle income countries (LMICs), with Lebanon as an exemplar. This “cocktail” of carcinogenic exposures warrants the pressing need to understand the complex etiology of lung malignancies developing in LMICs such as Lebanon.

Meeting report: the road to science-based policy - ESOF through the eyes of young scientists

Communication and common understanding between politicians, scientists, and the society can lead to evidence-based science policy, a core principle that guides high caliber research and open innovation for a sustainable future.

Parthenolide: from plant shoots to cancer roots

Parthenolide (PTL), a sesquiterpene lactone (SL) originally purified from the shoots of feverfew (Tanacetum parthenium), has shown potent anticancer and anti-inflammatory activities. It is currently being tested in cancer clinical trials. Structure-activity relationship (SAR) studies of parthenolide revealed key chemical properties required for biological activities and epigenetic mechanisms, and led to the derivatization of an orally bioavailable analog, dimethylamino-parthenolide (DMAPT). Parthenolide is the first small molecule found to be selective against cancer stem cells (CSC), which it achieves by targeting specific signaling pathways and killing cancer from its roots. In this review, we highlight the exciting journey of parthenolide, from plant shoots to cancer roots.