Search

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

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.

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.

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.