Search

Supplementary Fig S8 and legend

Supplementary Figures S10-S20 and figure legends

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

Supplementary Methods

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

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.

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.

Soluble guanylyl cyclase (sGC) is the recognized receptor for the nitric oxide (NO) messenger molecule. Binding of NO to the heme moiety of sGC activates the enzyme several hundred fold. In addition to NO several NO-independent allosteric regulators can activate the enzyme in heme-dependent (YC-1, BAY 41–2272) or heme-independent (BAY 58–2667) fashions. One of the early recognized NO-independent activators of sGC is protoporphyrin IX, a tetrapyrrole compound that actively substitutes for heme and binds to the heme-deficient sGC.

Although soluble guanylyl cyclase (sGC) functions in an environment in which O(2), NO, and CO are potential ligands for its heme moiety, the enzyme displays a high affinity for only its physiological ligand, NO, but has a limited affinity for CO and no affinity for O(2). Recent studies of a truncated version of the sGC beta(1)-subunit containing the heme-binding domain (Boon, E. M., Huang, S H., and Marletta, M. A. (2005) Nat. Chem. Biol., 1, 53-59) showed that introduction of the hydrogen-bonding tyrosine into the distal heme pocket changes the ligand specificity of the heme moiety and results in an oxygen-binding sGC. The hypothesis that the absence of hydrogen-bonding residues in the distal heme pocket is sufficient to provide oxygen discrimination by sGC was put forward. We tested this hypothesis in a context of a complete sGC heterodimer containing both the intact alpha(1)- and beta(1)-subunits. We found that the I145Y substitution in the full-length beta-subunit of the sGC heterodimer did not produce an oxygen-binding enzyme. However, this substitution impeded the association of NO and destabilized the NO.heme complex. The tyrosine in the distal heme pocket also impeded both the binding and dissociation of the CO ligand. We propose that the mechanism of oxygen exclusion by sGC not only involves the lack of hydrogen bonding in the distal heme pocket, but also depends on structural elements from other domains of sGC.