Enhancing NSCLC recurrence prediction with PET/CT habitat imaging, ctDNA, and integrative radiogenomics-blood insights

While we recognize the prognostic importance of clinicopathological measures and circulating tumor DNA (ctDNA), the independent contribution of quantitative image markers to prognosis in non-small cell lung cancer (NSCLC) remains underexplored. In our multi-institutional study of 394 NSCLC patients, we utilize pre-treatment computed tomography (CT) and 18F-fluorodeoxyglucose positron emission tomography (FDG-PET) to establish a habitat imaging framework for assessing regional heterogeneity within individual tumors. This framework identifies three PET/CT subtypes, which maintain prognostic value after adjusting for clinicopathologic risk factors including tumor volume. Additionally, these subtypes complement ctDNA in predicting disease recurrence. Radiogenomics analysis unveil the molecular underpinnings of these imaging subtypes, highlighting downregulation in interferon alpha and gamma pathways in the high-risk subtype. In summary, our study demonstrates that these habitat imaging subtypes effectively stratify NSCLC patients based on their risk levels for disease recurrence after initial curative surgery or radiotherapy, providing valuable insights for personalized treatment approaches.

Synthetic PET from CT improves diagnosis and prognosis for lung cancer: Proof of concept

[18F]Fluorodeoxyglucose positron emission tomography (FDG-PET) and computed tomography (CT) are indispensable components in modern medicine. Although PET can provide additional diagnostic value, it is costly and not universally accessible, particularly in low-income countries. To bridge this gap, we have developed a conditional generative adversarial network pipeline that can produce FDG-PET from diagnostic CT scans based on multi-center multi-modal lung cancer datasets (n = 1,478). Synthetic PET images are validated across imaging, biological, and clinical aspects. Radiologists confirm comparable imaging quality and tumor contrast between synthetic and actual PET scans. Radiogenomics analysis further proves that the dysregulated cancer hallmark pathways of synthetic PET are consistent with actual PET. We also demonstrate the clinical values of synthetic PET in improving lung cancer diagnosis, staging, risk prediction, and prognosis. Taken together, this proof-of-concept study testifies to the feasibility of applying deep learning to obtain high-fidelity PET translated from CT.

Circulating microRNA Panel for Prediction of Recurrence and Survival in Early-Stage Lung Adenocarcinoma

Early-stage lung adenocarcinoma (LUAD) patients remain at substantial risk for recurrence and disease-related death, highlighting the unmet need of biomarkers for the assessment and identification of those in an early stage who would likely benefit from adjuvant chemotherapy. To identify circulating miRNAs useful for predicting recurrence in early-stage LUAD, we performed miRNA microarray analysis with pools of pretreatment plasma samples from patients with stage I LUAD who developed recurrence or remained recurrence-free during the follow-up period. Subsequent validation in 85 patients with stage I LUAD resulted in the development of a circulating miRNA panel comprising miR-23a-3p, miR-320c, and miR-125b-5p and yielding an area under the curve (AUC) of 0.776 in predicting recurrence. Furthermore, the three-miRNA panel yielded an AUC of 0.804, with a sensitivity of 45.8% at 95% specificity in the independent test set of 57 stage I and II LUAD patients. The miRNA panel score was a significant and independent factor for predicting disease-free survival (p < 0.001, hazard ratio [HR] = 1.64, 95% confidence interval [CI] = 1.51-4.22) and overall survival (p = 0.001, HR = 1.51, 95% CI = 1.17-1.94). This circulating miRNA panel is a useful noninvasive tool to stratify early-stage LUAD patients and determine an appropriate treatment plan with maximal efficacy.

Deep Learning-Based H-Score Quantification of Immunohistochemistry-Stained Images

Immunohistochemistry (IHC) is a well-established and commonly used staining method for clinical diagnosis and biomedical research. In most IHC images, the target protein is conjugated with a specific antibody and stained using diaminobenzidine (DAB), resulting in a brown coloration, whereas hematoxylin serves as a blue counterstain for cell nuclei. The protein expression level is quantified through the H-score, calculated from DAB staining intensity within the target cell region. Traditionally, this process requires evaluation by 2 expert pathologists, which is both time consuming and subjective. To enhance the efficiency and accuracy of this process, we have developed an automatic algorithm for quantifying the H-score of IHC images. To characterize protein expression in specific cell regions, a deep learning model for region recognition was trained based on hematoxylin staining only, achieving pixel accuracy for each class ranging from 0.92 to 0.99. Within the desired area, the algorithm categorizes DAB intensity of each pixel as negative, weak, moderate, or strong staining and calculates the final H-score based on the percentage of each intensity category. Overall, this algorithm takes an IHC image as input and directly outputs the H-score within a few seconds, significantly enhancing the speed of IHC image analysis. This automated tool provides H-score quantification with precision and consistency comparable to experienced pathologists but at a significantly reduced cost during IHC diagnostic workups. It holds significant potential to advance biomedical research reliant on IHC staining for protein expression quantification.

A distinct neurogenomic response to a trade-off between social challenge and opportunity in male sticklebacks (Gasterosteus aculeatus)

Animals frequently make adaptive decisions about what to prioritize when faced with multiple, competing demands simultaneously. However, the proximate mechanisms of decision-making in the face of competing demands are not well understood. We explored this question using brain transcriptomics in a classic model system: threespined sticklebacks, where males face conflict between courtship and territorial defence. We characterized the behaviour and brain gene expression profiles of males confronted by a trade-off between courtship and territorial defence by comparing them to males not confronted by this trade-off. When faced with the trade-off, males behaviourally prioritized defence over courtship, and this decision was reflected in their brain gene expression profiles. A distinct set of genes and biological processes was recruited in the brain when males faced a trade-off and these responses were largely non-overlapping across two brain regions. Combined, these results raise new questions about the interplay between the neural and molecular mechanisms involved in decision-making.

Automated Cellular-Level Dual Global Fusion of Whole-Slide Imaging for Lung Adenocarcinoma Prognosis

Histopathologic whole-slide images (WSI) are generally considered the gold standard for cancer diagnosis and prognosis. Survival prediction based on WSI has recently attracted substantial attention. Nevertheless, it remains a central challenge owing to the inherent difficulties of predicting patient prognosis and effectively extracting informative survival-specific representations from WSI with highly compounded gigapixels. In this study, we present a fully automated cellular-level dual global fusion pipeline for survival prediction. Specifically, the proposed method first describes the composition of different cell populations on WSI. Then, it generates dimension-reduced WSI-embedded maps, allowing for efficient investigation of the tumor microenvironment. In addition, we introduce a novel dual global fusion network to incorporate global and inter-patch features of cell distribution, which enables the sufficient fusion of different types and locations of cells. We further validate the proposed pipeline using The Cancer Genome Atlas lung adenocarcinoma dataset. Our model achieves a C-index of 0.675 (±0.05) in the five-fold cross-validation setting and surpasses comparable methods. Further, we extensively analyze embedded map features and survival probabilities. These experimental results manifest the potential of our proposed pipeline for applications using WSI in lung adenocarcinoma and other malignancies.

Immune cellular patterns of distribution affect outcomes of patients with non-small cell lung cancer

Studying the cellular geographic distribution in non-small cell lung cancer is essential to understand the roles of cell populations in this type of tumor. In this study, we characterize the spatial cellular distribution of immune cell populations using 23 makers placed in five multiplex immunofluorescence panels and their associations with clinicopathologic variables and outcomes. Our results demonstrate two cellular distribution patterns-an unmixed pattern mostly related to immunoprotective cells and a mixed pattern mostly related to immunosuppressive cells. Distance analysis shows that T-cells expressing immune checkpoints are closer to malignant cells than other cells. Combining the cellular distribution patterns with cellular distances, we can identify four groups related to inflamed and not-inflamed tumors. Cellular distribution patterns and distance are associated with survival in univariate and multivariable analyses. Spatial distribution is a tool to better understand the tumor microenvironment, predict outcomes, and may can help select therapeutic interventions.

Abstract 3744: A circulating microRNA panel predicts recurrence and survival in early-stage lung adenocarcinoma

Background: Early-stage lung adenocarcinoma (LUAD) patients have substantial risk for recurrence and disease-related death. Cisplatin-based adjuvant chemotherapy remains the standard for care of LUAD patients who have undergone surgical resection with a high risk of recurrence. However, adjuvant chemotherapy is associated with increased risk of toxicity including chemotherapy-related death, with only a modest survival benefit. Therefore, there is an unmet need of biomarkers for assessment and identification of those in an early stage who would likely benefit from adjuvant chemotherapy. Circulating miRNAs are stably present in blood and potentially reflect different expressions in cancerous and non-cancerous tissues, making them attractive biomarkers. The purpose of this study was to identify circulating miRNAs useful for predicting recurrence in early-stage LUAD.

Materials and Methods: miRNA microarray analysis was performed with pooled pretreatment plasma samples from stage I LUAD patients who developed recurrence within two years after curative surgery or remained recurrence free over a six-year follow-up period, as well as from healthy controls. miRNA biomarker candidates were assayed in two independent plasma sample sets from 85 stage I LUAD (validation set) and from 57 stage I and II LUAD (test set) patients.

Results: Based on miRNA microarray data and previous reports, predictive performance of miR-23a-3p, miR-23b-3p, miR-191-5p, miR-185-5p, miR-151a-3p, miR-320c, miR-21-5p, miR-125b-5p, miR-30d-5p, and miR-197-3p was evaluated in the validation set. Plasma levels of miR-23a-3p, miR-185-5p, miR-320c, miR-21-5p, miR-125b-5p, miR-30d-5p, and miR-197-3p were significantly higher in those with recurrence as compared to those without. A miRNA panel comprised of miR-23a-3p, miR-320c, and miR-125b-5p was developed based on a logistic regression, with yielding an AUC of 0.776 (95% confidence interval [CI] = 0.660 to 0.893). The three-miRNA panel with fixed coefficients yielded an AUC of 0.804 (95% CI = 0.688 to 0.920) with a sensitivity of 45.8% at 95% specificity in the test set. The miRNA panel score was a significant and independent factor for predicting disease-free survival (DFS; P &lt; 0.001, HR = 1.64, 95% CI = 1.51-4.22) and overall survival (OS; P = 0.001, HR = 1.51, 95% CI = 1.17-1.94).

Conclusion: This circulating miRNA panel may serve as a noninvasive blood test for predicting DFS and OS in early-stage LUAD patients. Our findings provide rationale for further investigation to stratify early-stage LUAD patients using blood-based biomarkers to increase the ability to provide more personalized care.

Citation Format: Mei-Chee Tai, Leonidas E. Bantis, Gargy Parhy, Taketo Kato, Ichidai Tanaka, Chi-Wan Chow, Junya Fujimoto, Carmen Behrens, Tetsunari Hase, Koji Kawaguchi, Johannes F. Fahrmann, Edwin Ostrin, Kohei Yokoi, Toyofumi F. Chen-Yoshikawa, Yoshinori Hasegawa, Samir M. Hanash, Ignacio I. Wistuba, Ayumu Taguchi. A circulating microRNA panel predicts recurrence and survival in early-stage 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 3744.

Abstract 4460: Spatial profiling of immune biomarkers in resected treatment-naïve early stage lung adenocarcinoma

Background: Recently, neoadjuvant immunotherapy plus chemotherapy has been approved for treatment of resectable non-small cell lung carcinoma (NSCLC). Defining the immune landscape of these tumors and its spatial distribution will help to understand lung cancer biology. Here, we analyzed the distribution of immune-related biomarkers in tumor-defined regions and its associations with clinicopathological variables in resected lung adenocarcinomas using high-plex profiling approaches.

Methodology: Thirty-three FFPE tumor tissues from surgically resected treatment-naïve lung adenocarcinoma stage I/II were used to construct a tissue microarray from the MD Anderson ICON cohort. We used three 1-mm core per patient [2 from central tumor (CT), and 1 from invasive margin (IM)] and performed the GeoMx Digital Spatial Profiling protein protocol to assess 49 immune biomarkers. Pancytokeratin (panCK; epithelial), CD45 (immune) and SYTO 13 (nuclear) were utilized as morphology biomarkers. Regions of interests were placed in cores containing tumor, and segmented in ‘”Tumor (Tu)” (PanCK+) and the “tumor microenvironment (TME)” (PanCK-). Digital counts were normalized using background correction. Statistical analysis was performed using linear mixed model. A p value equal or less than 0.05 was considered significant.

Results: We first compared the relative counts of immune biomarkers in the TME in CT and IM. IM had higher CD3, CD8, CD45RO, as well as CD163 and STING (P ranges 0.006 to 0.035), while CT had higher PD1 (P 0.025). Then we analyzed differential biomarker expression by sex and smoking status. Females had higher counts of Immune related biomarkers: CD45, CD3, CD20, immune checkpoints: PD-L1, VISTA, CTLA4, LAG3, and ICOS, and myeloid: CD68, CD11c, CD163, and B2M (P ranges 0.0005 to 0.046). Smokers had higher counts of CD66b (P 0.007) and B2M (P 0.039) while never smokers had higher counts of HLA-DR, CD34, FoxP3, OX40L, Tim-3, and B7-H3 (P ranges 0.001 to 0.049). Finally, we analyzed Tu segments. IM had higher CD66b, VISTA, CD163, OX40L, HLA-DR, GZMB, STING, and CD8 (P ranges 0.007 to 0.045) than CT. Female patients had higher CD45, CD68, CD11c and CD163 (P ranges 0.008 to 0.045), and males had higher SMA (P 0.005). Smokers had higher CD66b (P 0.007), B2-microglobulin (P 0.043), and never smokers had higher HLA-DR, STING, CD34, CD44, FoxP3 and CD25 (P ranges 0.006 to 0.043).

Conclusions: In this study, biomarker analysis of treatment-naïve adenocarcinoma in CT areas and IM indicates a higher immune response in the IM and presence of inhibitory signaling inside the tumor. Our data also showed that tumors from females have higher immune response than tumors from males, which is concordant with previous studies. Distinct profiling by smoking status was also observed. Further analysis of gene expression analysis of this set is ongoing.

Citation Format: Sharia D. Hernandez, Wei Lu, Alejandra G. Serrano, Claudio J. Arrechedera, Beatriz Sanchez-Espiridion, Nejla Ozirmak, Max Molina, Larisa Kostousov, Sean Barnes, Khaja Khan, Ximing Tang, Junya Fujimoto, Edwin R. Parra, Gabriela Raso, Stephanie T. Schmidt, Carmen Behrens, John Heymach, Jianjun Zhang, Ken Chen, Boris Sepesi, Tina Cascone, Don Gibbons, Ignacio I. Wistuba, Cara Haymaker, Luisa M. Solis. Spatial profiling of immune biomarkers in resected treatment-naïve early stage 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 4460.

Features of tumor-microenvironment images predict targeted therapy survival benefit in patients with EGFR-mutant lung cancer

Tyrosine kinase inhibitors (TKIs) targeting epidermal growth factor receptor (EGFR) are effective for many patients with lung cancer with EGFR mutations. However, not all patients are responsive to EGFR TKIs, including even those harboring EGFR-sensitizing mutations. In this study, we quantified the cells and cellular interaction features of the tumor microenvironment (TME) using routine H&E-stained biopsy sections. These TME features were used to develop a prediction model for survival benefit from EGFR TKI therapy in patients with lung adenocarcinoma and EGFR-sensitizing mutations in the Lung Cancer Mutation Consortium 1 (LCMC1) and validated in an independent LCMC2 cohort. In the validation data set, EGFR TKI treatment prolonged survival in the predicted-to-benefit group but not in the predicted-not-to-benefit group. Among patients treated with EGFR TKIs, the predicted-to-benefit group had prolonged survival outcomes compared with the predicted not-to-benefit group. The EGFR TKI survival benefit positively correlated with tumor-tumor interaction image features and negatively correlated with tumor-stroma interaction. Moreover, the tumor-stroma interaction was associated with higher activation of the hepatocyte growth factor/MET-mediated PI3K/AKT signaling pathway and epithelial-mesenchymal transition process, supporting the hypothesis of fibroblast-involved resistance to EGFR TKI treatment.

Immunogenomic intertumor heterogeneity across primary and metastatic sites in a patient with lung adenocarcinoma

Background

Lung cancer is the leading cause of cancer death, partially owing to its extensive heterogeneity. The analysis of intertumor heterogeneity has been limited by an inability to concurrently obtain tissue from synchronous metastases unaltered by multiple prior lines of therapy.

Methods

In order to study the relationship between genomic, epigenomic and T cell repertoire heterogeneity in a rare autopsy case from a 32-year-old female never-smoker with left lung primary late-stage lung adenocarcinoma (LUAD), we did whole-exome sequencing (WES), DNA methylation and T cell receptor (TCR) sequencing to characterize the immunogenomic landscape of one primary and 19 synchronous metastatic tumors.

Results

We observed heterogeneous mutation, methylation, and T cell patterns across distinct metastases. Only TP53 mutation was detected in all tumors suggesting an early event while other cancer gene mutations were later events which may have followed subclonal diversification. A set of prevalent T cell clonotypes were completely excluded from left-side thoracic tumors indicating distinct T cell repertoire profiles between left-side and non left-side thoracic tumors. Though a limited number of predicted neoantigens were shared, these were associated with homology of the T cell repertoire across metastases. Lastly, ratio of methylated neoantigen coding mutations was negatively associated with T-cell density, richness and clonality, suggesting neoantigen methylation may partially drive immunosuppression.

Conclusions

Our study demonstrates heterogeneous genomic and T cell profiles across synchronous metastases and how restriction of unique T cell clonotypes within an individual may differentially shape the genomic and epigenomic landscapes of synchronous lung metastases.

Immunogenomic profiling of lung adenocarcinoma reveals poorly differentiated tumors are associated with an immunogenic tumor microenvironment

OBJECTIVES

Pathologists have routinely observed distinct histologic patterns of growth in early-stage lung adenocarcinoma (LUAD), which have been suggested to be associated with prognosis. Herein, we investigated the relationship between LUAD patterns of growth, as defined by the updated international association for the study of lung cancer (IASLC) grading criteria, and differences in the tumor immune microenvironment to identify predictors of response to immunotherapy.

METHODS

174 resected stage I-III LUAD tumors were classified by histologic pattern of growth (i.e. solid, micropapillary, acinar, papillary, and lepidic) and then grouped as well differentiated, moderately differentiated, and poorly differentiated. Comprehensive multiplatform analysis including whole exome sequencing, gene expression profiling, immunohistochemistry, CIBERSORT, and T-cell receptor sequencing was performed and groups were compared for differences in genomic drivers, immune cell infiltrate, clonality, and survival. Finally, multivariate analysis was performed adjusting for pathologic stage and smoking status.

RESULTS

Poorly differentiated tumors demonstrated a strong association with smoking relative to moderately differentiated or well differentiated tumors. However, unlike in prior reports, poorly differentiated tumors were not associated with a worse survival after curative-intent resection. Genomic analysis revealed that poorly differentiated tumors are associated with high tumor mutation burden but showed no association with oncogenic drivers. Immune analyses revealed that poorly differentiated tumors are associated with increased T-cell clonality, expression of PD-L1, and infiltration by cytotoxic CD8 T-cells, activated CD4 T-cells, and pro-inflammatory (M1) macrophages. Finally, multivariate analysis controlling for stage and smoking status confirmed independence of immune differences between IASLC grade groups.

CONCLUSIONS

Poorly differentiated tumors, as defined by the updated IASLC grading criteria, are associated with a distinct immunogenic tumor microenvironment that predicts for therapeutic response to immune agents, including checkpoint inhibitors, and should be included in the clinical trial design of immunotherapy studies in early-stage lung adenocarcinoma.

Targeting de novo lipogenesis and the Lands cycle induces ferroptosis in KRAS-mutant lung cancer

Mutant KRAS (KM), the most common oncogene in lung cancer (LC), regulates fatty acid (FA) metabolism. However, the role of FA in LC tumorigenesis is still not sufficiently characterized. Here, we show that KMLC has a specific lipid profile, with high triacylglycerides and phosphatidylcholines (PC). We demonstrate that FASN, the rate-limiting enzyme in FA synthesis, while being dispensable in EGFR-mutant or wild-type KRAS LC, is required for the viability of KMLC cells. Integrating lipidomic, transcriptomic and functional analyses, we demonstrate that FASN provides saturated and monounsaturated FA to the Lands cycle, the process remodeling oxidized phospholipids, such as PC. Accordingly, blocking either FASN or the Lands cycle in KMLC, promotes ferroptosis, a reactive oxygen species (ROS)- and iron-dependent cell death, characterized by the intracellular accumulation of oxidation-prone PC. Our work indicates that KM dictates a dependency on newly synthesized FA to escape ferroptosis, establishing a targetable vulnerability in KMLC.

Shared Nearest Neighbors Approach and Interactive Browser for Network Analysis of a Comprehensive Non-Small-Cell Lung Cancer Data Set

PURPOSE

Advances in biological measurement technologies are enabling large-scale studies of patient cohorts across multiple omics platforms. Holistic analysis of these data can generate actionable insights for translational research and necessitate new approaches for data integration and mining.

METHODS

We present a novel approach for integrating data across platforms on the basis of the shared nearest neighbors algorithm and use it to create a network of multiplatform data from the immunogenomic profiling of non-small-cell lung cancer project.

RESULTS

Benchmarking demonstrates that the shared nearest neighbors-based network approach outperforms a traditional gene-gene network in capturing established interactions while providing new ones on the basis of the interplay between measurements from different platforms. When used to examine patient characteristics of interest, our approach provided signatures associated with and new leads related to recurrence and TP53 oncogenotype.

CONCLUSION

The network developed offers an unprecedented, holistic view into immunogenomic profiling of non-small-cell lung cancer, which can be explored through the accompanying interactive browser that we built.

Abstract 6152: Distinct immune gene programs associated with host tumor immunity, neoadjuvant chemotherapy and chemoimmunotherapy in resectable NSCLC

Background: Our understanding of the immunopathology of early-stage NSCLC is still limited. While neoadjuvant immunotherapeutic strategies have recently shown anti-tumor effects in resectable NSCLC, their mechanisms remain inadequately understood. Here, we explore immune programs that inform of tumor immunity and response to neoadjuvant chemotherapy and chemoimmunotherapy in localized NSCLC.

Methods: Targeted immune gene sequencing using the HTG Precision Immuno-Oncology panel was performed in localized NSCLCs from three cohorts based on treatment: naïve (n=190), neoadjuvant chemotherapy (n=38) and neoadjuvant chemoimmunotherapy (n=21). Three tumor immune microenvironment (TIME) phenotypes (inflamed, cold, excluded) were derived based on CD8+ T cell infiltration. Signatures of immune cell abundance and immune genes were statistically compared based on tumoral PD-L1 expression, immune phenotypes, associated with pathological response, and were cross-compared across the three cohorts.

Results: PD-L1 positive tumors exhibited increased signature scores for various lymphoid and myeloid cell subsets (both, p&lt;0.05). TIME phenotypes exhibited disparate frequencies by stage, PD-L1 expression, and mutational burden. Inflamed NSCLCs displayed overall significantly heightened levels of immune signatures with the excluded group representing an intermediate state. A signature of cytotoxic T cells was associated with favorable survival in neoadjuvant chemotherapy-treated NSCLCs (p&lt;0.05). Major pathological response to chemoimmunotherapy was positively associated with CD8 T cells (p&lt;0.05) and Th1 cells were significantly reduced post-chemoimmunotherapy (p&lt;0.001). Among the three cohorts, chemoimmunotherapy-treated NSCLCs exhibited highest scores for various immune cell subsets including T effector and B cells (both, p&lt;0.05).

Conclusions: Our findings highlight immune gene programs that may underlie host tumor immunity and response to neoadjuvant chemotherapy and chemoimmunotherapy in early-stage NSCLC.

Citation Format: Pedro Rocha, Jiexin Zhang, Raquel Laza-Briviesca, Alberto Cruz-Bermúdez, Katsuhiro Yoshimura, Carmen Behrens, Apar Pataer, Edwin Parra-Cuentas, Cara Haymaker, Junya Fujimoto, Stephen Swisher, John Heymach, Don L. Gibbons, J Jack Lee, Boris Sepesi, Tina Cascone, Luisa M. Solis, Mariano Provencio, Ignacio I. Wistuba, Humam Kadara. Distinct immune gene programs associated with host tumor immunity, neoadjuvant chemotherapy and chemoimmunotherapy in resectable NSCLC [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 6152.

Mutational Activation of the NRF2 Pathway Upregulates Kynureninase Resulting in Tumor Immunosuppression and Poor Outcome in Lung Adenocarcinoma

Activation of the NRF2 pathway through gain-of-function mutations or loss-of-function of its suppressor KEAP1 is a frequent finding in lung cancer. NRF2 activation has been reported to alter the tumor microenvironment. Here, we demonstrated that NRF2 alters tryptophan metabolism through the kynurenine pathway that is associated with a tumor-promoting, immune suppressed microenvironment. Specifically, proteomic profiles of 47 lung adenocarcinoma (LUAD) cell lines (11 KEAP1 mutant and 36 KEAP1 wild-type) revealed the tryptophan-kynurenine enzyme kynureninase (KYNU) as a top overexpressed protein associated with activated NRF2. The siRNA-mediated knockdown of NFE2L2, the gene encoding for NRF2, or activation of the NRF2 pathway through siRNA-mediated knockdown of KEAP1 or via chemical induction with the NRF2-activator CDDO-Me confirmed that NRF2 is a regulator of KYNU expression in LUAD. Metabolomic analyses confirmed KYNU to be enzymatically functional. Analysis of multiple independent gene expression datasets of LUAD, as well as a LUAD tumor microarray demonstrated that elevated KYNU was associated with immunosuppression, including potent induction of T-regulatory cells, increased levels of PD1 and PD-L1, and resulted in poorer survival. Our findings indicate a novel mechanism of NRF2 tumoral immunosuppression through upregulation of KYNU.

Distinct immune gene programs associated with host tumor immunity, neoadjuvant chemotherapy and chemoimmunotherapy in resectable NSCLC

PURPOSE

Our understanding of the immunopathology of resectable NSCLC is still limited. Here, we explore immune programs that inform of tumor immunity and response to neoadjuvant chemotherapy and chemoimmunotherapy in localized NSCLC.

EXPERIMENTAL DESIGN

Targeted immune gene sequencing using the HTG Precision Immuno-Oncology panel was performed in localized NSCLCs from three cohorts based on treatment: naïve (n=190), neoadjuvant chemotherapy (n=38) and neoadjuvant chemoimmunotherapy (n=21). Tumor immune microenvironment (TIME) phenotypes based on the location of CD8+ T cells (inflamed, cold, excluded), tumoral PD-L1 expression (&lt;1% and {greater than or equal to}1%), and tumor infiltrating lymphocytes (TILs). Immune programs and signatures were statistically analyzed based on tumoral PD-L1 expression, immune phenotypes, pathological response and were cross-compared across the three cohorts.

RESULTS

PD-L1 positive tumors exhibited increased signature scores for various lymphoid and myeloid cell subsets (p&lt;0.05). TIME phenotypes exhibited disparate frequencies by stage, PD-L1 expression, and mutational burden. Inflamed and PD-L1+/TILs+ NSCLCs displayed overall significantly heightened levels of immune signatures, with the excluded group representing an intermediate state. A cytotoxic T cell signature was associated with favorable survival in neoadjuvant chemotherapy-treated NSCLCs (p&lt;0.05). Pathological response to chemoimmunotherapy was positively associated with higher expression of genes involved in immune activation, chemotaxis, as well as T and NK cells (p&lt;0.05 for all). Among the three cohorts, chemoimmunotherapy-treated NSCLCs exhibited highest scores for various immune cell subsets including T effector and B cells (p&lt;0.05).

CONCLUSIONS

Our findings highlight immune gene programs that may underlie host tumor immunity and response to neoadjuvant chemotherapy and chemoimmunotherapy in resectable NSCLC.

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

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

AXL targeting restores PD-1 blockade sensitivity of STK11/LKB1 mutant NSCLC through expansion of TCF1+ CD8 T cells

Mutations in STK11/LKB1 in non-small cell lung cancer (NSCLC) are associated with poor patient responses to immune checkpoint blockade (ICB), and introduction of a Stk11/Lkb1 (L) mutation into murine lung adenocarcinomas driven by mutant Kras and Trp53 loss (KP) resulted in an ICB refractory syngeneic KPL tumor. Mechanistically this occurred because KPL mutant NSCLCs lacked TCF1-expressing CD8 T cells, a phenotype recapitulated in human STK11/LKB1 mutant NSCLCs. Systemic inhibition of Axl results in increased type I interferon secretion from dendritic cells that expanded tumor-associated TCF1+PD-1+CD8 T cells, restoring therapeutic response to PD-1 ICB in KPL tumors. This was observed in syngeneic immunocompetent mouse models and in humanized mice bearing STK11/LKB1 mutant NSCLC human tumor xenografts. NSCLC-affected individuals with identified STK11/LKB1 mutations receiving bemcentinib and pembrolizumab demonstrated objective clinical response to combination therapy. We conclude that AXL is a critical targetable driver of immune suppression in STK11/LKB1 mutant NSCLC.

The histologic phenotype of lung cancers is associated with transcriptomic features rather than genomic characteristics

Histology plays an essential role in therapeutic decision-making for lung cancer patients. However, the molecular determinants of lung cancer histology are largely unknown. We conduct whole-exome sequencing and microarray profiling on 19 micro-dissected tumor regions of different histologic subtypes from 9 patients with lung cancers of mixed histology. A median of 68.9% of point mutations and 83% of copy number aberrations are shared between different histologic components within the same tumors. Furthermore, different histologic components within the tumors demonstrate similar subclonal architecture. On the other hand, transcriptomic profiling reveals shared pathways between the same histologic subtypes from different patients, which is supported by the analyses of the transcriptomic data from 141 cell lines and 343 lung cancers of different histologic subtypes. These data derived from mixed histologic subtypes in the setting of identical genetic background and exposure history support that the histologic fate of lung cancer cells is associated with transcriptomic features rather than the genomic profiles in most tumors.

The molecular determinants of lung cancer histologic subtypes are not well understood. Here the authors analyze lung cancers of mixed histology and find that histologic subtypes are associated with transcriptomic features rather than genomic profiles in most tumors.

Abstract P009: A shared nearest neighbors approach for integrated, multi-platform networks and its application to the exploration of multiomics data from early-stage non-small cell lung cancers

Background: The ImmunogenomiC prOfiling of Non-small cell lung cancer (NSCLC) Project (ICON) represents an ambitious undertaking to comprehensively characterize immuno-genomic diversity in NSCLC across diverse platforms. The depth and breadth of this cohort presented a unique opportunity to develop a specialized method for multi-platform data integration and exploration, which can be broadly applied to forthcoming large-scale patient profiling studies. Such a holistic approach can unlock insights for therapeutic targets, biomarkers, and treatment plans by providing a more complete view of phenomena driving disease pathogenesis and evolution. Purpose: We developed a novel shared nearest neighbors (SNN) approach to create an integrated network of ICON’s multi-platform data and identified collections of closely related measurements within the resulting network tied to noteworthy patient characteristics, including recurrence and oncogenotype. Methods: The ICON dataset is derived from tumor and normal lung tissue samples collected from 150 patients at time of resection as well as blood samples collected then and at intervals during the year following. Tissue samples underwent RNA-sequencing (RNA-seq), whole exome sequencing, T-cell receptor sequencing, multiplex immunofluorescence for immune cells, and reverse phase protein array profiling; flow cytometry for immune cells was performed on tissue and blood samples. From these data, the ICON data network was built using an integrative approach based on the SNN algorithm in which genes were linked on the basis of their shared top correlates in orthogonal datasets. Results: The ICON data network currently includes over 20,000 genes linked by over 500,000 connections derived from correlations between RNA-seq and orthogonal platforms. We captured established associations between cancer-related genes and examined these along with new ones in the network. To do so, we used the InfoMap algorithm to extract more interpretable sub-networks, termed modules, from the ICON data network. Single sample gene set enrichment scores for each module were used in multivariate analysis to highlight modules linked to clinical characteristics of interest. As an example, we found modules significantly tied to disease recurrence. The most notable of these was strongly associated with metabolic pathways, and other modules associated with platelets and ion channels were also identified. The metabolic pathway module is being explored as a prognostic biomarker, underscoring the opportunites enabled by mining the network. Conclusions: Through the framework developed, we identified modules in the ICON data network significantly associated with important patient characteristics like recurrence and oncogenotype. We are validating the gene sets identified as potential biomarkers and are developing an interactive application to facilitate further mining of the network. Taken together, our SNN network-building approach enables the integration and exploration of patient data from diverse platforms.

Citation Format: Stephanie T. Schmidt, Neal Akhave, Alexandre Reuben, Tina Cascone, Jianhua Zhang, Jun Li, Junya Fujimoto, Lauren A. Byers, Beatriz Sanchez-Espiridion, Lixia Diao, Jing Wang, Lorenzo Federico, Marie-Andree Forget, Daniel J McGrail, Annikka Weissferdt, Shiaw-Yih Lin, Younghee Lee, Natalie Vokes, Carmen Behrens, Ignacio I. Wistuba, Andrew Futreal, Ara Vaporciyan, Boris Sepesi, John V. Heymach, Chantale Bernatchez, Cara Haymaker, Jianjun Zhang, Christopher A. Bristow, Timothy P. Heffernan, Marcelo V. Negrao, Don L. Gibbons. A shared nearest neighbors approach for integrated, multi-platform networks and its application to the exploration of multiomics data from early-stage non-small cell lung cancers [abstract]. In: Proceedings of the AACR-NCI-EORTC Virtual International Conference on Molecular Targets and Cancer Therapeutics; 2021 Oct 7-10. Philadelphia (PA): AACR; Mol Cancer Ther 2021;20(12 Suppl):Abstract nr P009.

Cold and heterogeneous T cell repertoire is associated with copy number aberrations and loss of immune genes in small-cell lung cancer

Small-cell lung cancer (SCLC) is speculated to harbor complex genomic intratumor heterogeneity (ITH) associated with high recurrence rate and suboptimal response to immunotherapy. Here, using multi-region whole exome/T cell receptor (TCR) sequencing as well as immunohistochemistry, we reveal a rather homogeneous mutational landscape but extremely cold and heterogeneous TCR repertoire in limited-stage SCLC tumors (LS-SCLCs). Compared to localized non-small cell lung cancers, LS-SCLCs have similar predicted neoantigen burden and genomic ITH, but significantly colder and more heterogeneous TCR repertoire associated with higher chromosomal copy number aberration (CNA) burden. Furthermore, copy number loss of IFN-γ pathway genes is frequently observed and positively correlates with CNA burden. Higher mutational burden, higher T cell infiltration and positive PD-L1 expression are associated with longer overall survival (OS), while higher CNA burden is associated with shorter OS in patients with LS-SCLC.

602 AXL targeting with bemcentinb restores PD-1 blockade sensitivity of STK11/LKB1 mutant NSCLC through innate immune cell mediated expansion of TCF1+ CD8 T cells

Background

Mutations in tumor suppressor STK11/LKB1 are associated with negative predictive and prognostic impact in NSCLC patients receiving immune checkpoint inhibitors (CPI) in several published cohorts, although there have been some conflicting reports on the association of such mutations with patient outcomes in this setting [1–9]. STK11/LKB1 tumors are characterized by a suppressive tumor micro-environment devoid of cytotoxic T cells, and we hypothesized that targeting the receptor tyrosine kinase AXL, a known driver of an innate immune suppressive microenvironment, would restore sensitivity to PD-1 in syngeneic pre-clinical models as well as in patients harboring STK11/LKB1 mutated NSCLC.

Methods

Stk11/Lkb1 (L) mutation was introduced by CRISPR technology into murine lung adenocarcinomas driven by mutant Kras and Trp53 loss (KP). Sensitivity towards anti-PD-1 was evaluated in the absence and presence of the small molecule AXL inhibitor bemcentinib in the KPL model and in a human NSCLC xenograft model carrying a STK11/LKB1 mutation. The immune tumor landscape was mapped following introduction of the Stk11/Lkb1 mutation and therapeutic intervention with anti-PD-1/pembrolizumab and bemcentinib. FFPE fine-needle aspirate biopsies of target lesions were acquired from patients at screening immediately prior to enrollment in BGBC008, a PhII single-arm, 2-stage study with bemcentinib (200mg/d) and pembrolizumab (200 mg/q3wk) for previously-treated stage IV lung adenocarcinoma patients who were CPI naïve or CPI refractory. Patients were assessed for response according to RECIST1.1 criteria at scheduled scan intervals.

Results

Introduction of a STK11/LKB1 mutation into murine lung adenocarcinomas driven by mutant Kras and Trp53 loss resulted in a PD-1 refractory syngeneic KPL tumor. Mechanistically this occurred because KPL mutant NSCLCs lacked TCF1-expressing CD8 T cells, a phenotype that was recapitulated in human STK11/LKB1 mutant NSCLCs. Systemic inhibition of AXL with bemcentinib resulted in increased type I interferon secretion from dendritic cells resulting in expansion of tumor-associated TCF1+PD-1+CD8 T cells and restored therapeutic response to PD-1. This effect was observed in a syngeneic immunocompetent mouse model and in humanized mice bearing STK11/LKB1 mutant NSCLC human tumor xenografts.In an ongoing clinical trial (NCT03184571), 3 evaluable NSCLC patients with identified STK11/LKB1 mutations demonstrated objective clinical response/clinical benefit to the combination of AXL inhibitor bemcentinib and pembrolizumab

Conclusions

In these models, AXL is a critical targetable driver of immune suppression in STK11/LKB1 mutant NSCLC contributing to CPI resistance. Our results show that inhibition of AXL rescues this deficit and represents a new clinical strategy in combination with anti-PD-1 therapy in NSCLC patients carrying a STK11/LKB1 mutation

Acknowledgements

The authors would like to thank all patients and their caretakers for participating in this trial.

Trial Registration

Patients treated with bemcentinib and pembrolizumab combination therapy were enrolled in the BGBC008 clinical trial (BerGenBio ASA and Merck & Co., Inc., Kenilworth NJ, USA, NCT03184571)

References

Gu M, Xu T, Chang P. KRAS/LKB1 and KRAS/TP53 co-mutations create divergent immune signatures in lung adenocarcinomas. Ther Adv Med Oncol. 2021;13:17588359211006950.

Cho BC, Lopes G, Kowalski DM. Relationship between STK11 and KEAP1 mutational status and efficacy in KEYNOTE-042: pembrolizumab monotherapy as first-line therapy for PD-L1 positive advanced NSCLC. Cancer Res. 2020;80(16 Supplement):CT084.

Aredo JV, Padda SK, Kunder CA. Impact of KRAS mutation subtype and concurrent pathogenic mutations on non-small cell lung cancer outcomes. Lung Cancer. 2019;133:144–150.

Kwack WG, Shin SY, Lee SH. Primary Resistance to Immune Checkpoint Blockade in an STK11/TP53/KRAS-Mutant Lung Adenocarcinoma with High PD-L1 Expression. Oncol Targets Ther. 2020;13:8901–8905.

Shire NJ, Klein AB, Golozar A. STK11 (LKB1) mutations in metastatic NSCLC: Prognostic value in the real world. PLoS One. 2020;15(9):e0238358. 6. Skoulidis F, Goldberg ME, Greenawalt DM. STK11/LKB1 Mutations and PD-1 Inhibitor Resistance in KRAS-Mutant Lung Adenocarcinoma. Cancer Discov. 2018;8(7):822–835. 7. Wang H, Guo J, Shang X. Less immune cell infiltration and worse prognosis after immunotherapy for patients with lung adenocarcinoma who harbored STK11 mutation. Int. Immunopharmacol. 2020;84:106574. 8. Kitajima S, Ivanova E, Gou S. Suppression of STING Associated with LKB1 Loss in KRAS-Driven Lung Cancer. Cancer Discov. 2019;9(1):34–45. 9. Mograbi B, Heeke S, Hofman P. The Importance of STK11/LKB1 Assessment in Non-Small Cell Lung Carcinomas. Diagnostics (Basel). 2021;11(2):196.

Ethics Approval

This study was approved by the following ethical committees: Use of human cord blood: UT Southwestern (UTSW) Parkland Hospital, STU 112010-047Animal studies: UTSW Medical Center, Institutional Animal Care and Use Committee, APN 2015-100921Clinical study: London Bridge Research Ethics Committee (UK): 17/LO/0418; REC-South East (Norway): 2017/473; Drug Research Ethics Committee of the University Hospital Clinic of Barcelona (Spain): BGBC008/MK-3475_PN-531; Medical College of Wisconsin Institutional Review Board #4 (USA): PRO00029453

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.

SRGN-Triggered Aggressive and Immunosuppressive Phenotype in a Subset of TTF-1-Negative Lung Adenocarcinomas

BACKGROUND

About 20% of lung adenocarcinoma (LUAD) is negative for the lineage-specific oncogene Thyroid transcription factor 1 (TTF-1) and exhibits worse clinical outcome with a low frequency of actionable genomic alterations. To identify molecular features associated with TTF-1-negative LUAD, we compared the transcriptomic and proteomic profiles of LUAD cell lines. SRGN, a chondroitin sulfate proteoglycan Serglycin, was identified as a markedly overexpressed gene in TTF-1-negative LUAD. We therefore investigated the roles and regulation of SRGN in TTF-1-negative LUAD.

METHODS

Proteomic and metabolomic analyses of 41 LUAD cell lines were done using mass spectrometry. The function of SRGN was investigated in 3 TTF-1-negative and 4 TTF-1-positive LUAD cell lines and in a syngeneic mouse model (n = 5 to 8 mice per group). Expression of SRGN in was evaluated in 94 and 105 surgically resected LUAD tumor specimens using immunohistochemistry. All statistical tests were two-sided.

RESULTS

SRGN was markedly overexpressed at mRNA and protein levels in TTF-1-negative LUAD cell lines (P < .001 for both mRNA and protein levels). Expression of SRGN in LUAD tumor tissue was associated with poor outcome (hazard ratio = 4.22, 95% confidential interval = 1.12 to 15.86; likelihood ratio test, P = .03), and with higher expression of Programmed cell death 1 ligand 1 (PD-L1) in tumor cells and higher infiltration of Programmed cell death protein 1 (PD-1)-positive lymphocytes. SRGN regulated expression of PD-L1, as well as proinflammatory cytokines including Interleukin-6 (IL-6), Interleukin-8 (IL-8), and C-X-C motif chemokine 1 (CXCL1) in LUAD cell lines, and increased migratory and invasive properties of LUAD cells and fibroblasts, and enhanced angiogenesis. SRGN was induced by DNA de-methylation resulting from Nicotinamide N-methyltransferase (NNMT)-mediated impairment of methionine metabolism.

CONCLUSION

Our findings suggest that SRGN plays a pivotal role in tumor-stromal interaction and reprogramming into an aggressive and immunosuppressive tumor microenvironment in TTF-1-negative LUAD.

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.

Immune evolution from preneoplasia to invasive lung adenocarcinomas and underlying molecular features

The mechanism by which anti-cancer immunity shapes early carcinogenesis of lung adenocarcinoma (ADC) is unknown. In this study, we characterize the immune contexture of invasive lung ADC and its precursors by transcriptomic immune profiling, T cell receptor (TCR) sequencing and multiplex immunofluorescence (mIF). Our results demonstrate that anti-tumor immunity evolved as a continuum from lung preneoplasia, to preinvasive ADC, minimally-invasive ADC and frankly invasive lung ADC with a gradually less effective and more intensively regulated immune response including down-regulation of immune-activation pathways, up-regulation of immunosuppressive pathways, lower infiltration of cytotoxic T cells (CTLs) and anti-tumor helper T cells (Th), higher infiltration of regulatory T cells (Tregs), decreased T cell clonality, and lower frequencies of top T cell clones in later-stages. Driver mutations, chromosomal copy number aberrations (CNAs) and aberrant DNA methylation may collectively impinge host immune responses and facilitate immune evasion, promoting the outgrowth of fit subclones in preneoplasia into dominant clones in invasive ADC.

Female Gender Predicts Augmented Immune Infiltration in Lung Adenocarcinoma

INTRODUCTION

Immune infiltration in lung adenocarcinomas (LUADs) has been associated with response to immune checkpoint inhibitors. Clinical features underlying differential responses of patients with LUADs to immunotherapy are not well understood. Here, we analyzed the association between LUAD immune infiltration and clinicopathologic variables.

MATERIALS AND METHODS

Intratumoral CD3, CD8, and CD68 cell densities (tumor-associated immune cells [TAICs]) were immunohistochemically assessed in 146 surgically resected LUADs. LUADs were classified into 2 groups, low and high TAICs, based on the median values of cell densities for CD3, CD8, and CD68. Somatic mutation burden and driver gene mutation status were analyzed in a subset of the cases (n = 92). We statistically analyzed the association between the TAIC groups and various clinicopathologic and molecular variables by using the χ²/Fisher and Wilcoxon sum tests and multivariable logistic regression models.

RESULTS

Patient gender, tumor size, and STK11 mutations were significantly associated with TAIC levels in LUAD. Female patients exhibited significantly elevated TAIC levels (P = .005) compared with male patients. Tumor size was inversely associated with TAIC levels (P = .012). STK11 mutated tumors were associated with lower TAICs (P = .008). Higher TAICs were consistently observed in female patients with LUADs after adjusting for stage, tumor size, and age. Multivariable regression models confirmed female gender as an independent variable associated with TAIC levels in LUAD (P = .0141).

CONCLUSION

Immune infiltration in LUADs was significantly higher in female patients, warranting further exploration into the association between this clinical variable and immunotherapeutic response in LUAD.

Characterization of the Immune Landscape of EGFR-Mutant NSCLC Identifies CD73/Adenosine Pathway as a Potential Therapeutic Target

INTRODUCTION

Lung adenocarcinomas harboring EGFR mutations do not respond to immune checkpoint blockade therapy and their EGFR wildtype counterpart. The mechanisms underlying this lack of clinical response have been investigated but remain incompletely understood.

METHODS

We analyzed three cohorts of resected lung adenocarcinomas (Profiling of Resistance Patterns of Oncogenic Signaling Pathways in Evaluation of Cancer of Thorax, Immune Genomic Profiling of NSCLC, and The Cancer Genome Atlas) and compared tumor immune microenvironment of EGFR-mutant tumors to EGFR wildtype tumors, to identify actionable regulators to target and potentially enhance the treatment response.

RESULTS

EGFR-mutant NSCLC exhibited low programmed death-ligand 1, low tumor mutational burden, decreased number of cytotoxic T cells, and low T cell receptor clonality, consistent with an immune-inert phenotype, though T cell expansion ex vivo was preserved. In an analysis of 75 immune checkpoint genes, the top up-regulated genes in the EGFR-mutant tumors (NT5E and ADORA1) belonged to the CD73/adenosine pathway. Single-cell analysis revealed that the tumor cell population expressed CD73, both in the treatment-naive and resistant tumors. Using coculture systems with EGFR-mutant NSCLC cells, T regulatory cell proportion was decreased with CD73 knockdown. In an immune-competent mouse model of EGFR-mutant lung cancer, the CD73/adenosine pathway was markedly up-regulated and CD73 blockade significantly inhibited tumor growth.

CONCLUSIONS

Our work revealed that EGFR-mutant NSCLC has an immune-inert phenotype. We identified the CD73/adenosine pathway as a potential therapeutic target for EGFR-mutant NSCLC.

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.

PI4KIIIβ is a therapeutic target in chromosome 1q-amplified lung adenocarcinoma

Heightened secretion of protumorigenic effector proteins is a feature of malignant cells. Yet, the molecular underpinnings and therapeutic implications of this feature remain unclear. Here, we identify a chromosome 1q region that is frequently amplified in diverse cancer types and encodes multiple regulators of secretory vesicle biogenesis and trafficking, including the Golgi-dedicated enzyme phosphatidylinositol (PI)-4-kinase IIIβ (PI4KIIIβ). Molecular, biochemical, and cell biological studies show that PI4KIIIβ-derived PI-4-phosphate (PI4P) synthesis enhances secretion and accelerates lung adenocarcinoma progression by activating Golgi phosphoprotein 3 (GOLPH3)-dependent vesicular release from the Golgi. PI4KIIIβ-dependent secreted factors maintain 1q-amplified cancer cell survival and influence prometastatic processes in the tumor microenvironment. Disruption of this functional circuitry in 1q-amplified cancer cells with selective PI4KIIIβ antagonists induces apoptosis and suppresses tumor growth and metastasis. These results support a model in which chromosome 1q amplifications create a dependency on PI4KIIIβ-dependent secretion for cancer cell survival and tumor progression.

Evolution of DNA methylome from precancerous lesions to invasive lung adenocarcinomas

The evolution of DNA methylome and methylation intra-tumor heterogeneity (ITH) during early carcinogenesis of lung adenocarcinoma has not been systematically studied. We perform reduced representation bisulfite sequencing of invasive lung adenocarcinoma and its precursors, atypical adenomatous hyperplasia, adenocarcinoma in situ and minimally invasive adenocarcinoma. We observe gradual increase of methylation aberrations and significantly higher level of methylation ITH in later-stage lesions. The phylogenetic patterns inferred from methylation aberrations resemble those based on somatic mutations suggesting parallel methylation and genetic evolution. De-convolution reveal higher ratio of T regulatory cells (Tregs) versus CD8 + T cells in later-stage diseases, implying progressive immunosuppression with neoplastic progression. Furthermore, increased global hypomethylation is associated with higher mutation burden, copy number variation burden and AI burden as well as higher Treg/CD8 ratio, highlighting the potential impact of methylation on chromosomal instability, mutagenesis and tumor immune microenvironment during early carcinogenesis of lung adenocarcinomas.

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.

Interplay between estrogen and Stat3/NF-κB-driven immunomodulation in lung cancer

K-ras mutant lung adenocarcinoma (LUAD) is the most common type of lung cancer, displays abysmal prognosis and is tightly linked to tumor-promoting inflammation, which is increasingly recognized as a target for therapeutic intervention. We have recently shown a gender-specific role for epithelial Stat3 signaling in the pathogenesis of K-ras mutant LUAD. The absence of epithelial Stat3 in male K-ras mutant mice (LR/Stat3Δ/Δ mice) promoted tumorigenesis and induced a nuclear factor-kappaB (NF-κB)-driven pro-tumor immune response while reducing tumorigenesis and enhancing anti-tumor immunity in female counterparts. In the present study, we manipulated estrogen and NF-κB signaling to study the mechanisms underlying this intriguing gender-disparity. In LR/Stat3Δ/Δ females, estrogen deprivation by bilateral oophorectomy resulted in higher tumor burden, an induction of NF-κB-driven immunosuppressive response, and reduced anti-tumor cytotoxicity, whereas estrogen replacement reversed these changes. On the other hand, exogenous estrogen in males successfully inhibited tumorigenesis, attenuated NF-κB-driven immunosuppression and boosted anti-tumor immunity. Mechanistically, genetic targeting of epithelial NF-κB activity resulted in reduced tumorigenesis and enhanced the anti-tumor immune response in LR/Stat3Δ/Δ males, but not females. Our data suggest that estrogen exerts a context-specific anti-tumor effect through inhibiting NF-κB-driven tumor-promoting inflammation and provide insights into developing novel personalized therapeutic strategies for K-ras mutant LUAD.

Multiomics profiling of primary lung cancers and distant metastases reveals immunosuppression as a common characteristic of tumor cells with metastatic plasticity

BACKGROUND

Metastasis is the primary cause of cancer mortality accounting for 90% of cancer deaths. Our understanding of the molecular mechanisms driving metastasis is rudimentary.

RESULTS

We perform whole exome sequencing (WES), RNA sequencing, methylation microarray, and immunohistochemistry (IHC) on 8 pairs of non-small cell lung cancer (NSCLC) primary tumors and matched distant metastases. Furthermore, we analyze published WES data from 35 primary NSCLC and metastasis pairs, and transcriptomic data from 4 autopsy cases with metastatic NSCLC and one metastatic lung cancer mouse model. The majority of somatic mutations are shared between primary tumors and paired distant metastases although mutational signatures suggest different mutagenesis processes in play before and after metastatic spread. Subclonal analysis reveals evidence of monoclonal seeding in 41 of 42 patients. Pathway analysis of transcriptomic data reveals that downregulated pathways in metastases are mainly immune-related. Further deconvolution analysis reveals significantly lower infiltration of various immune cell types in metastases with the exception of CD4+ T cells and M2 macrophages. These results are in line with lower densities of immune cells and higher CD4/CD8 ratios in metastases shown by IHC. Analysis of transcriptomic data from autopsy cases and animal models confirms that immunosuppression is also present in extracranial metastases. Significantly higher somatic copy number aberration and allelic imbalance burdens are identified in metastases.

CONCLUSIONS

Metastasis is a molecularly late event, and immunosuppression driven by different molecular events, including somatic copy number aberration, may be a common characteristic of tumors with metastatic plasticity.

Neoadjuvant Chemotherapy Increases Cytotoxic T Cell, Tissue Resident Memory T Cell, and B Cell Infiltration in Resectable NSCLC

INTRODUCTION

The combination of programmed cell death protein-1 or programmed death-ligand 1 immune checkpoint blockade and chemotherapy has revolutionized the treatment of advanced NSCLC, but the mechanisms underlying this synergy remain incompletely understood. In this study, we explored the relationships between neoadjuvant chemotherapy and the immune microenvironment (IME) of resectable NSCLC to identify novel mechanisms by which chemotherapy may enhance the effect of immune checkpoint blockade.

METHODS

Genomic, transcriptomic, and immune profiling data of 511 patients treated with neoadjuvant chemotherapy followed by surgery (NCT) versus upfront surgery (US) were compared with determined differential characteristics of the IMEs derived from whole-exome sequencing (NCT = 18; US = 73), RNA microarray (NCT = 45; US = 202), flow cytometry (NCT = 17; US = 39), multiplex immunofluorescence (NCT = 10; US = 72), T-cell receptor sequencing (NCT = 16 and US = 63), and circulating cytokines (NCT = 18; US = 73).

RESULTS

NCT was associated with increased infiltration of cytotoxic CD8⁺ T cells and CD20⁺ B cells. Moreover, NCT was associated with increases in CD8⁺CD103⁺ and CD4⁺CD103⁺PD-1⁺TIM3^- tissue resident memory T cells. Gene expression profiling supported memory function of CD8⁺ and CD4⁺ T cells. However, NCT did not affect T-cell receptor clonality, richness, or tumor mutational burden. Finally, NCT was associated with decreased plasma BDNF (TrkB) at baseline and week 4 after surgery.

CONCLUSIONS

Our study supports that, in the context of resectable NSCLC, neoadjuvant chemotherapy promotes antitumor immunity through T and B cell recruitment in the IME and through a phenotypic change toward cytotoxic and memory CD8⁺ and CD4⁺ memory helper T cells.

The TRIXS end-station for femtosecond time-resolved resonant inelastic x-ray scattering experiments at the soft x-ray free-electron laser FLASH

We present the experimental end-station TRIXS dedicated to time-resolved soft x-ray resonant inelastic x-ray scattering (RIXS) experiments on solid samples at the free-electron laser FLASH. Using monochromatized ultrashort femtosecond XUV/soft x-ray photon pulses in combination with a synchronized optical laser in a pump-probe scheme, the TRIXS setup allows measuring sub-picosecond time-resolved high-resolution RIXS spectra in the energy range from 35 eV to 210 eV, thus spanning the M-edge (M₁ and M_2,3) absorption resonances of 3d transition metals and N_4,5-edges of rare earth elements. A Kirkpatrick-Baez refocusing mirror system at the first branch of the plane grating monochromator beamline (PG1) provides a focus of (6 × 6) μm² (FWHM) at the sample. The RIXS spectrometer reaches an energy resolution of 35-160 meV over the entire spectral range. The optical laser system based on a chirped pulse optical parametric amplifier provides approximately 100 fs (FWHM) long photon pulses at the fundamental wavelength of 800 nm and a fluence of 120 mJ/cm² at a sample for optical pump-XUV probe measurements. Furthermore, optical frequency conversion enables experiments at 400 nm or 267 nm with a fluence of 80 and 30 mJ/cm², respectively. Some of the first (pump-probe) RIXS spectra measured with this setup are shown. The measured time resolution for time-resolved RIXS measurements has been characterized as 287 fs (FWHM) for the used energy resolution.

Abstract 399: Tertiary lymphoid structures features associate with outcome in non-small cell lung carcinoma

Introduction: Tertiary lymphoid structures (TLS) are ectopic lymphoid structures organized in a nodular pattern secondary to chronic inflammation. The presence of TLS has been associated with efficacious response to immune checkpoint blockade in a wide array of tumor types, including Non-Small Cell Lung Carcinoma (NSCLC). At early stages of NSCLC, TLS are observed in up to 70% of primary tumors where they are also associated with effective anti-tumor immune-responses and response to anti-PD-1 therapy. Despite complete surgical resection, however, up to 50% of patients with early stage NSCLC will eventually relapse. A systematic histomorphologic analysis of TLS and their association with relapse has not been well characterized.

Design: Serial sections from Formalin-Fixed Paraffin-Embedded (FFPE) tissue specimens from 33 patients with stage I NSCLC were obtained (23 adenocarcinomas and 10 squamous cell carcinomas). Slides were stained with hematoxylin and eosin (H&E), immunohistochemistry for Vimentin, and 2 multiplex immunofluorescence (mIF) panels (Panel 1: cytokeratin (CK), CD3, CD8, CD68, PD-1, PD-L1 and DAPI Panel 2: CK, CD3 CD8, CD45RO, FOXP3, Granzyme B and DAPI). TLS were classified based on the H&E and Vimentin analysis into lymphoid aggregates (LA), immature TLS (iTLS), and mature TLS (mTLS). Morphometric analysis (number, area, and distance to nearest malignant cell) of intratumoral (IT) and peritumoral (PT) TLS was performed. mIF slides were scanned and IT representative areas were selected for cell densities and percentage quantification of immune-phenotypes. Morphometric analysis of TLSs was correlated with clinical outcomes, and tumor infiltrating immune cells.

Results: Patients who recurred had fewer IT mTLS (1.8 vs 6; p=0.02) and smaller area of mTLS (64541.7μm2 vs 149870.5 μm2; p=0.004) compared with patients who did not recurr.The cell density of IT antigen-experienced cytotoxic T-lymphocytes (CTLs), regulatory T-cells, memory CTLs, and memory CTLs expressing FOXP3 were inversely correlated with the mTLS area (r=-0.6, p≤0.02). Antigen-experienced CTLs (r=-0.61; p≤0.05) and non-CTLs (r=-0.63; p≤0.05) were inversely correlated with the number of IT mTLS.

Conclusion: Detailed morphometric analysis of mTLS offers relevant prognostic information for recurrence at stage I of NSCLC. mTLSs are associated with reduced IT infiltration by PD-1+ TILs. Taken together, germinal center development in mTLSs might convey a protective immune response due to immuno-dominant neoantigen presentation. Supported in part by CPRIT RP160668 grant and UT Lung SPORE.

Citation Format: Rossana Natalia Lazcano Segura, Andre Catao, Luisa M. Solis, Mei Jiang, Auriole Tamegnon, Junya Fujimoto, Jaime Rodriguez-Canales, Chi-Wan B. Chow, Carmen Behrens, Neda Kalhor, Annika Weissferdt, John Heymach, Stephen Swisher, Boris Sepesi, Jack Lee, Cesar Moran, Andrew Futreal, Jianjun Zhang, Edwin R. Parra, Ignacio I. Wistuba, Michael T. Tetzlaff, Alejandro Francisco-Cruz. Tertiary lymphoid structures features associate with outcome in non-small cell lung carcinoma [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 399.

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.

Computational Staining of Pathology Images to Study the Tumor Microenvironment in Lung Cancer

The spatial organization of different types of cells in tumor tissues reveals important information about the tumor microenvironment (TME). To facilitate the study of cellular spatial organization and interactions, we developed Histology-based Digital-Staining, a deep learning-based computation model, to segment the nuclei of tumor, stroma, lymphocyte, macrophage, karyorrhexis, and red blood cells from standard hematoxylin and eosin-stained pathology images in lung adenocarcinoma. Using this tool, we identified and classified cell nuclei and extracted 48 cell spatial organization-related features that characterize the TME. Using these features, we developed a prognostic model from the National Lung Screening Trial dataset, and independently validated the model in The Cancer Genome Atlas lung adenocarcinoma dataset, in which the predicted high-risk group showed significantly worse survival than the low-risk group (P = 0.001), with a HR of 2.23 (1.37-3.65) after adjusting for clinical variables. Furthermore, the image-derived TME features significantly correlated with the gene expression of biological pathways. For example, transcriptional activation of both the T-cell receptor and programmed cell death protein 1 pathways positively correlated with the density of detected lymphocytes in tumor tissues, while expression of the extracellular matrix organization pathway positively correlated with the density of stromal cells. In summary, we demonstrate that the spatial organization of different cell types is predictive of patient survival and associated with the gene expression of biological pathways. SIGNIFICANCE: These findings present a deep learning-based analysis tool to study the TME in pathology images and demonstrate that the cell spatial organization is predictive of patient survival and is associated with gene expression.See related commentary by Rodriguez-Antolin, p. 1912.

Genomic assessment distinguishes intrapulmonary metastases from synchronous primary lung cancers

Background

Multiple synchronous lung tumors (MSLT), particularly within a single lobe, represent a diagnostic and treatment challenge. While histologic assessment was once the only method to possibly distinguish multiple primary lung cancers, there is a growing interest in identifying unique genomic features or mutations to best characterize these processes.

Methods

In order to differentiate multiple primary lung malignancies from intrapulmonary metastases in patients with MSLT, we performed whole exome sequencing (WES) on 10 tumor samples from 4 patients with MSLT.

Results

Shared mutations between tumors from the same patient varied from 0-91%. Patient 3 shared no common mutations; however, in Patients 2 and 4, identical mutations were identified among all tumors from each patient, suggesting that the three tumors identified in Patient 3 represent separate primary lung cancers, while those of Patients 1, 2 and 4 signify hematogenous and lymphatic spread.

Conclusions

A high proportion of shared mutations between different lung tumors is likely indicative of intrapulmonary metastatic disease, while tumors with distinct genomic profiles likely represent multiple primary malignancies driven by distinct molecular events. Application of genomic profiling in the clinical setting may prove to be important to precise management of patients with MSLT.

eIF5B drives integrated stress response-dependent translation of PD-L1 in lung cancer

Cancer cells express high levels of PD-L1, a ligand of the PD-1 receptor on T cells, allowing tumors to suppress T cell activity. Clinical trials utilizing antibodies that disrupt the PD-1/PD-L1 checkpoint have yielded remarkable results, with anti-PD-1 immunotherapy approved as first-line therapy for lung cancer patients. We used CRISPR-based screening to identify regulators of PD-L1 in human lung cancer cells, revealing potent induction of PD-L1 upon disruption of heme biosynthesis. Impairment of heme production activates the integrated stress response (ISR), allowing bypass of inhibitory upstream open reading frames in the PD-L1 5' UTR, resulting in enhanced PD-L1 translation and suppression of anti-tumor immunity. We demonstrated that ISR-dependent PD-L1 translation requires the translation initiation factor eIF5B. eIF5B overexpression, which is frequent in lung adenocarcinomas and associated with poor prognosis, is sufficient to induce PD-L1. These findings illuminate mechanisms of immune checkpoint activation and identify targets for therapeutic intervention.

STING Pathway Expression Identifies NSCLC With an Immune-Responsive Phenotype

INTRODUCTION

Although the combination of anti-programmed cell death-1 or anti-programmed cell death ligand-1 (PD-L1) with platinum chemotherapy is a standard of care for NSCLC, clinical responses vary. Even though predictive biomarkers (which include PD-L1 expression, tumor mutational burden, and inflamed immune microenvironment) are validated for immunotherapy, their relevance to chemoimmunotherapy combinations is less clear. We have recently reported that activation of the stimulator of interferon genes (STING) innate immune pathway enhances immunotherapy response in SCLC. Here, we hypothesize that STING pathway activation may predict and underlie predictive correlates of antitumor immunity in NSCLC.

METHODS

We analyzed transcriptomic and proteomic profiles in two NSCLC cohorts from our institution (treatment-naive patients in the Profiling of Resistance Patterns and Oncogenic Signaling Pathways in Evaluation of Cancers of the Thorax study and relapsed patients in the Biomarker-Integrated Approaches of Targeted Therapy for Lung Cancer Elimination study) and The Cancer Genome Atlas (N = 1320). Tumors were stratified by STING activation on the basis of protein or mRNA expression of cyclic GMP-AMP synthase, phospho-STING, and STING-mediated chemokines (chemokine ligand 5 [CCL5] and C-X-C motif chemokine 10 [CXCL10]). STING activation in patient tumors and in platinum-treated preclinical NSCLC models was correlated with biomarkers of immunotherapy response.

RESULTS

STING activation is associated with higher levels of intrinsic DNA damage, targetable immune checkpoints, and chemokines in treatment-naive and relapsed lung adenocarcinoma. We observed that tumors with lower STING and immune gene expression show higher frequency of serine-threonine kinase 11 (STK11) mutations; however, we identified a subset of these tumors that are TP53 comutated and display high immune- and STING-related gene expression. Treatment with cisplatin increases STING pathway activation and PD-L1 expression in multiple NSCLC preclinical models, including adeno- and squamous cell carcinoma.

CONCLUSIONS

STING pathway activation in NSCLC predicts features of immunotherapy response and is enhanced by cisplatin treatment. This suggests a possible predictive biomarker and mechanism for improved response to chemoimmunotherapy combinations.

A Phase I/II Study of Neoadjuvant Cisplatin, Docetaxel, and Nintedanib for Resectable Non-Small Cell Lung Cancer

PURPOSE

Nintedanib enhances the activity of chemotherapy in metastatic non-small cell lung cancer (NSCLC). In this phase I/II study, we assessed safety and efficacy of nintedanib plus neoadjuvant chemotherapy, using major pathologic response (MPR) as primary endpoint.

PATIENTS AND METHODS

Eligible patients had stage IB (≥4 cm)-IIIA resectable NSCLC. A safety run-in phase was followed by an expansion phase with nintedanib 200 mg orally twice a day (28 days), followed by three cycles of cisplatin (75 mg/m²), docetaxel (75 mg/m²) every 21 days plus nintedanib, followed by surgery. With 33 planned patients, the study had 90% power to detect an MPR increase from 15% to 35%.

RESULTS

Twenty-one patients (stages I/II/III, N = 1/8/12) were treated. One of 15 patients treated with nintedanib 200 mg achieved MPR [7%, 95% confidence interval (CI), 0.2%-32%]. Best ORR in 20 evaluable patients was 30% (6/20, 95% CI, 12%-54%). Twelve-month recurrence-free survival and overall survival were 66% (95% CI, 47%-93%) and 91% (95% CI, 79%-100%), respectively. Most frequent treatment-related grade 3-4 toxicities were transaminitis and electrolyte abnormalities. On the basis of an interim analysis the study was discontinued for futility. Higher levels of CD3⁺ and cytotoxic CD3⁺CD8⁺ T cells were found in treated tumors of patients who were alive than in those who died (652.8 vs. 213.4 cells/mm², P = 0.048; 142.3 vs. 35.6 cells/mm², P = 0.018).

CONCLUSIONS

Although tolerated, neoadjuvant nintedanib plus chemotherapy did not increase MPR rate compared with chemotherapy historical controls. Additional studies of the combination in this setting are not recommended. Posttreatment levels of tumor-infiltrating T cells were associated with patient survival. Use of MPR facilitates the rapid evaluation of neoadjuvant therapies.See related commentary by Blakely and McCoach, p. 3499.