Search

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.

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.

The molecular events involved in the development of bronchial premalignant lesions (PMLs), and their progression to lung squamous cell carcinoma, are not well understood. Prior work characterized lung PML molecular subtypes by identifying co-expressed gene modules associated with histologic severity and progression/persistence. The proliferative subtype was enriched for PMLs with dysplasia. Genes related to interferon signaling and antigen processing/presentation (Module 9) were decreased in progressive/persistent lesions within this subtype, suggesting early immune suppression is related to PML progression. However, the mechanisms that drive these alterations are unclear. We investigated the role of microRNAs (miRNAs) in regulating gene expression associated with PML outcomes. mRNA and miRNA were extracted and sequenced from longitudinally collected endobronchial biopsies from patients with PMLs (148 samples, 30 patients). miRNAs targeting each gene co-expression module were identified based on target gene enrichment and the degree of negative correlation between the miRNA expression and its targets. Expression association with outcome within the proliferative subtype was tested with a mixed effects model adjusting for batch and patient as random effect. Cell type specificity of miRNA expression was tested based on cell type specific sequencing data from FANTOM5 and cell type marker correlation analysis. Genes regulated by NLRC5 were identified with ChIPseq data from Ludigs et al. Target gene suppression was confirmed by transfecting SW900 cells with miR-149-5p. miR-149-5p level in PML biopsies was examined by miRNA in situ hybridization (miR-ISH). miR-149-5p is identified as potential regulator of Module 9 gene expression and is significantly up-regulated in progressive/persistent PMLs. Its expression is highly enriched in epithelial cells in FANTOM5 and positively correlates with basal cell markers within PMLs. Predicted targets of miR-149-5p are down-regulated in the progressive PMLs in both our and data from Merrick et al. MHC-I and related gene expressions are down-regulated in progressing/persistent PMLs. These genes are regulated by the transcriptional coactivator NLRC5 which is a predicted target gene of miR-149-5p. We find that overexpressing miR-149-5p in SW900 cells decreases the expression levels of both NLRC5 and NLRC5 regulated genes. Additionally, miR-ISH targeting miR-149-5p in proliferative biopsy samples confirm its expression in epithelium compartment and association with outcome. Our data suggest epithelial miR-149-5p might be a key regulator of gene expression contributing to PML progression. We hypothesize that by suppressing NLRC5, miR-149-5p inhibits MHC-I gene expression of epithelial cells, promoting early immune depletion and lesion progression. miR-149-5p might therefore be a therapeutic target for preventing PML progression. Citation Format: Boting Ning, Roxana M. Pfefferkorn, Gang Liu, Sherry Zhang, Hanqiao Liu, Christopher Stevenson, Sarah A. Mazzilli, Avrum E. Spira, Marc E. Lenburg, Jennifer E. Beane. The role of epithelial miR-149 in immune modulation and progression of bronchial premalignant lesions [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 1488.

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.