High-Resolution Profiling of Lung Adenocarcinoma Identifies Expression Subtypes with Specific Biomarkers and Clinically Relevant Vulnerabilities

RNA sequencing identifies lung cancer lineage and facilitates drug repositioning

Recent breakthrough therapies have improved survival rates in non-small cell lung cancer (NSCLC), but a paradigm for prospective confirmation is still lacking. Patientdatasets were mainly downloaded from TCGA, CPTAC and GEO. We conducted downstream analysis by collecting metagenes and generated 42-gene subtype classifiers to elucidate biological pathways. Subsequently, scRNA, eRNA, methylation, mutation, and copy number variation were depicted from a phenotype perspective. Enhancing the clinical translatability of molecular subtypes, preclinical models including CMAP, CCLE, and GDSC were utilized for drug repositioning. Importantly, we verified the presence of previously described three phenotypes including bronchioid, neuroendocrine, and squamoid. Poor prognosis was seen in squamoid and neuroendocrine clusters for treatment-naive and immunotherapy populations. The neuroendocrine cluster was dominated by STK11 mutations and 14q13.3 amplifications, whose related methylated loci are predictive of immunotherapy. And the greatest therapeutic potential lies in the bronchioid cluster. We further estimated the relative cell abundance of the tumor microenvironment (TME), specific cell types could be reflected among three clusters. Meanwhile, the higher portion of immune cell infiltration belonged to bronchioid and squamoid, not the neuroendocrine cluster. In drug repositioning, MEK inhibitors resisted bronchioid but were squamoid-sensitive. To conceptually validate compounds/targets, we employed RNA-seq and CCK-8/western blot assays. Our results indicated that dinaciclib and alvocidib exhibited similar activity and sensitivity in the neuroendocrine cluster. Also, a lineage factor named KLF5 recognized by inferred transcriptional factors activity could be suppressed by verteporfin.

Dissecting transcriptome signals of anti-PD-1 response in lung adenocarcinoma

Immune checkpoint blockades are actively adopted in diverse cancer types including metastatic melanoma and lung cancer. Despite of durable response in 20-30% of patients, we still lack molecular markers that could predict the patient responses reliably before treatment. Here we present a composite model for predicting anti-PD-1 response based on tumor mutation burden (TMB) and transcriptome sequencing data of 85 lung adenocarcinoma (LUAD) patients who received anti-PD-(L)1 treatment. We found that TMB was a good predictor (AUC = 0.81) for PD-L1 negative patients (n = 20). For PD-L1 positive patients (n = 65), we built an ensemble model of 100 XGBoost learning machines where gene expression, gene set activities and cell type composition were used as input features. The transcriptome-based models showed excellent accuracy (AUC > 0.9) and highlighted the contribution of T cell activities. Importantly, nonresponder patients with high prediction score turned out to have high CTLA4 expression, which suggested that neoadjuvant CTLA4 combination therapy might be effective for these patients. Our data and analysis results provide valuable insights into developing biomarkers and strategies for treating LUAD patients using immune checkpoint inhibitors.

Systematic dissection of tumor-normal single-cell ecosystems across a thousand tumors of 30 cancer types

The complexity of the tumor microenvironment poses significant challenges in cancer therapy. Here, to comprehensively investigate the tumor-normal ecosystems, we perform an integrative analysis of 4.9 million single-cell transcriptomes from 1070 tumor and 493 normal samples in combination with pan-cancer 137 spatial transcriptomics, 8887 TCGA, and 1261 checkpoint inhibitor-treated bulk tumors. We define a myriad of cell states constituting the tumor-normal ecosystems and also identify hallmark gene signatures across different cell types and organs. Our atlas characterizes distinctions between inflammatory fibroblasts marked by AKR1C1 or WNT5A in terms of cellular interactions and spatial co-localization patterns. Co-occurrence analysis reveals interferon-enriched community states including tertiary lymphoid structure (TLS) components, which exhibit differential rewiring between tumor, adjacent normal, and healthy normal tissues. The favorable response of interferon-enriched community states to immunotherapy is validated using immunotherapy-treated cancers (n = 1261) including our lung cancer cohort (n = 497). Deconvolution of spatial transcriptomes discriminates TLS-enriched from non-enriched cell types among immunotherapy-favorable components. Our systematic dissection of tumor-normal ecosystems provides a deeper understanding of inter- and intra-tumoral heterogeneity.

IBPGNET: lung adenocarcinoma recurrence prediction based on neural network interpretability

Lung adenocarcinoma (LUAD) is the most common histologic subtype of lung cancer. Early-stage patients have a 30-50% probability of metastatic recurrence after surgical treatment. Here, we propose a new computational framework, Interpretable Biological Pathway Graph Neural Networks (IBPGNET), based on pathway hierarchy relationships to predict LUAD recurrence and explore the internal regulatory mechanisms of LUAD. IBPGNET can integrate different omics data efficiently and provide global interpretability. In addition, our experimental results show that IBPGNET outperforms other classification methods in 5-fold cross-validation. IBPGNET identified PSMC1 and PSMD11 as genes associated with LUAD recurrence, and their expression levels were significantly higher in LUAD cells than in normal cells. The knockdown of PSMC1 and PSMD11 in LUAD cells increased their sensitivity to afatinib and decreased cell migration, invasion and proliferation. In addition, the cells showed significantly lower EGFR expression, indicating that PSMC1 and PSMD11 may mediate therapeutic sensitivity through EGFR expression.

Inflamed immune phenotype predicts favorable clinical outcomes of immune checkpoint inhibitor therapy across multiple cancer types

BACKGROUND

The inflamed immune phenotype (IIP), defined by enrichment of tumor-infiltrating lymphocytes (TILs) within intratumoral areas, is a promising tumor-agnostic biomarker of response to immune checkpoint inhibitor (ICI) therapy. However, it is challenging to define the IIP in an objective and reproducible manner during manual histopathologic examination. Here, we investigate artificial intelligence (AI)-based immune phenotypes capable of predicting ICI clinical outcomes in multiple solid tumor types.

METHODS

Lunit SCOPE IO is a deep learning model which determines the immune phenotype of the tumor microenvironment based on TIL analysis. We evaluated the correlation between the IIP and ICI treatment outcomes in terms of objective response rates (ORR), progression-free survival (PFS), and overall survival (OS) in a cohort of 1,806 ICI-treated patients representing over 27 solid tumor types retrospectively collected from multiple institutions.

RESULTS

We observed an overall IIP prevalence of 35.2% and significantly more favorable ORRs (26.3% vs 15.8%), PFS (median 5.3 vs 3.1 months, HR 0.68, 95% CI 0.61 to 0.76), and OS (median 25.3 vs 13.6 months, HR 0.66, 95% CI 0.57 to 0.75) after ICI therapy in IIP compared with non-IIP patients, respectively (p<0.001 for all comparisons). On subgroup analysis, the IIP was generally prognostic of favorable PFS across major patient subgroups, with the exception of the microsatellite unstable/mismatch repair deficient subgroup.

CONCLUSION

The AI-based IIP may represent a practical, affordable, clinically actionable, and tumor-agnostic biomarker prognostic of ICI therapy response across diverse tumor types.

Transcriptional analysis of landmark molecular pathways in lung adenocarcinoma results in a clinically relevant classification with potential therapeutic implications

Lung adenocarcinoma (LUAD) is a molecularly heterogeneous disease. In addition to genomic alterations, cancer transcriptional profiling can be helpful to tailor cancer treatment and to estimate each patient's outcome. Transcriptional activity levels of 50 molecular pathways were inferred in 4573 LUAD patients using Gene Set Variation Analysis (GSVA) method. Seven LUAD subtypes were defined and independently validated based on the combined behavior of the studied pathways: AD (adenocarcinoma subtype) 1-7. AD1, AD4, and AD5 subtypes were associated with better overall survival. AD1 and AD4 subtypes were enriched in epidermal growth factor receptor (EGFR) mutations, whereas AD2 and AD6 showed higher tumor protein p53 (TP53) alteration frequencies. AD2 and AD6 subtypes correlated with higher genome instability, proliferation-related pathway expression, and specific sensitivity to chemotherapy, based on data from LUAD cell lines. LUAD subtypes were able to predict immunotherapy response in addition to CD274 (PD-L1) gene expression and tumor mutational burden (TMB). AD2 and AD4 subtypes were associated with potential resistance and response to immunotherapy, respectively. Thus, analysis of transcriptomic data could improve patient stratification beyond genomics and single biomarkers (i.e., PD-L1 and TMB) and may lay the foundation for more personalized treatment avenues, especially in driver-negative LUAD.

High-resolution transcriptomics analysis of CXCL13+ EPSTI1+ CDK1+ cells with a specific focus on lung adenocarcinoma

Background

Programmed cell death ligand 1 (PD-L1) blocking therapy has transformed the treatment of lung adenocarcinoma (LUAD), which has significantly changed the landscape of immunotherapy. We aimed to explore specific cell subpopulations to understand tumor progression and identify markers of response to PD-L1 blocking therapy.

Methods

Bulk, fluorescence-activated cell sorting (FACS), and single-cell RNA (scRNA) sequencing were used to profile CXCL13, EPSTI1, and CDK1. The gene set variation analysis (GSVA) R package was utilized for score calculation, and prognostic analyses included receiver operating characteristic (ROC) curves, Cox proportional hazard models, and meta-analysis. Additionally, we analyzed tumor microenvironment (TME), genomics, compound perturbations, and clinical indicators. The high-dimensional analysis captured the intrinsic characteristics of the subpopulation. Furthermore, subpopulation differential genes were used for enrichment analysis of transcription factors and compounds.

Results

Literature and website analyses supported the essential role of CXCL13, CDK1, and EPSTI1 in immunotherapy. This led us to focus specifically on LUAD by representing a pan-cancer profile of immune-sensitive genes. Logically, the high-characteristic population may consist of samples positive for CXCL13, EPSTI1, and CDK1. The three-gene signature was a favorable indicator of immunotherapy response in the Stand Up to Cancer-Mark Foundation (SU2C-MARK) LUAD cohort but showed a poor prognosis before treatment in the Lung Cancer Explorer (LCE) database. Further mechanistic exploration revealed specific mutations associated with the three-gene signature in SU2C-MARK LUAD, such as STK11. In The Cancer Genome Atlas (TCGA)-LUAD cohort, the high-scoring group exhibited a higher tumor mutational burden (TMB) and global methylation but a lower fraction genome altered (FGA) and estimated tumor purity. Moreover, dasatinib demonstrated sensitivity in the high-scoring group. The co-localization of the CXCL13, EPSTI1, and CDK1 subpopulation was validated through spatial transcriptome and immunohistochemical databases. Assessment of the subpopulation depicted high-resolution intercellular communication. Maintenance of specific pathways, such as TNF, CD74, and CD44, contributed to immunotherapy sensitivity. Finally, the subpopulation-enriched targets and drugs were confirmed through ConnectivityMap (CMAP) analysis and multi-omics, respectively.

Conclusions

In this study, positive samples for CXCL13, EPSTI1, and CDK1 exhibited poor prognostic significance in treatment-naïve LUAD cases but demonstrated benefits from PD-L1 blockade and dasatinib therapies.

Evolutionary proteogenomic landscape from pre-invasive to invasive lung adenocarcinoma

Lung adenocarcinoma follows a stepwise progression from pre-invasive to invasive. However, there remains a knowledge gap regarding molecular events from pre-invasive to invasive. Here, we conduct a comprehensive proteogenomic analysis comprising whole-exon sequencing, RNA sequencing, and proteomic and phosphoproteomic profiling on 98 pre-invasive and 99 invasive lung adenocarcinomas. The deletion of chr4q12 contributes to the progression from pre-invasive to invasive adenocarcinoma by downregulating SPATA18, thus suppressing mitophagy and promoting cell invasion. Proteomics reveals diverse enriched pathways in normal lung tissues and pre-invasive and invasive adenocarcinoma. Proteomic analyses identify three proteomic subtypes, which represent different stages of tumor progression. We also illustrate the molecular characterization of four immune clusters, including endothelial cells, B cells, DCs, and immune depression subtype. In conclusion, this comprehensive proteogenomic study characterizes the molecular architecture and hallmarks from pre-invasive to invasive lung adenocarcinoma, guiding the way to a deeper understanding of the tumorigenesis and progression of this disease.

High-Resolution Profiling of Head and Neck Squamous Cells Carcinoma Identifies Specific Biomarkers and Expression Subtypes of Clinically Relevant Vulnerabilities

BACKGROUND

Head and neck squamous cell carcinoma (HNSC) is the seventh most common cancer worldwide. Although there are several options for the treatment of HNSC, there is still a lack of better biomarkers to accurately predict the response to treatment and thus be more able to correctly treat the therapeutic modality.

METHODS

First, we typed cases from the TCGA-HNSC cohort into subtypes by a Bayesian non-negative matrix factorization (BayesNMF)-based consensus clustering approach. Subsequently, genomic and proteomic data from HNSC cell lines were integrated to identify biomarkers of response to targeted therapies and immunotherapies. Finally, associations between HNSC subtypes and CD8 T-cell-associated effector molecules, common immune checkpoint genes, were compared to assess the potential of HNSC subtypes as clinically predictive immune checkpoint blockade therapy.

RESULTS

The 500 HNSC cases from TCGA were put through a consensus clustering approach to identify six HNSC expression subtypes. In addition, subtypes with unique proteomics and dependency profiles were defined based on HNSC cell line histology and proteomics data. Subtype 4 (S4) exhibits hyperproliferative and hyperimmune properties, and S4-associated cell lines show specific vulnerability to ADAT2, EIF5AL1, and PAK2. PD-L1 and CASP1 inhibitors have therapeutic potential in S4, and we have also demonstrated that S4 is more responsive to immune checkpoint blockade therapy.

CONCLUSION

Overall, our HNSC typing approach identified robust tumor-expressing subtypes, and data from multiple screens also revealed subtype-specific biology and vulnerabilities. These HNSC expression subtypes and their biomarkers will help develop more effective therapeutic strategies.

The circadian clock circuitry deconvolutes colorectal cancer and lung adenocarcinoma heterogeneity in a dynamic time-related framework

Increasing evidence imputes cancer progression and resistance to therapy to intra-tumor molecular heterogeneity set off by cancer cell plasticity. Re-activation of developmental programs strictly linked to epithelial-to-mesenchymal transition and gaining of stem cells properties are crucial in this setting. Many biological processes involved in cancer onset and progression show rhythmic fluctuations driven by the circadian clock circuitry. Novel cancer patient stratification tools taking into account the temporal dimension of these biological processes are definitely needed. Lung cancer and colorectal cancer (CRC) are the leading causes of cancer death worldwide. Here, by developing an innovative computational approach we named Phase-Finder, we show that the molecular heterogeneity characterizing the two deadliest cancers, CRC and lung adenocarcinoma (LUAD), rather than a merely stochastic event is the readout of specific cancer molecular states which correlate with time-qualified patterns of gene expression. We performed time-course transcriptome analysis of CRC and LUAD cell lines and upon computing circadian genes expression-based correlation matrices we derived pseudo-time points to infer time-qualified patterns in the transcriptomic analysis of real-world data (RWD) from large cohorts of CRC and LUAD patients. Our temporal classification of CRC and LUAD cohorts was able to effectively render time-specific patterns in cancer phenotype switching determining dynamical distribution of molecular subtypes impacting patient prognosis.

Integrative splicing-quantitative-trait-locus analysis reveals risk loci for non-small-cell lung cancer

Splicing quantitative trait loci (sQTLs) have been demonstrated to contribute to disease etiology by affecting alternative splicing. However, the role of sQTLs in the development of non-small-cell lung cancer (NSCLC) remains unknown. Thus, we performed a genome-wide sQTL study to identify genetic variants that affect alternative splicing in lung tissues from 116 individuals of Chinese ancestry, which resulted in the identification of 1,385 sQTL-harboring genes (sGenes) containing 378,210 significant variant-intron pairs. A comprehensive characterization of these sQTLs showed that they were enriched in actively transcribed regions, genetic regulatory elements, and splicing-factor-binding sites. Moreover, sQTLs were largely distinct from expression quantitative trait loci (eQTLs) and showed significant enrichment in potential risk loci of NSCLC. We also integrated sQTLs into NSCLC GWAS datasets (13,327 affected individuals and 13,328 control individuals) by using splice-transcriptome-wide association study (spTWAS) and identified alternative splicing events in 19 genes that were significantly associated with NSCLC risk. By using functional annotation and experiments, we confirmed an sQTL variant, rs35861926, that reduced the risk of lung adenocarcinoma (rs35861926-T, OR = 0.88, 95% confidence interval [CI]: 0.82-0.93, p = 1.87 × 10-5) by promoting FARP1 exon 20 skipping to downregulate the expression level of the long transcript FARP1-011. Transcript FARP1-011 promoted the migration and proliferation of lung adenocarcinoma cells. Overall, our study provided informative lung sQTL resources and insights into the molecular mechanisms linking sQTL variants to NSCLC risk.

Pan-cancer analysis of post-translational modifications reveals shared patterns of protein regulation

Post-translational modifications (PTMs) play key roles in regulating cell signaling and physiology in both normal and cancer cells. Advances in mass spectrometry enable high-throughput, accurate, and sensitive measurement of PTM levels to better understand their role, prevalence, and crosstalk. Here, we analyze the largest collection of proteogenomics data from 1,110 patients with PTM profiles across 11 cancer types (10 from the National Cancer Institute's Clinical Proteomic Tumor Analysis Consortium [CPTAC]). Our study reveals pan-cancer patterns of changes in protein acetylation and phosphorylation involved in hallmark cancer processes. These patterns revealed subsets of tumors, from different cancer types, including those with dysregulated DNA repair driven by phosphorylation, altered metabolic regulation associated with immune response driven by acetylation, affected kinase specificity by crosstalk between acetylation and phosphorylation, and modified histone regulation. Overall, this resource highlights the rich biology governed by PTMs and exposes potential new therapeutic avenues.

Computed Tomography-based Prognostication in Lung Adenocarcinomas through Histopathological Feature Learning: A Retrospective Multicenter Study

Rationale: Modeling imaging surrogates for well-validated histopathological risk factors would enable prognostication in early-stage lung adenocarcinomas. Objectives: We aimed to develop and validate computed tomography (CT)-based deep learning (DL) models for the prognostication of early-stage lung adenocarcinomas through learning histopathological features and to investigate the models' reproducibility using retrospective, multicenter datasets. Methods: Two DL models were trained to predict visceral pleural invasion and lymphovascular invasion, respectively, using preoperative chest CT scans from 1,426 patients with stage I-IV lung adenocarcinomas. The averaged model output was defined as the composite score and evaluated for the prognostic discrimination and its added value to clinicopathological factors in temporal (n = 610) and external test sets (n = 681) of stage I lung adenocarcinomas. The study outcomes were freedom from recurrence (FFR) and overall survival (OS). Interscan and interreader reproducibility were analyzed in 31 patients with lung cancer who underwent same-day repeated CT scans. Results: For the temporal test set, the time-dependent area under the receiver operating characteristic curve was 0.76 (95% confidence interval [CI], 0.71-0.81) for 5-year FFR and 0.67 (95% CI, 0.59-0.75) for 5-year OS. For the external test set, the area under the curve was 0.69 (95% CI, 0.63-0.75) for 5-year OS. The discrimination performance remained stable in 10-year follow-up for both outcomes. The prognostic value of the composite score was independent of and complementary to the clinical factors (adjusted per-percent hazard ratio for FFR [temporal test], 1.04 [95% CI, 1.03-1.05; P < 0.001]; OS [temporal test], 1.03 [95% CI, 1.02-1.04; P < 0.001]; OS [external test], 1.03 [95% CI, 1.02-1.04; P < 0.001]). The likelihood ratio tests indicated added value of the composite score (all P < 0.05). The interscan and interreader reproducibility were excellent (Pearson's correlation coefficient, 0.98 for both). Conclusions: The CT-based composite score obtained from DL of histopathological features predicted survival in early-stage lung adenocarcinomas with high reproducibility.

The role of molecular subtypes and immune infiltration characteristics based on disulfidptosis-associated genes in lung adenocarcinoma

Lung adenocarcinoma (LUAD) is the most common type of lung cancer which accounts for about 40% of all lung cancers. Early detection, risk stratification and treatment are important for improving outcomes for LUAD. Recent studies have found that abnormal accumulation of cystine and other disulfide occurs in the cell under glucose starvation, which induces disulfide stress and increases the content of disulfide bond in actin cytoskeleton, resulting in cell death, which is defined as disulfidptosis. Because the study of disulfidptosis is in its infancy, its role in disease progression is still unclear. In this study, we detected the expression and mutation of disulfidptosis genes in LUAD using a public database. Clustering analysis based on disulfidptosis gene was performed and differential genes of disulfidptosis subtype were analyzed. 7 differential genes of disulfidptosis subtype were used to construct a prognostic risk model, and the causes of prognostic differences were investigated by immune-infiltration analysis, immune checkpoint analysis, and drug sensitivity analysis. qPCR was used to verify the expression of 7 key genes in lung cancer cell line (A549) and normal bronchial epithelial cell line (BEAS-2B). Since G6PD had the highest risk factor of lung cancer, we further verified the protein expression of G6PD in lung cancer cells by western blot, and confirmed through colony formation experiment that interference with G6PD was able to significantly inhibit the proliferation ability of lung cancer cells. Our results provide evidence for the role of disulfidptosis in LUAD and provide new ideas for individualized precision therapy of LUAD.