Genome-wide identification of differentially methylated promoters and enhancers associated with response to anti-PD-1 therapy in non-small cell lung cancer

Genome-wide multimediator analyses using the generalized Berk-Jones statistics with the composite test

MOTIVATION

Mediation analysis is performed to evaluate the effects of a hypothetical causal mechanism that marks the progression from an exposure, through mediators, to an outcome. In the age of high-throughput technologies, it has become routine to assess numerous potential mechanisms at the genome or proteome scales. Alongside this, the necessity to address issues related to multiple testing has also arisen. In a sparse scenario where only a few genes or proteins are causally involved, conventional methods for assessing mediation effects lose statistical power because the composite null distribution behind this experiment cannot be attained. The power loss hence decreases the true mechanisms identified after multiple testing corrections. To fairly delineate a uniform distribution under the composite null, Huang (Genome-wide analyses of sparse mediation effects under composite null hypotheses. Ann Appl Stat 2019a;13:60-84; AoAS) proposed the composite test to provide adjusted P-values for single-mediator analyses.

RESULTS

Our contribution is to extend the method to multimediator analyses, which are commonly encountered in genomic studies and also flexible to various biological interests. Using the generalized Berk-Jones statistics with the composite test, we proposed a multivariate approach that favors dense and diverse mediation effects, a decorrelation approach that favors sparse and consistent effects, and a hybrid approach that captures the edges of both approaches. Our analysis suite has been implemented as an R package MACtest. The utility is demonstrated by analyzing the lung adenocarcinoma datasets from The Cancer Genome Atlas and Clinical Proteomic Tumor Analysis Consortium. We further investigate the genes and networks whose expression may be regulated by smoking-induced epigenetic aberrations.

AVAILABILITY AND IMPLEMENTATION

An R package MACtest is available on https://github.com/roqe/MACtest.

Potential biomarkers for immunotherapy in non-small-cell lung cancer

For individuals with advanced or metastatic non-small cell lung cancer (NSCLC), the primary treatment is platinum-based doublet chemotherapy. Immune checkpoint inhibitors (ICIs), primarily PD-1/PD-L1 and CTLA-4, have been found to be effective in patients with NSCLC who have no EGFR/ALK mutations. Furthermore, ICIs are considered a standard therapy. The quantity of fresh immunogenic antigens discovered by cytotoxic T cells was measured by PD-L1 expression and tumor mutational burden (TMB), which were the first biomarkers assessed in clinical trials. However, immunotherapy did not have response efficacy markers similar to targeted therapy, highlighting the significance of newly developed biomarkers. This investigation aims to review the research on immunotherapy for NSCLC, focusing primarily on the impact of biomarkers on efficacy prediction to determine whether biomarkers may be utilized to evaluate the effectiveness of immunotherapy.

A tumor-associated endothelial signature score model in immunotherapy and prognosis across pan-cancers

Background: The tumor-associated endothelial cell (TAE) component plays a vital role in tumor immunity. However, systematic tumor-associated endothelial-related gene assessment models for predicting cancer immunotherapy (CIT) responses and survival across human cancers have not been explored. Herein, we investigated a TAE gene risk model to predict CIT responses and patient survival in a pan-cancer analysis. Methods: We analyzed publicly available datasets of tumor samples with gene expression and clinical information, including gastric cancer, metastatic urothelial cancer, metastatic melanoma, non-small cell lung cancer, primary bladder cancer, and renal cell carcinoma. We further established a binary classification model to predict CIT responses using the least absolute shrinkage and selection operator (LASSO) computational algorithm. Results: The model demonstrated a high predictive accuracy in both training and validation cohorts. The response rate of the high score group to immunotherapy in the training cohort was significantly higher than that of the low score group, with CIT response rates of 51% and 27%, respectively. The survival analysis showed that the prognosis of the high score group was significantly better than that of the low score group (all p < 0·001). Tumor-associated endothelial gene signature scores positively correlated with immune checkpoint genes, suggesting that immune checkpoint inhibitors may benefit patients in the high score group. The analysis of TAE scores across 33 human cancers revealed that the TAE model could reflect immune cell infiltration and predict the survival of cancer patients. Conclusion: The TAE signature model could represent a CIT response prediction model with a prognostic value in multiple cancer types.

Identification of immunotherapy biomarkers for improving the clinical outcome of homologous recombination deficiency patients with lung adenocarcinoma

Homologous recombination deficiency (HRD) is a common molecular signature of genomic instability and has been shown to be a biomarker for targeted therapies. However, there is a lack of studies on the role of HRD changes in lung adenocarcinoma (LUAD) transcriptomics. HRD scores were determined using single nucleotide polymorphism (SNP) array data from LUAD patients from The Cancer Genome Atlas (TCGA) database. Transcriptional data from patients with different scores were analyzed to identify biomarkers associated with HRD. Candidate biomarkers were validated using Gene Expression Omnibus (GEO)-sourced datasets and an immunotherapy cohort. According to the bulk transcriptome and clinical characteristics of 912 LUAD patients and Single-cell RNA-seq of 9 LUAD patients from TCGA and GEO databases, we observed increased MS4A6A expression in HRD tumors; high MS4A6A expression predicted improved survival outcomes. Furthermore, a comprehensive analysis of the tumor immune microenvironment (TIME) revealed a positive correlation between MS4A6A expression and neoantigen loading and immune cell infiltration. Additionally, the immunotherapy cohort confirmed the possibility of using MS4A6A as a biomarker. Collectively, we suggest that MS4A6A is associated with HRD and provide a new perspective toward identifying promising biomarkers for immunotherapy.

Immune-related interaction perturbation networks unravel biological peculiars and clinical significance of glioblastoma

The immune system is an interacting network of plentiful molecules that could better characterize the relationship between immunity and cancer. This study aims to investigate the behavioral patterns of immune-related interaction perturbation networks in glioblastoma. An immune-related interaction-perturbation framework was introduced to characterize four heterogeneous subtypes using RNA-seq data of TCGA/CGGA glioblastoma tissues and GTEx normal brain tissues. The stability and robustness of the four subtypes were validated in public datasets and our in-house cohort. In the four subtypes, C1 was an inflammatory subtype with high immune infiltration, low tumor purity, and potential response to immunotherapy; C2, an invasive subtype, was featured with dismal prognosis, telomerase reverse transcriptase promoter mutations, moderate levels of immunity, and stromal constituents, as well as sensitivity to receptor tyrosine kinase signaling inhibitors; C3 was a proliferative subtype with high tumor purity, immune-desert microenvironment, sensitivity to phosphatidylinositol 3'-kinase signaling inhibitor and DNA replication inhibitors, and potential resistance to immunotherapy; C4, a synaptogenesis subtype with the best prognosis, exhibited high synaptogenesis-related gene expression, prevalent isocitrate dehydrogenase mutations, and potential sensitivity to radiotherapy and chemotherapy. Overall, this study provided an attractive platform from the perspective of immune-related interaction perturbation networks, which might advance the tailored management of glioblastoma.

Identification of the relationship between single-cell N6-methyladenosine regulators and the infiltrating immune cells in esophageal carcinoma

Background

N6-methyladenosine (m⁶A) RNA methylation plays a crucial role in important genomic processes in a variety of malignancies. However, the characterization of m⁶A with infiltrating immune cells in the tumor microenvironment (TME) in esophageal squamous carcinoma (ESCC) remains unknown.

Methods

The single-cell transcriptome data from five ESCC patients in our hospital were analyzed, and TME clusters associated with prognosis and immune checkpoint genes were investigated. Cell isolation and qPCR were conducted to validate the gene characterization in different cells.

Results

According to distinct biological processes and marker genes, macrophages, T cells, and B cells clustered into three to four different subgroups. In addition, we demonstrated that m⁶A RNA methylation regulators were strongly related to the clinical and biological features of ESCC. Analysis of transcriptome data revealed that m⁶A-mediated TME cell subsets had high predictive value and showed a close relationship with immune checkpoint genes. The validation results from qPCR demonstrated the characteristics of essential genes. CellChat analysis revealed that RNA from TME cells m⁶A methylation-associated cell subtypes had substantial and diversified interactions with cancer cells. Further investigation revealed that MIF- (CD74+CXCR4) and MIF- (CD74⁺CD44) ligand-receptor pairings facilitated communication between m⁶A-associated subtypes of TME cells and cancer cells.

Conclusion

Overall, our study demonstrated for the first time the function of m⁶A methylation-mediated intercellular communication in the microenvironment of tumors in controlling tumor development and anti-tumor immune regulation in ESCC.

Aggrephagy-related patterns in tumor microenvironment, prognosis, and immunotherapy for acute myeloid leukemia: a comprehensive single-cell RNA sequencing analysis

Acute myeloid leukemia (AML) is a complex mixed entity composed of malignant tumor cells, immune cells and stromal cells, with intra-tumor and inter-tumor heterogeneity. Single-cell RNA sequencing enables a comprehensive study of the highly complex tumor microenvironment, which is conducive to exploring the evolutionary trajectory of tumor cells. Herein, we carried out comprehensive analyses of aggrephagy-related cell clusters based on single-cell sequencing for patients with acute myeloid leukemia. A total of 11 specific cell types (T, NK, CMP, Myeloid, GMP, MEP, Promono, Plasma, HSC, B, and Erythroid cells) using t-SNE dimension reduction analysis. Several aggrephagy-related genes were highly expressed in the 11 specific cell types. Using Monocle analysis and NMF clustering analysis, six aggrephagy-related CD8+ T clusters, six aggrephagy-related NK clusters, and six aggrephagy-related Mac clusters were identified. We also evaluated the ligand-receptor links and Cell-cell communication using CellChat package and CellChatDB database. Furthermore, the transcription factors (TFs) of aggrephagy-mediated cell clusters for AML were assessed through pySCENIC package. Prognostic analysis of the aggrephagy-related cell clusters based on R package revealed the differences in prognosis of aggrephagy-mediated cell clusters. Immunotherapy of the aggrephagy-related cell clusters was investigated using TIDE algorithm and public immunotherapy cohorts. Our study revealed the significance of aggrephagy-related patterns in tumor microenvironment, prognosis, and immunotherapy for AML.

A risk score combining co-expression modules related to myeloid cells and alternative splicing associates with response to PD-1/PD-L1 blockade in non-small cell lung cancer

Background

Comprehensive analysis of transcriptomic profiles of non-small cell lung cancer (NSCLC) may provide novel evidence for biomarkers associated with response to PD-1/PD-L1 immune checkpoint blockade (ICB).

Methods

We utilized weighted gene co-expression network analysis (WGCNA) to analyze transcriptomic data from two NSCLC datasets from Gene Expression Omnibus (GSE135222 and GSE126044) that involved patients received ICB treatment. We evaluated the correlation of co-expression modules with ICB responsiveness and functionally annotated ICB-related modules using pathway enrichment analysis, single-cell RNA sequencing, flow cytometry and alternative splicing analysis. We built a risk score using Lasso-COX regression based on hub genes from ICB-related modules. We investigated the alteration of tumor microenvironment between high- and low- risk groups and the association of the risk score with previously established predictive biomarkers.

Results

Our results identified a black with positive correlation and a blue module with negative correlation to ICB responsiveness. The black module was enriched in pathway of T cell activation and antigen processing and presentation, and the genes assigned to it were consistently expressed on myeloid cells. We observed decreased alternative splicing events in samples with high signature scores of the blue module. The Lasso-COX analysis screened out three genes (EVI2B, DHX9, HNRNPM) and constructed a risk score from the hub genes of the two modules. We validated the predictive value of the risk score for poor response to ICB therapy in an in-house NSCLC cohort and a pan-cancer cohort from the KM-plotter database. The low-risk group had more immune-infiltrated microenvironment, with higher frequencies of precursor exhausted CD8⁺ T cells, tissue-resident CD8⁺ T cells, plasmacytoid dendritic cells and type 1 conventional dendritic cells, and a lower frequency of terminal exhausted CD8⁺ T cells, which may explain its superior response to ICB therapy. The significant correlation of the risk score to gene signature of tertiary lymphoid structure also implicated the possible mechanism of this predictive biomarker.

Conclusions

Our study identified two co-expression modules related to ICB responsiveness in NSCLC and developed a risk score accordingly, which could potentially serve as a predictive biomarker for ICB response.

Targeting immunosuppressive Ly6C+ classical monocytes reverses anti-PD-1/CTLA-4 immunotherapy resistance

Introduction

Despite significant clinical advancement with the use of immune checkpoint blockade (ICB) in non-small cell lung cancer (NSCLC) there are still a major subset of patients that develop adaptive/acquired resistance. Understanding resistance mechanisms to ICB is critical to developing new therapeutic strategies and improving patient survival. The dynamic nature of the tumor microenvironment and the mutational load driving tumor immunogenicity limit the efficacy to ICB. Recent studies indicate that myeloid cells are drivers of ICB resistance. In this study we sought to understand which immune cells were contributing to resistance and if we could modify them in a way to improve response to ICB therapy.

Results

Our results show that combination anti-PD-1/CTLA-4 produces an initial antitumor effect with evidence of an activated immune response. Upon extended treatment with anti-PD-1/CTLA-4 acquired resistance developed with an increase of the immunosuppressive populations, including T-regulatory cells, neutrophils and monocytes. Addition of anti-Ly6C blocking antibody to anti-PD-1/CTLA-4 was capable of completely reversing treatment resistance and restoring CD8 T cell activity in multiple KP lung cancer models and in the autochthonous lung cancer KrasLSL-G12D/p53fl/fl model. We found that there were higher classical Ly6C+ monocytes in anti-PD-1/CTLA-4 combination resistant tumors. B7 blockade illustrated the importance of dendritic cells for treatment efficacy of anti-Ly6C/PD-1/CTLA-4. We further determined that classical Ly6C+ monocytes in anti-PD-1/CTLA-4 resistant tumors are trafficked into the tumor via IFN-γ and the CCL2-CCR2 axis. Mechanistically we found that classical monocytes from ICB resistant tumors were unable to differentiate into antigen presenting cells and instead differentiated into immunosuppressive M2 macrophages or myeloid-derived suppressor cells (MDSC). Classical Ly6C+ monocytes from ICB resistant tumors had a decrease in both Flt3 and PU.1 expression that prevented differentiation into dendritic cells/macrophages.

Conclusions

Therapeutically we found that addition of anti-Ly6C to the combination of anti-PD-1/CTLA-4 was capable of complete tumor eradication. Classical Ly6C+ monocytes differentiate into immunosuppressive cells, while blockade of classical monocytes drives dendritic cell differentiation/maturation to reinvigorate the anti-tumor T cell response. These findings support that immunotherapy resistance is associated with infiltrating monocytes and that controlling the differentiation process of monocytes can enhance the therapeutic potential of ICB.

Mendelian randomization analysis identified tumor necrosis factor as being associated with severe COVID-19

Background: Observational studies have shown that anti-tumor necrosis factor (TNF) therapy may be beneficial for patients with coronavirus disease 2019 (COVID-19). Nevertheless, because of the methodological restrictions of traditional observational studies, it is a challenge to make causal inferences. This study involved a two-sample Mendelian randomization analysis to investigate the causal link between nine TNFs and COVID-19 severity using publicly released genome-wide association study summary statistics. Methods: Summary statistics for nine TNFs (21,758 cases) were obtained from a large-scale genome-wide association study. Correlation data between single-nucleotide polymorphisms and severe COVID-19 (18,152 cases vs. 1,145,546 controls) were collected from the COVID-19 host genetics initiative. The causal estimate was calculated by inverse variance-weighted (IVW), MR-Egger, and weighted median methods. Sensitivity tests were conducted to assess the validity of the causal relationship. Results: Genetically predicted TNF receptor superfamily member 6 (FAS) positively correlated with the severity of COVID-19 (IVW, odds ratio = 1.10, 95% confidence interval = 1.01-1.19, p = 0.026), whereas TNF receptor superfamily member 5 (CD40) was protective against severe COVID-19 (IVW, odds ratio = 0.92, 95% confidence interval = 0.87-0.97, p = 0.002). Conclusion: Genetic evidence from this study supports that the increased expression of FAS is associated with the risk of severe COVID-19 and that CD40 may have a potential protective effect against COVID-19.

A prognostic and immunotherapeutic predictive model based on the cell-originated characterization of tumor microenvironment in lung adenocarcinoma

Tumor microenvironment (TME) plays a crucial role in predicting prognosis and response to therapy in lung cancer. Our study established a prognostic and immunotherapeutic predictive model, the tumor immune cell score (TICS), by differentiating cell origins in lung adenocarcinoma (LUAD) based on the transcriptomic data of 2,510 patients in 14 independent cohorts, including 12 public datasets and two in-house cohorts. The high TICS was associated with prolonged overall survival (OS), especially in the early-stage LUAD. For the advanced-stage LUAD, high TICS predicted a superior OS in patients who were treated with immunotherapy instead of chemotherapy or TKI. The result suggested that TICS could serve as an indicator for the prognostic stratification management of patients in the early-stage LUAD, and as a potential guide for therapeutic decision-marking in the advanced-stage LUAD. Our findings provided an insight into prognosis stratification and potential guidance for treatment strategy selection.

Depiction of neuroendocrine features associated with immunotherapy response using a novel one-class predictor in lung adenocarcinoma

BACKGROUND

Tumours with no evidence of neuroendocrine transformation histologically but harbouring neuroendocrine features are collectively referred to as non-small cell lung cancer (NSCLC) with neuroendocrine differentiation (NED). Investigating the mechanisms underlying NED is conducive to designing appropriate treatment options for NSCLC patients.

METHODS

In the present study, we integrated multiple lung cancer datasets to identify neuroendocrine features using a one-class logistic regression (OCLR) machine learning algorithm trained on small cell lung cancer (SCLC) cells, a pulmonary neuroendocrine cell type, based on the transcriptome of NSCLC and named the NED index (NEDI). Single-sample gene set enrichment analysis, pathway enrichment analysis, ESTIMATE algorithm analysis, and unsupervised subclass mapping (SubMap) were performed to assess the altered pathways and immune characteristics of lung cancer samples with different NEDI values.

RESULTS

We developed and validated a novel one-class predictor based on the expression values of 13,279 mRNAs to quantitatively evaluate neuroendocrine features in NSCLC. We observed that a higher NEDI correlated with better prognosis in patients with LUAD. In addition, we observed that a higher NEDI was significantly associated with reduced immune cell infiltration and immune effector molecule expression. Furthermore, we found that etoposide-based chemotherapy might be more effective in the treatment of LUAD with high NEDI values. Moreover, we noted that tumours with low NEDI values had better responses to immunotherapy than those with high NEDI values.

CONCLUSIONS

Our findings improve the understanding of NED and provide a useful strategy for applying NEDI-based risk stratification to guide decision-making in the treatment of LUAD.

A novel methionine metabolism-related signature predicts prognosis and immunotherapy response in lung adenocarcinoma

Recent research revealed methionine metabolism as a key mediator of tumor initiation and immune evasion. However, the relationship between methionine metabolism and tumor microenvironment (TME) in lung adenocarcinoma (LUAD) remains unknown. Here, we comprehensively analyzed the genomic alterations, expression patterns, and prognostic values of 68 methionine-related regulators (MRGs) in LUAD. We found that most MRGs were highly prognostic based on 30 datasets including 5024 LUAD patients. Three distinct MRG modification patterns were identified, which showed significant differences in clinical outcomes and TME characteristics: The C2 subtype was characterized by higher immune score, while the C3 subtype had more malignant cells and worse survival. We developed a MethScore to measure the level of methionine metabolism in LUAD. MethScore was positively correlated with T-cell dysfunction and tumor-associated macrophages (TAMs), indicating a dysfunctional TME phenotype in the high MethScore group. In addition, two immunotherapy cohorts confirmed that patients with a lower MethScore exhibited significant clinical benefits. Our study highlights the important role of methionine metabolism in modeling the TME. Evaluating methionine modification patterns will enhance our understanding of TME characteristics and can guide more effective immunotherapy strategies.

Reshaping the tumour immune microenvironment in solid tumours via tumour cell and immune cell DNA methylation: from mechanisms to therapeutics

In recent years, the tumour microenvironment (TME) of solid tumours has attracted more and more attention from researchers, especially those non-tumour components such as immune cells. Infiltration of various immune cells causes tumour immune microenvironment (TIME) heterogeneity, and results in different therapeutic effects. Accumulating evidence showed that DNA methylation plays a crucial role in remodelling TIME and is associated with the response towards immune checkpoint inhibitors (ICIs). During carcinogenesis, DNA methylation profoundly changes, specifically, there is a global loss of DNA methylation and increased DNA methylation at the promoters of suppressor genes. Immune cell differentiation is disturbed, and exclusion of immune cells from the TME occurs at least in part due to DNA methylation reprogramming. Therefore, pharmaceutical interventions targeting DNA methylation are promising. DNA methyltransferase inhibitors (DNMTis) enhance antitumor immunity by inducing transcription of transposable elements and consequent viral mimicry. DNMTis upregulate the expression of tumour antigens, mediate immune cells recruitment and reactivate exhausted immune cells. In preclinical studies, DNMTis have shown synergistic effect when combined with immunotherapies, suggesting new strategies to treat refractory solid tumours.

Identification of a single cell-based signature for predicting prognosis risk and immunotherapy response in patients with glioblastoma

This study constructed a novel gene pair signature based on bulk and single-cell sequencing samples in relative expression order within the samples. The subsequent analysis included glioma samples from Xiangya Hospital. Gene pair signatures possessed a solid ability to predict the prognosis of glioblastoma and pan-cancer. Samples having different malignant biological hallmarks were distinguished by the algorithm, with the high gene pair score group featuring classic copy number variations, oncogenic mutations, and extensive hypomethylation, mediating poor prognosis. The increased gene pair score group with a poorer prognosis demonstrated significant enrichment in tumor and immune-related signaling pathways while presenting immunological diversity. The remarkable infiltration of M2 macrophages in the high gene pair score group was validated by multiplex immunofluorescence, suggesting that combination therapies targeting adaptive and innate immunity may serve as a therapeutic option. Overall, a gene pair signature applicable to predict prognosis hopefully provides a reference to guide clinical practice.

Deciphering the immune heterogeneity dominated by natural killer cells with prognostic and therapeutic implications in hepatocellular carcinoma

Belonging to type 1 innate lymphoid cells (ILC1), natural killer (NK) cells play an important role not only in fighting microbial infections but also in anti-tumor response. Hepatocellular carcinoma (HCC) represents an inflammation-related malignancy and NK cells are enriched in the liver, making them an essential component of the HCC immune microenvironment. In this study, we performed single-cell RNA-sequencing (scRNA-seq) analysis to identify the NK cell marker genes (NKGs) and uncovered 80 prognosis-related ones by the TCGA-LIHC dataset. Based on prognostic NKGs, HCC patients were categorized into two subtypes with distinct clinical outcomes. Subsequently, we conducted LASSO-COX and stepwise regression analysis on prognostic NKGs to establish a five-gene (UBB, CIRBP, GZMH, NUDC, and NCL) prognostic signature-NKscore. Different mutation statuses of the two risk groups stratified by NKscore were comprehensively characterized. Besides, the established NKscore-integrated nomogram presented enhanced predictive performance. Single sample gene set enrichment analysis (ssGSEA) analysis was used to uncover the landscape of the tumor immune microenvironment (TIME) and the high-NKscore risk group was characterized with an immune-exhausted phenotype while the low-NKscore risk group held relatively strong anti-cancer immunity. T cell receptor (TCR) repertoire, tumor inflammation signature (TIS), and Immunophenoscore (IPS) analyses revealed differences in immunotherapy sensitivity between the two NKscore risk groups. Taken together, we developed a novel NK cell-related signature to predict the prognosis and immunotherapy efficacy for HCC patients.

B cell-dependent subtypes and treatment-based immune correlates to survival in stage 3 and 4 lung adenocarcinomas

Lung cancer is the leading cause of cancer-related deaths worldwide. Surgery and chemoradiation are the standard of care in early stages of non-small cell lung cancer (NSCLC), while immunotherapy is the standard of care in late-stage NSCLC. The immune composition of the tumor microenvironment (TME) is recognized as an indicator for responsiveness to immunotherapy, although much remains unknown about its role in responsiveness to surgery or chemoradiation. In this pilot study, we characterized the NSCLC TME using mass cytometry (CyTOF) and bulk RNA sequencing (RNA-Seq) with deconvolution of RNA-Seq being performed by Kassandra, a recently published deconvolution tool. Stratification of patients based on the intratumoral abundance of B cells identified that the B-cell rich patient group had increased expression of CXCL13 and greater abundance of PD1+ CD8 T cells. The presence of B cells and PD1+ CD8 T cells correlated positively with the presence of intratumoral tertiary lymphoid structures (TLS). We then assessed the predictive and prognostic utility of these cell types and TLS within publicly available stage 3 and 4 lung adenocarcinoma (LUAD) RNA-Seq datasets. As previously described by others, pre-treatment expression of intratumoral 12-chemokine TLS gene signature is associated with progression free survival (PFS) in patients who receive treatment with immune checkpoint inhibitors (ICI). Notably and unexpectedly pre-treatment percentages of intratumoral B cells are associated with PFS in patients who receive surgery, chemotherapy, or radiation. Further studies to confirm these findings would allow for more effective patient selection for both ICI and non-ICI treatments.

An immune indicator based on BTK and DPEP2 identifies hot and cold tumors and clinical treatment outcomes in lung adenocarcinoma

In lung adenocarcinoma (LUAD), immune heterogeneity of hot and cold tumors has been recognized as one of the major factors affecting immunotherapy and other common treatments. However, there is still a lack of biomarkers that can effectively identify the immunophenotype of cold and hot tumors. First, the immune signatures were obtained based on literature mining, including macrophage/monocyte, IFN-γ response, TGF-β response, IL12 response, lymphocyte activation, and ECM/Dve/immune response. Subsequently, LUAD patients were further clustered into different immune phenotypes based on these immune signatures. Next, the key genes related to the immune phenotypes were screened by WGCNA analysis, univariate analysis, and lasso-cox analysis, and the risk signature was established via the key genes. In additional, we compared the clinicopathological characteristics, drug sensitivity, the abundance of immune infiltration, and the efficacy of immunotherapy and commonly used therapies between patients in the high- and low-risk groups in LUAD. LUAD patients were divided into immune hot phenotype and immune cold phenotype groups. The clinical presentation showed that patients with the immune hot phenotype had higher immunoactivity (including higher MHC, CYT, immune, stromal, ESTIMATE scores, higher abundance of immune cell infiltration, higher abundance of TIL, and enrichment of immune-enriched subtypes) and better survival outcomes than those with the immune cold phenotype. Subsequently, WGCNA analysis, univariate analysis, and lasso-cox analysis identified the genes highly associated with the immune phenotype: BTK and DPEP2. The risk signature, consisting of BTK and DPEP2, is highly correlated with the immune phenotype. High-risk scores were enriched in patients with immune cold phenotype and low-risk scores were enriched in patients with immune hot phenotype. Compared to the high-risk group, the low-risk group had better clinical performance, higher drug sensitivity, and a higher degree of immunoactivity, as well as better efficacy in receiving immunotherapy and common adjuvant therapy. This study developed an immune indicator consisting of BTK and DPEP2 based on the heterogeneity of hot and cold Immunophenotypes of the tumor microenvironment. This indicator has good efficacy in predicting prognosis and assessing the efficacy of immunotherapy, chemotherapy, and radiotherapy. It has the potential to facilitate personalized and precise treatment of LUAD in the future.

Integrative analysis revealed that distinct cuprotosis patterns reshaped tumor microenvironment and responses to immunotherapy of colorectal cancer

Background

Cuprotosis is a novel form of programmed cell death that involves direct targeting of key enzymes in the tricarboxylic acid (TCA) cycle by excess copper and may result in mitochondrial metabolic dysfunction. However, whether cuprotosis may mediate the tumor microenvironment (TME) and immune regulation in colorectal cancer (CRC) remains unclear.

Methods

Ten cuprotosis-related genes were selected and unsupervised consensus clustering was performed to identify the cuprotosis patterns and the correlated TME characteristics. Using principal component analysis, a COPsig score was established to quantify cuprotosis patterns in individual patients. The top 9 most important cuprotosis signature genes were analyzed using single-cell transcriptome data.

Results

Three distinct cuprotosis patterns were identified. The TME cell infiltration characteristics of three patterns were associated with immune-excluded, immune-desert, and immune-inflamed phenotype, respectively. Based on individual cuprotosis patterns, patients were assigned into high and low COPsig score groups. Patients with a higher COPsig score were characterized by longer overall survival time, lower immune cell as well as stromal infiltration, and greater tumor mutational burden. Moreover, further analysis demonstrated that CRC patients with a higher COPsig score were more likely to respond to immune checkpoint inhibitors and 5-fluorouracil chemotherapy. Single-cell transcriptome analysis indicated that cuprotosis signature genes recruited tumor-associated macrophages to TME through the regulation of TCA and the metabolism of glutamine and fatty acid, thus influencing the prognosis of CRC patients.

Conclusion

This study indicated that distinct cuprotosis patterns laid a solid foundation to the explanation of heterogeneity and complexity of individual TME, thus guiding more effective immunotherapy as well as adjuvant chemotherapy strategies.

CD83 expression characterizes precursor exhausted T cell population

T cell exhaustion is a main obstacle against effective cancer immunotherapy. Exhausted T cells include a subpopulation that maintains proliferative capacity, referred to as precursor exhausted T cells (T_PEX). While functionally distinct and important for antitumor immunity, T_PEX possess some overlapping phenotypic features with the other T-cell subsets within the heterogeneous tumor-infiltrating T-lymphocytes (TIL). Here we explore surface marker profiles unique to T_PEX using the tumor models treated by chimeric antigen receptor (CAR)-engineered T cells. We find that CD83 is predominantly expressed in the CCR7⁺PD1⁺ intratumoral CAR-T cells compared with the CCR7^-PD1⁺ (terminally differentiated) and CAR-negative (bystander) T cells. The CD83⁺CCR7⁺ CAR-T cells exhibit superior antigen-induced proliferation and IL-2 production compared with the CD83^- T cells. Moreover, we confirm selective expression of CD83 in the CCR7⁺PD1⁺ T-cell population in primary TIL samples. Our findings identify CD83 as a marker to discriminate T_PEX from terminally exhausted and bystander TIL.

Homologous recombination deficiency in triple-negative breast cancer: Multi-scale transcriptomics reveals distinct tumor microenvironments and limitations in predicting immunotherapy response

BACKGROUND

Triple-negative breast cancer (TNBC) is the most aggressive subtype of breast cancer and has the highest proportion of homologous recombination deficiency (HRD). HRD has been considered a biomarker of response to immune checkpoint inhibitors (ICIs), but the reality is more complicated. A comprehensive comparison of the tumor microenvironment (TME) in HRD and non-HRD TNBC samples may be helpful.

METHODS

Datasets from single-cell, spatial, and bulk RNA-sequencing were collected to explore the role of HRD in the development of TME at multiple scales. Based on the findings in the TME, machine learning algorithms were used to construct a response prediction model in eleven ICI therapy cohorts.

RESULTS

A more exhausted phenotype of T cells and a more tolerogenic phenotype of dendritic cells were found in the non-HRD group. HRD reprograms the predominant phenotype of cancer-associated fibroblasts (CAFs) from myofibroblastic CAFs to inflammatory-like CAFs. As interactions between myofibroblastic CAFs and other cells, DPP4-chemokines associated with reduced immune cell recruitment were unique in the non-HRD group. The prediction model based on DPP4-related genes had acceptable performance in predicting response, prognosis, and immune cell content. Higher HRD scores in bulk RNA-sequencing samples indicated more activated immune cell function, but not higher immune cell content, which may be affected by factors such as antigen-presenting capacity.

CONCLUSIONS

Based on multi-scale transcriptomics, our findings comprehensively reveal differences in the TME between HRD and non-HRD samples. Combining HRD with the prediction model or other methods for assessing immune cell content, may better predict response to ICIs in TNBC.

SECTM1 is upregulated in immuno-hot tumors and predicts immunotherapeutic efficacy in multiple cancers

Immune checkpoint inhibitors (ICIs) have transformed the management of advanced cancers. However, many patients could not benefit from ICIs therapy, and thus several biomarkers for therapeutic prediction have been uncovered. In this research, more than ten public and in-house cohorts were used to explore the predictive value and immunological correlations of secreted and transmembrane 1 (SECTM1) in cancers. SECTM1 expression was enhanced in tumors from patients with well immunotherapeutic responses in multiple cancers. In addition, SECTM1 was immuno-correlated in pan-cancer and enhanced in immuno-hot tumors. In vitro assays revealed that SECTM1 was upregulated by the IFN-γ/STAT1 signaling. Moreover, analysis of in-house immunotherapy cohorts suggested both tumor-expressed and circulating SECTM1 are promising biomarkers to predict therapeutic responses. Overall, this study reveals that SECTM1 is a biomarker of benefit to ICIs in cancer patients. Further studies including large-scale patients are needed to establish its utilization as a biomarker of benefit to ICIs.

Circulating Biomarkers for Prediction of Immunotherapy Response in NSCLC

Non-small cell lung cancer (NSCLC) constitutes the majority of the lung cancer population and the prognosis is poor. In recent years, immunotherapy has become the standard of care for advanced NSCLC patients as numerous trials demonstrated that immune checkpoint inhibitors (ICI) are more efficacious than conventional chemotherapy. However, only a minority of NSCLC patients benefit from this treatment. Therefore, there is an unmet need for biomarkers that could accurately predict response to immunotherapy. Liquid biopsy allows repeated sampling of blood-based biomarkers in a non-invasive manner for the dynamic monitoring of treatment response. In this review, we summarize the efforts and progress made in the identification of circulating biomarkers that predict immunotherapy benefit for NSCLC patients. We also discuss the challenges with future implementation of circulating biomarkers into clinical practice.

Spatial multi-omics revealed the impact of tumor ecosystem heterogeneity on immunotherapy efficacy in patients with advanced non-small cell lung cancer treated with bispecific antibody

BACKGROUND

Immunotherapy for malignant tumors has made great progress, but many patients do not benefit from it. The complex intratumoral heterogeneity (ITH) hindered the in-depth exploration of immunotherapy. Conventional bulk sequencing has masked intratumor complexity, preventing a more detailed discovery of the impact of ITH on treatment efficacy. Hence, we initiated this study to explore ITH at the multi-omics spatial level and to seek prognostic biomarkers of immunotherapy efficacy considering the presence of ITH.

METHODS

Using the segmentation strategy of digital spatial profiling (DSP), we obtained differential information on tumor and stromal regions at the proteomic and transcriptomic levels. Based on the consideration of ITH, signatures constructed by candidate proteins in different regions were used to predict the efficacy of immunotherapy.

RESULTS

Eighteen patients treated with a bispecific antibody (bsAb)-KN046 were enrolled in this study. The tumor and stromal areas of the same samples exhibited distinct features. Signatures consisting of 11 and 18 differentially expressed DSP markers from the tumor and stromal areas, respectively, were associated with treatment response. Furthermore, the spatially resolved signature identified from the stromal areas showed greater predictive power for bsAb immunotherapy response (area under the curve=0.838). Subsequently, our stromal signature was validated in an independent cohort of patients with non-small cell lung cancer undergoing immunotherapy.

CONCLUSION

We deciphered ITH at the spatial level and demonstrated for the first time that genetic information in the stromal region can better predict the efficacy of bsAb treatment.

TRIAL REGISTRATION NUMBER

NCT03838848.

Immune Pathways with Aging Characteristics Improve Immunotherapy Benefits and Drug Prediction in Human Cancer

(1) Background: Perturbation of immune-related pathways can make substantial contributions to cancer. However, whether and how the aging process affects immune-related pathways during tumorigenesis remains largely unexplored. (2) Methods: Here, we comprehensively investigated the immune-related genes and pathways among 25 cancer types using genomic and transcriptomic data. (3) Results: We identified several pathways that showed aging-related characteristics in various cancers, further validated by conventional aging-related gene sets. Genomic analysis revealed high mutation burdens in cytokines and cytokines receptors pathways, which were strongly correlated with aging in diverse cancers. Moreover, immune-related pathways were found to be favorable prognostic factors in melanoma. Furthermore, the expression level of these pathways had close associations with patient response to immune checkpoint blockade therapy in melanoma and non-small cell lung cancer. Applying a net-work-based method, we predicted immune- and aging-related genes in pan-cancer and utilized these genes for potential immunotherapy drug discovery. Mapping drug target data to our top-ranked genes identified potential drug targets, FYN, JUN, and SRC. (4) Conclusions: Taken together, our systematic study helped interpret the associations among immune-related pathways, aging, and cancer and could serve as a resource for promoting clinical treatment.

The prognostic impact of tumor-infiltrating B lymphocytes in patients with solid malignancies: A systematic review and meta-analysis

This study reviewed the prognostic effect of tumor-infiltrating B lymphocytes (TIBLs) on solid malignancies, to determine the potential role of TIBLs in predicting cancer patient's prognosis and their response to immunotherapy. A total of 45 original papers involving 11,099 individual patients were included in this meta-analysis covering 7 kinds of cancer. The pooled results suggested that high levels of TIBLs were correlated with favorable OS in lung, esophageal, gastric, colorectal, liver, and breast cancer; improved RFS in lung cancer; and improved DFS in gastrointestinal neoplasms. Additionally, TIBLs were significantly correlated with negative lymphatic invasion in gastric cancer, small tumor size in hepatocellular carcinoma, and negative distant metastasis in colorectal cancer. Additionally, TIBLs were reported as a discriminative feature of patients treated with immunotherapy with improved survival. We concluded that TIBLs play a favorable prognostic role among the common solid malignancie, providing theoretical evidence for further prognosis prediction for solid tumors.

Genome-wide DNA methylation profile analysis identifies an individualized predictive signature for melanoma immune response

PURPOSE

The current evaluation methods for tumor infiltrating lymphocytes (TILs), particularly CD8 + TILs, mainly rely on semiquantitative immunohistochemistry with high variability. We aimed to construct an individualized DNA methylation-based signature for CD8 + TILs (CD8 + MeTIL) that may characterize melanoma immune microenvironment and guide therapeutic selection.

METHODS

The transcriptome profiles and DNA methylation data of 457 melanoma patients from The Cancer Genome Atlas (TCGA) database were analyzed. Differential methylation analysis between groups with high and low CD8 + TILs was performed to select differentially methylated positions (DMPs) and define CD8 + MeTIL. The prognostic value of CD8 + MeTIL and its predictive value for immunotherapy response were investigated using multiple melanoma cohorts.

RESULTS

We successfully constructed the CD8 + MeTIL signature based on four DMPs. The survival analyses showed that higher CD8 + MeTIL score was associated with worse survival outcomes in TCGA-SKCM and GSE144487 cohorts. The ROC curve for the predictive analysis revealed that the survival prediction of CD8 + MeTIL score was superior compared with CD8 + TILs (CIBERSORT) and CD8B mRNA expression. Furthermore, we founded that tumors with higher CD8 + MeTIL score were marked with immunosuppressive characteristics, including low immune score and downregulated immune-related pathways. More importantly, the CD8 + MeTIL score showed a potential predictive value for the benefit from immunotherapy in two published cohorts. When combined CD8 + MeTIL with PD-L1 expression, the patient classification showed significantly different immunotherapy response rates and long-term survival outcomes.

CONCLUSIONS

The CD8 + MeTIL signature might be as a novel method to evaluate CD8 + TILs and guide immunotherapy approaches.

SUMOylation patterns and signature characterize the tumor microenvironment and predict prognosis in lung adenocarcinoma

Background: Recent studies have revealed that SUMOylation modifications are involved in various biological processes, including cancer development and progression. However, the precise role of SUMOylation in lung adenocarcinoma (LUAD), especially in the tumor immune microenvironment, is not yet clear. Methods: We identified SUMOylation patterns by unsupervised consensus clustering based on the expression of SUMOylation regulatory genes. The tumor microenvironment in lung adenocarcinoma was analyzed using algorithms such as GSVA and ssGSEA. Key genes of SUMOylation patterns were screened for developing a SUMOylation scoring model to assess immunotherapy and chemotherapy responses in lung adenocarcinoma patients. Experiments were conducted to validate the differential expression of model genes in lung adenocarcinoma. Finally, we constructed a nomogram based on the SUMOylation score to assess the prognosis of individual lung adenocarcinoma patients. Results: Two patterns of SUMOylation were identified, namely, SUMO-C1, which showed anti-tumor immune phenotype, and SUMO-C2, which showed immunosuppressive phenotype. Different genomic subtypes were also identified; subtype gene-T1 exhibited a reciprocal restriction between the immune microenvironment and stromal microenvironment. High SUMOylation scores were indicative of poor lung adenocarcinoma prognosis. SUMOylation score was remarkably negatively correlated with the infiltration of anti-tumor immune cells, and significantly positively correlated with immune cells promoting immune escape and immune suppression. In addition, patients with low scores responded better to immunotherapy. Therefore, the developed nomogram has a high prognostic predictive value. Conclusion: The SUMOylation patterns can well discriminate the tumor microenvironment features of lung adenocarcinoma, especially the immune cell infiltration status. The SUMOylation score can further assess the relationship between SUMOylation and immune cell crosstalk and has significant prognostic value and can be used to predict immunotherapy and chemotherapy response in patients with lung adenocarcinoma.

Receptor-ligand pair typing and prognostic risk model of response or resistance to immune checkpoint inhibitors in lung adenocarcinoma

Introduction

Currently, programmed cell death-1 (PD-1)-targeted treatment is ineffective for a sizable minority of patients, and drug resistance still cannot be overcome.

Methods

To explore the mechanisms of immunotherapy and identify new therapeutic opportunities in lung adenocarcinoma (LUAD), data from patients who did and did not respond to the anti-PD-1 treatment were evaluated using single-cell RNA sequencing, and bulk RNA sequencing were collected.

Results

We investigated the gene expression that respond or not respond to immunotherapy in diverse cell types and revealed transcriptional characteristics at the single-cell level. To ultimately explore the molecular response or resistance to anti-PD-1 therapy, cell-cell interactions were carried out to identify the different LRIs (ligand-receptor interactions) between untreated patients vs. no-responders, untreated patients vs. responders, and responders vs. non-responders. Next, two molecular subgroups were proposed based on 73 LRI genes, and subtype 1 had a poor survival status and was likely to be the immunosuppressive tumor subtype. Furthermore, based on the LASSO Cox regression analysis results, we found that TNFSF13, AXL, KLRK1, FAS, PROS1, and CDH1 can be distinct prognostic biomarkers, immune infiltration levels, and responses to immunotherapy in LUAD.

Discussion

Altogether, the effects of immunotherapy were connected to LRIs scores, indicating that potential medications targeting these LRIs could contribute to the clinical benefit of immunotherapy. Our integrative omics analysis revealed the mechanisms underlying the anti-PD-1 therapy response and offered abundant clues for potential strategies to improve precise diagnosis and immunotherapy.

Identification of genes and cellular response factors related to immunotherapy response in mismatch repair-proficient colorectal cancer: a bioinformatics analysis

Background

Mismatch repair-proficient (pMMR) colorectal cancers (CRCs) are thought to be primarily resistant to immune checkpoint inhibitor (ICI) monotherapy. However, recent clinical trials have reported that early-to-mid stage (non-metastatic) CRC responds well to ICI monotherapy. We hypothesized that the efficacy of immunotherapy is linked to a series of gene expression profiles that can characterize the pMMR CRC disease stage.

Methods

Using The Cancer Genome Atlas (TCGA) CRC data sets, we first investigated transcriptomic features that continuously changed (were continuously upregulated or downregulated) with pMMR CRC disease-stage progression. We defined these gene sets as stage-associated genes. The deconvolution algorithm then enriched these genes with the dynamic changes in the cell type populations of the CRC tumor microenvironment (TME). Finally, the stage-associated genes were cross-referenced to the current transcriptome profile data on ICI treatment of pMMR CRC, which revealed the gene set specifying an effective pMMR tumor response.

Results

In total, 774 genes were found to increase in expression and 845 genes to decrease in expression as the stage increased. Using deconvolution methods, we discovered 2 major disease stage-associated alterations in the cellular composition of pMMR CRCs, including changes in cell types involved in host immune responses and tumor cell metastasis. The central memory CD8⁺ T cell population decreased as the pMMR CRC disease stage increased, but the endothelial cell populations associated with proliferation and metastasis increased. Using a different cell type annotation set (LM22), we discovered that as the disease progressed, M1 macrophages and CD8⁺ T cells decreased in the TME. In mismatch repair-deficient patients with CRC, however, such a decrease was not observed. Finally, we identified 27 signature genes that can be used to assess ICI efficacy in treatment-naïve patients with pMMR CRC.

Conclusions

The current study sought to identify the underlying molecular mechanisms, pathways, and cell landscapes that explain why early-to-mid stage pMMR CRC responds well to ICI treatment. This analysis might be valuable for the selection of patients who might benefit from immunotherapeutic strategies.

ACAP1 Deficiency Predicts Inferior Immunotherapy Response in Solid Tumors

BACKGROUND

ACAP1 plays a key role in endocytic recycling, which is essential for the normal function of lymphocytes. However, the expression and function of ACAP1 in lymphocytes have rarely been studied.

METHODS

Large-scale genomic data, including multiple bulk RNA-sequencing datasets, single-cell sequencing datasets, and immunotherapy cohorts, were exploited to comprehensively characterize ACAP1 expression, regulation, and function. Gene set enrichment analysis (GSEA) was used to uncover the pathways associated with ACAP1 expression. Eight algorithms, including TIMER, CIBERSORT, CIBERSORT-ABS, QUANTISEQ, xCELL, MCPCOUNTER, EPIC, and TIDE, were applied to estimate the infiltrating level of immune cells. Western blotting, qPCR, and ChIP-PCR were used to validate the findings from bioinformatic analyses. A T-cell co-culture killing assay was used to investigate the function of ACAP1 in lymphocytes.

RESULTS

ACAP1 was highly expressed in immune-related tissues and cells and minimally in other tissues. Moreover, single-cell sequencing analysis in tumor samples revealed that ACAP1 is expressed primarily in tumor-infiltrating lymphocytes (TILs), including T, B, and NK cells. ACAP1 expression is negatively regulated by promoter DNA methylation, with its promoter hypo-methylated in immune cells but hyper-methylated in other cells. Furthermore, SPI1 binds to the ACAP1 promoter and positively regulates its expression in immune cells. ACAP1 levels positively correlate with the infiltrating levels of TILs, especially CD8+ T cells, across a broad range of solid cancer types. ACAP1 deficiency is associated with poor prognosis and immunotherapeutic response in multiple cancer types treated with checkpoint blockade therapy (ICT). Functionally, the depletion of ACAP1 by RNA interference significantly impairs the T cell-mediated killing of tumor cells.

CONCLUSIONS

Our study demonstrates that ACAP1 is essential for the normal function of TILs, and its deficiency indicates an immunologically "cold" status of tumors that are resistant to ICT.

Analysis on methylation and expression of PSMB8 and its correlation with immunity and immunotherapy in lung adenocarcinoma

Aim: To find biomarkers for immunity and immunotherapy in lung adenocarcinoma (LUAD) through multiomics analysis. Materials & methods: The multiomics data of patients with LUAD were downloaded from the TCGA and GEO databases. CIBERSORT, quanTIseq, ESTIMATEScore, k-means clustering, gene set enrichment analysis, gene set variation analysis, immunophenoscore and logistic regression were used in this study. Results: PSMB8 HypoMet-HighExp group patients have more active immune-related pathways, more antitumor immune cells, less protumor immune cells, higher immunophenoscore and longer progression-free survival of immune checkpoint inhibitor therapy than HyperMet-LowExp group. In multivariate analysis, PSMB8 showed an independent value. Conclusion: The combination of DNA methylation and mRNA expression of PSMB8 could independently distinguish types of tumor immune microenvironment and predict programmed cell death protein 1/programmed cell death-ligand 1 inhibitors' effects in patients with LUAD.

DLX2 Is a Potential Immune-Related Prognostic Indicator Associated with Remodeling of Tumor Microenvironment in Lung Squamous Cell Carcinoma: An Integrated Bioinformatical Analysis

BACKGROUND

It is still an unmet clinical need to identify potent biomarkers for immunotherapy on patients with lung squamous cell carcinoma (LUSC).

METHODS

In this study, we explored the differentially expressed genes (DEGs) that were simultaneously correlated with four pathways (i.e. CD8+ αβT cell proliferation/differentiation/activation pathways and ferroptosis pathway) and possibly related to the remodeling of tumor microenvironment via the TCGA-LUSC dataset. Besides, four GEO datasets (GSE157009, GSE157010, GSE19188, and GSE126045) and IMvigor210 dataset were utilized for confirmation and validation.

RESULTS

The co-downregulated DEG DLX2 was selected for further analysis. Function enrichment analysis revealed that low-expression of DLX2 was closely related to various immune-related pathways like T/B/NK cell mediated immunity, interferon gamma/alpha response, and various autoimmune disease. DLX2-downregulated group was enriched in more immune-activating cells and lower tumor immune dysfunction and exclusion (TIDE) score. Via the Cancer Immunome Atlas (TCIA) database, lower expression of DLX2 was also found to be associated with better IPS score of PD-1/PD-L1 blockade (p < 0.001) as well as CTLA-4 combined with PD-1/PD-L1 blockade (p < 0.001). Furthermore, patients in DLX2-low group were found to have significant longer median OS than those in DLX2-high group in IMvigor210 dataset (10.8 vs 7.4 months; hazard ratio [HR]=0.74, 95% confidence interval [95%CI] 0.57-0.96; p = 0.024).

CONCLUSIONS

Our study on an integrated bioinformatical analysis implied that DLX2 could be served as a promising indicator for remodeling tumor microenvironment status and predicting ICI response of patients with LUSC.

Integrated multi-omics approach revealed cellular senescence landscape

Cellular senescence is a complex multifactorial biological phenomenon that plays essential roles in aging, and aging-related diseases. During this process, the senescent cells undergo gene expression altering and chromatin structure remodeling. However, studies on the epigenetic landscape of senescence using integrated multi-omics approaches are limited. In this research, we performed ATAC-seq, RNA-seq and ChIP-seq on different senescent types to reveal the landscape of senescence and identify the prime regulatory elements. We also obtained 34 key genes and deduced that NAT1, PBX1 and RRM2, which interacted with each other, could be the potential markers of aging and aging-related diseases. In summary, our work provides the landscape to study accessibility dynamics and transcriptional regulations in cellular senescence. The application of this technique in different types of senescence allows us to identify the regulatory elements responsible for the substantial regulation of transcription, providing the insights into molecular mechanisms of senescence.

Molecular subtypes based on cuproptosis-related genes and immune profiles in lung adenocarcinoma

Background: Recent studies have identified several molecular subtypes of lung adenocarcinoma (LUAD) that have different prognoses to help predict the efficacy of immunotherapy. However, the prognostic prediction is less than satisfactory. Alterations in intracellular copper levels may affect the tumor immune microenvironment and are linked to cancer progression. Previous studies have identified some genes related to cuproptosis. The characteristics of the cuproptosis molecular subtypes have not been thoroughly studied in LUAD. Methods: The transcriptomic data and clinical information of 632 LUAD patients were used to investigate the LUAD molecular subtypes that are associated with the cuproptosis-related genes (CRGs), the tumor immune microenvironment, and stemness. The cuproptosis score was constructed using univariate Cox regression and the minor absolute shrinkage and selection operator (LASSO) to quantify the prognostic characteristics. Results: Three different molecular subtypes related to cuproptosis, with different prognoses, were identified in LUAD. Cluster A had the highest cuproptosis score and the worst prognosis. Patients in the high cuproptosis score group had a higher somatic mutation frequency and stemness scores. Patients in the low cuproptosis score group had more immune infiltration and better prognosis. Conclusion: Molecular subtypes of LUAD based on CRGs reflect the differences in LUAD patients. The cuproptosis score can be used as a promising biomarker, which is of great significance to distinguish the relationship between cuproptosis and the immune microenvironment. The cuproptosis signature based on the cuproptosis score and clinical characteristics of individual patients will be useful for guiding immunotherapy in LUAD.

Integrative Analysis of Bioinformatics and Machine Learning Algorithms Identifies a Novel Diagnostic Model Based on Costimulatory Molecule for Predicting Immune Microenvironment Status in Lung Adenocarcinoma

Costimulatory molecules are an indispensable signal for activating immune cells. However, the features of many costimulatory molecule genes (CMGs) in lung adenocarcinoma (LUAD) are poorly understood. This study systematically explored expression patterns of CMGs in the tumor immune microenvironment (TIME) status of patients with LUAD. Their expression profiles were downloaded from The Cancer Genome Atlas and the Gene Expression Omnibus databases. Two robust TIME subtypes ("hot" and "cold") were classified by K-means clustering and estimation of stromal and immune cells in malignant tumor tissues using expression data. The "hot" subtype presented higher infiltration in activated immune cells and enrichments in the immune cell receptor signaling pathway and adaptive immune response. Three CMGs (CD80, LTB, and TNFSF8) were screened as final diagnostic markers by means of Least Absolute Shrinkage Selection Operator and Support Vector Machine-Recursive Feature Elimination algorithms. Accordingly, the diagnostic nomogram for predicting individualized TIME status showed satisfactory diagnostic accuracy in The Cancer Genome Atlas training cohort as well as GSE31210 and GSE180347 validation cohorts. Immunohistochemistry staining of 16 specimens revealed an apparently positive correlation between the expression of CMG biomarkers and pathologic response to immunotherapy. Thus, this diagnostic nomogram provided individualized predictions in TIME status of LUAD patients with good predictive accuracy, which could serve as a potential tool for identifying ideal candidates for immunotherapy.

Construction of a predictive model for immunotherapy efficacy in lung squamous cell carcinoma based on the degree of tumor-infiltrating immune cells and molecular typing

Background

To construct a predictive model of immunotherapy efficacy for patients with lung squamous cell carcinoma (LUSC) based on the degree of tumor-infiltrating immune cells (TIIC) in the tumor microenvironment (TME).

Methods

The data of 501 patients with LUSC in the TCGA database were used as a training set, and grouped using non-negative matrix factorization (NMF) based on the degree of TIIC assessed by single-sample gene set enrichment analysis (GSEA). Two data sets (GSE126044 and GSE135222) were used as validation sets. Genes screened for modeling by least absolute shrinkage and selection operator (LASSO) regression and used to construct a model based on immunophenotyping score (IPTS). RNA extraction and qPCR were performed to validate the prognostic value of IPTS in our independent LUSC cohort. The receiver operating characteristic (ROC) curve was constructed to determine the predictive value of the immune efficacy. Kaplan–Meier survival curve analysis was performed to evaluate the prognostic predictive ability. Correlation analysis and enrichment analysis were used to explore the potential mechanism of IPTS molecular typing involved in predicting the immunotherapy efficacy for patients with LUSC.

Results

The training set was divided into a low immune cell infiltration type (C1) and a high immune cell infiltration type (C2) by NMF typing, and the IPTS molecular typing based on the 17-gene model could replace the results of the NMF typing. The area under the ROC curve (AUC) was 0.82. In both validation sets, the IPTS of patients who responded to immunotherapy were significantly higher than those who did not respond to immunotherapy (P = 0.0032 and P = 0.0451), whereas the AUC was 0.95 (95% CI = 1.00–0.84) and 0.77 (95% CI = 0.58–0.96), respectively. In our independent cohort, we validated its ability to predict the response to cancer immunotherapy, for the AUC was 0.88 (95% CI = 1.00–0.66). GSEA suggested that the high IPTS group was mainly involved in immune-related signaling pathways.

Conclusions

IPTS molecular typing based on the degree of TIIC in the TME could well predict the efficacy of immunotherapy in patients with LUSC with a certain prognostic value.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12967-022-03565-7.

Prediction of prognosis, immunogenicity and efficacy of immunotherapy based on glutamine metabolism in lung adenocarcinoma

BACKGROUND

Glutamine (Gln) metabolism has been reported to play an essential role in cancer. However, a comprehensive analysis of its role in lung adenocarcinoma is still unavailable. This study established a novel system of quantification of Gln metabolism to predict the prognosis and immunotherapy efficacy in lung cancer. Further, the Gln metabolism in tumor microenvironment (TME) was characterized and the Gln metabolism-related genes were identified for targeted therapy.

METHODS

We comprehensively evaluated the patterns of Gln metabolism in 513 patients diagnosed with lung adenocarcinoma (LUAD) based on 73 Gln metabolism-related genes. Based on differentially expressed genes (DEGs), a risk model was constructed using Cox regression and Lasso regression analysis. The prognostic efficacy of the model was validated using an individual LUAD cohort form Shandong Provincial Hospital, an integrated LUAD cohort from GEO and pan-cancer cohorts from TCGA databases. Five independent immunotherapy cohorts were used to validate the model performance in predicting immunotherapy efficacy. Next, a series of single-cell sequencing analyses were used to characterize Gln metabolism in TME. Finally, single-cell sequencing analysis, transcriptome sequencing, and a series of in vitro experiments were used to explore the role of EPHB2 in LUAD.

RESULTS

Patients with LUAD were eventually divided into low- and high-risk groups. Patients in low-risk group were characterized by low levels of Gln metabolism, survival advantage, "hot" immune phenotype and benefit from immunotherapy. Compared with other cells, tumor cells in TME exhibited the most active Gln metabolism. Among immune cells, tumor-infiltrating T cells exhibited the most active levels of Gln metabolism, especially CD8 T cell exhaustion and Treg suppression. EPHB2, a key gene in the model, was shown to promote LUAD cell proliferation, invasion and migration, and regulated the Gln metabolic pathway. Finally, we found that EPHB2 was highly expressed in macrophages, especially M2 macrophages. It may be involved in the M2 polarization of macrophages and mediate the negative regulation of M2 macrophages in NK cells.

CONCLUSION

This study revealed that the Gln metabolism-based model played a significant role in predicting prognosis and immunotherapy efficacy in lung cancer. We further characterized the Gln metabolism of TME and investigated the Gln metabolism-related gene EPHB2 to provide a theoretical framework for anti-tumor strategy targeting Gln metabolism.

Statin shapes inflamed tumor microenvironment and enhances immune checkpoint blockade in non-small cell lung cancer

Immune checkpoint blockade (ICB) therapy has achieved breakthroughs in the treatment of advanced non-small cell lung cancer (NSCLC). Nevertheless, the low response due to immuno-cold tumor microenvironment (TME) largely limits the application of ICB therapy. Based on the glycolytic/cholesterol synthesis axis, a stratification framework for EGFR wild-type NSCLC was developed to summarize the metabolic features of immuno-cold and immuno-hot tumors. The cholesterol subgroup displays the worst prognosis in immuno-cold NSCLC with significant enrichment of the cholesterol gene signature, indicating targeting cholesterol synthesis is essential for the therapy for immuno-cold NSCLC. Statin, the inhibitor for cholesterol synthesis, can suppress the aggressiveness of NSCLC in vitro and in vivo and also drastically reverse immuno-cold to an inflamed phenotype in vivo which exhibited a higher response to ICB therapy. Moreover, both our in-house data and meta-analysis further support that statin can significantly enhance ICB efficacy. In terms of preliminary mechanisms, statin could transcriptionally inhibit PD-L1 expression and induce ferroptosis in NSCLC cells. Overall, we reveal the significance of cholesterol synthesis in NSCLC and demonstrate the improved therapeutic efficacy of ICB in combination with statin. These findings could provide a innovative clinical insight to treat NSCLC patients with immuno-cold tumors.

Single-cell transcriptomic analysis reveals circadian rhythm disruption associated with poor prognosis and drug-resistance in lung adenocarcinoma

Circadian rhythm disruption (CRD) represents a major contributor to tumor proliferation. Nonetheless, the role of CRD in the clinical prediction of cancer outcomes has not been well studied. In this study, we developed a computational algorithm, which was implemented in an open-source R package CRDscore, to define the intratumoral status of circadian disruption in three representative single-cell RNA-seq data sets of lung adenocarcinoma. We found that the malignant cells with high CRDscore were characterized by activation of glycolysis and epithelial-mesenchymal transition pathways. Furthermore, cell communication analysis indicated that CRD played a pivotal role in T cell exhaustion, which may be responsible for the poor prognosis of the malignancy. We then validated the findings with public bulk transcriptome datasets involving 22 cancer types. Cox regression analysis revealed that the CRDscore was a valuable prognostic biomarker. A model containing 23 circadian-related genes performed well in predicting immunotherapeutic outcomes in 14 independent cohorts. Importantly, decreased CRDscore was detect by RNA sequencing on H1299 cells with melatonin treatment. Meanwhile, the cells downregulated the expression level of SNAIL and TWIST, which contributed to an invasive phenotype. In conclusion, this study provides a novel computational framework for characterizing CRD status using single-cell transcriptomic data and further confirmed the molecular mechanisms underlying metabolic reprogramming and T cell exhaustion under CRD. The better understanding of the mechanisms may provide new possibilities for incorporating "anticancer approaches based on circadian clocks" into the treatment protocols of precision medicine.

Construction and Validation of an Epigenetic Regulator Signature as A Novel Biomarker For Prognosis, Immunotherapy, And Chemotherapy In Hepatocellular Carcinoma

Background

Epigenetic modification regulates various aspects of cancer biology, from tumor growth and invasion to immune microenvironment modulation. Whether epigenetic regulators (EGRs) can decide tumor malignant degree and risk of immune evasion in liver hepatocellular carcinoma (LIHC) remains unclear.

Method

An EGR signature called “EGRscore” was constructed based on bulk RNA-seq data of EGR in hepatocellular carcinoma (HCC). The correlation between EGRscore and overall survival (OS) was validated in HCC cohorts and other tumor cohorts. Mutation profiles, copy number alterations (CNAs), enriched pathways, and response to immunotherapy and chemotherapy were compared between EGRscore-high and EGRscore-low patients.

Results

We found that EGRscore was associated with OS in HCC as well as several tumors including glioma, uveal melanoma (UVM), and kidney tumors. A mechanism study demonstrated that the distinct mutation profile of TP53 was present in EGRscore-high and EGRscore-low patients. Meanwhile, EGRscore-low patients were characterized with immune cells that promote killing tumors. Furthermore, EGRscore was associated with genes regulating drug resistance in HCC. Finally, we indicated that EGRscore-low patients had higher response rates to immunotherapy and targeted therapy.

Conclusions

EGRscore could be used to distinguish OS, tumor progression, mutation pattern, and immune microenvironment. The present study contributes to improving hepatocellular carcinoma patient prognosis and predicting response to immunotherapy.

A novel prognostic signature of metastasis-associated genes and personalized therapeutic strategy for lung adenocarcinoma patients

Lung adenocarcinoma (LUAD) is a highly invasive and metastatic malignant tumor with high morbidity and mortality. This study aimed to construct a prognostic signature for LUAD patients based on metastasis-associated genes (MAGs). RNA expression profiles were downloaded from the Cancer Genome Atlas (TCGA) database. RRA method was applied to identify differentially expressed MAGs. A total of 192 significantly robust MAGs were determined among seven GEO datasets. MAGs were initially selected through the Lasso Cox regression analysis and 6 MAGs were included to construct a prognostic signature model. Transcriptome profile, patient prognosis, correlation between the risk score and clinicopathological features, immune cell infiltration characteristics, immunotherapy sensitivity and chemotherapy sensitivity differed between low- and high-risk groups after grouping according to median risk score. The reliability and applicability of the signature were further validated in the GSE31210, GSE50081 and GSE68465 cohort. CMap predicted 62 small molecule drugs on the base of the prognostic MAGs. Targeted drug staurosporine had hydrogen bonding with Gln-172 of SLC2A1, which is one of MAGs. Staurosporine could inhibit cell migration in A549 and H1299. We further verified mRNA and protein expression of 6 MAGs in A549 and H1299. The signature can serve as a promising prognostic tool and may provide a novel personalized therapeutic strategy for LUAD patients.

HLA class II molecule HLA-DRA identifies immuno-hot tumors and predicts the therapeutic response to anti-PD-1 immunotherapy in NSCLC

Background

Immune checkpoint blockade (ICB) only works well for a certain subset of patients with non-small cell lung cancer (NSCLC). Therefore, biomarkers for patient stratification are desired, which can suggest the most beneficial treatment.

Methods

In this study, three datasets (GSE126044, GSE135222, and GSE136961) of immunotherapy from the Gene Expression Omnibus (GEO) database were analyzed, and seven intersected candidates were extracted as potential biomarkers for ICB followed by validation with The Cancer Genome Atlas (TCGA) dataset and the in-house cohort data.

Results

Among these candidates, we found that human leukocyte antigen-DR alpha (HLA-DRA) was downregulated in NSCLC tissues and both tumor and immune cells expressed HLA-DRA. In addition, HLA-DRA was associated with an inflamed tumor microenvironment (TME) and could predict the response to ICB in NSCLC. Moreover, we validated the predictive value of HLA-DRA in immunotherapy using an in-house cohort. Furthermore, HLA-DRA was related to the features of inflamed TME in not only NSCLC but also in most cancer types.

Conclusion

Overall, HLA-DRA could be a promising biomarker for guiding ICB in NSCLC.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12885-022-09840-6.

Identification and Validation of Immune-Related Long Non-Coding RNA Signature for Predicting Immunotherapeutic Response and Prognosis in NSCLC Patients Treated With Immunotherapy

Background

Numerous studies have reported that long non-coding RNAs (lncRNAs) play important roles in immune-related pathways in cancer. However, immune-related lncRNAs and their roles in predicting immunotherapeutic response and prognosis of non-small cell lung cancer (NSCLC) patients treated with immunotherapy remain largely unexplored.

Methods

Transcriptomic data from NSCLC patients were used to identify novel lncRNAs by a custom pipeline. ImmuCellAI was utilized to calculate the infiltration score of immune cells. The marker genes of immunotherapeutic response-related (ITR)-immune cells were used to identify immune-related (IR)-lncRNAs. A co-expression network was constructed to determine their functions. LASSO and multivariate Cox analyses were performed on the training set to construct an immunotherapeutic response and immune-related (ITIR)-lncRNA signature for predicting the immunotherapeutic response and prognosis of NSCLC. Four independent datasets involving NSCLC and melanoma patients were used to validate the ITIR-lncRNA signature.

Results

In total, 7,693 novel lncRNAs were identified for NSCLC. By comparing responders with non-responders, 154 ITR-lncRNAs were identified. Based on the correlation between the marker genes of ITR-immune cells and lncRNAs, 39 ITIR-lncRNAs were identified. A co-expression network was constructed and the potential functions of 38 ITIR-lncRNAs were annotated, most of which were related to immune/inflammatory-related pathways. Single-cell RNA-seq analysis was performed to confirm the functional prediction results of an ITIR-lncRNA, LINC01272. Four-ITIR-lncRNA signature was identified and verified for predicting the immunotherapeutic response and prognosis of NSCLC. Compared with non-responders, responders had a lower risk score in both NSCLC datasets (P<0.05). NSCLC patients in the high-risk group had significantly shorter PFS/OS time than those in the low-risk group in the training and testing sets (P<0.05). The AUC value was 1 of responsiveness in the training set. In melanoma validation datasets, patients in the high-risk group also had significantly shorter OS/PFS time than those in the low-risk group (P<0.05). The ITIR-lncRNA signature was an independent prognostic factor (P<0.001).

Conclusion

Thousands of novel lncRNAs in NSCLC were identified and characterized. In total, 39 ITIR-lncRNAs were identified, 38 of which were functionally annotated. Four ITIR-lncRNAs were identified as a novel ITIR-lncRNA signature for predicting the immunotherapeutic response and prognosis in NSCLC patients treated with immunotherapy.

Single-cell N6-methyladenosine regulator patterns guide intercellular communication of tumor microenvironment that contribute to colorectal cancer progression and immunotherapy

BACKGROUND

N6-methyladenosine (m⁶A) RNA methylation plays a critical role in key genetic events for various cancers; yet, how m⁶A functions within the tumor microenvironment (TME) remains to be elucidated.

METHODS

A total of 65,362 single cells from single-cell RNA-seq data derived from 33 CRC tumor samples were analyzed by nonnegative matrix factorization (NMF) for 23 m⁶A RNA methylation regulators. CRC and Immunotherapy cohorts from public repository were used to determine the prognosis and immune response of TME clusters.

RESULTS

The fibroblasts, macrophages, T and B cells were respectively grouped into 4 to 5 subclusters and then classified according to various biological processes and different marker genes. Furthermore, it revealed that the m⁶A RNA methylation regulators might be significantly related to the clinical and biological features of CRC, as well as the pseudotime trajectories of main TME cell types. Bulk-seq analysis suggested that these m⁶A-mediated TME cell subclusters had significant prognostic value for CRC patients and distinguished immune response for patients who underwent ICB therapy, especially for the CAFs and macrophages. Notably, CellChat analysis revealed that RNA m⁶A methylation-associated cell subtypes of TME cells manifested diverse and extensive interaction with tumor epithelial cells. Further analysis showed that ligand-receptor pairs, including MIF -  (CD74 + CXCR4), MIF -  (CD74 + CD44), MDK-NCL and LGALS9 - CD45, etc. mediated the communication between m⁶A associated subtypes of TME cells and tumor epithelial cells.

CONCLUSIONS

Taken together, our study firstly revealed the m⁶A methylation mediated intercellular communication of the tumor microenvironment in the regulation of tumor growth and antitumor immunomodulatory processes.

Hyper-Methylated Hub Genes of T-Cell Receptor Signaling Predict a Poor Clinical Outcome in Lung Adenocarcinoma

Background

Immune checkpoint inhibitors (ICIs) emerge as the first-line treatment of lung adenocarcinoma (LUAD); selection of subpopulations acquiring clinical benefit is required. Associations between epigenetic modulation of tumor microenvironment (TME) and clinical outcome are far from clear. We focused on immune-related genes closely regulated by DNA methylation to identify the potential clinical outcome indicators.

Methods

We systematically calculated immunophenotype score (IMpS) and classified immunophenotypes based on seven TME features in three independent cohorts. The overlapping of differential expressed genes and methylated probes targeted genes was regarded as genes closely regulated by DNA methylation. Then, probe/gene pairs which highly correlated with each other and IMpS were identified and named as immune-related probe/gene pairs (mIMg). Prognostic mIMg were selected and verified in seven independent validation cohorts.

Results

Three immune phenotypes were clustered, and similar results were obtained in the three independent training cohorts. C2 displayed as an immunologically hot phenotype, whereas C3 corresponded with immunologically cold phenotype. Average methylation level was decreased from C2 to C3 (C2 > C1 > C3). Similarly, ICIs nonresponders showed global hypo-methylation compared with responders. Genes in mIMg were mainly enriched, especially in T-cell receptor activation, and repressed in noninflamed TME by hyper-methylation. Among mIMg, low expression and hyper-methylation of CD247, LCK, and PSTPIP1 were risk factors of overall survival (OS). ICIs nonresponders were more likely to be hyper-methylated in the three genes. By integrating with the oncogenes status, we demonstrated that EGFR wt and SRGN overexpressed patients were associated with chronic inflammation and immune evasion, showing an immunologically hot phenotype, which might lead to the short OS but derive clinical benefit from ICIs.

Conclusions

This study identifies hyper-methylation and concurrent repression of CD247, LCK, PSTPIP1 as immune negative indicators and risk factors for prognosis in LUAD. Moreover, EGFR/SRGN axis may participate in immune modification to influence ICIs response and clinical outcome in LUAD.

Integrated Analysis of MATH-Based Subtypes Reveals a Novel Screening Strategy for Early-Stage Lung Adenocarcinoma

Lung adenocarcinoma (LUAD) is a frequently diagnosed cancer type, and many patients have already reached an advanced stage when diagnosed. Thus, it is crucial to develop a novel and efficient approach to diagnose and classify lung adenocarcinoma at an early stage. In our study, we combined in silico analysis and machine learning to develop a new five-gene–based diagnosis strategy, which was further verified in independent cohorts and in vitro experiments. Considering the heterogeneity in cancer, we used the MATH (mutant-allele tumor heterogeneity) algorithm to divide patients with early-stage LUAD into two groups (C1 and C2). Specifically, patients in C2 had lower intratumor heterogeneity and higher abundance of immune cells (including B cell, CD4 T cell, CD8 T cell, macrophage, dendritic cell, and neutrophil). In addition, patients in C2 had a higher likelihood of immunotherapy response and overall survival advantage than patients in C1. Combined drug sensitivity analysis (CTRP/PRISM/CMap/GDSC) revealed that BI-2536 might serve as a new therapeutic compound for patients in C1. In order to realize the application value of our study, we constructed the classifier (to classify early-stage LUAD patients into C1 or C2 groups) with multiple machine learning and bioinformatic analyses. The 21-gene–based classification model showed high accuracy and strong generalization ability, and it was verified in four independent validation cohorts. In summary, our research provided a new strategy for clinicians to make a quick preliminary assisting diagnosis of early-stage LUAD and make patient classification at the intratumor heterogeneity level. All data, codes, and study processes have been deposited to Github and are available online.

A Prognostic Risk Score Based on Hypoxia-, Immunity-, and Epithelialto-Mesenchymal Transition-Related Genes for the Prognosis and Immunotherapy Response of Lung Adenocarcinoma

Background: Lung adenocarcinoma (LUAD), the most common subtype of non-small cell lung cancer (NSCLC), is associated with poor prognosis. However, current stage-based clinical methods are insufficient for survival prediction and decision-making. This study aimed to establish a novel model for evaluating the risk of LUAD based on hypoxia, immunity, and epithelial-mesenchymal transition (EMT) gene signatures.

Methods: In this study, we used data from TCGA-LUAD for the training cohort and GSE68465 and GSE72094 for the validation cohorts. Immunotherapy datasets GSE135222, GSE126044, and IMvigor210 were obtained from a previous study. Using bioinformatic and machine algorithms, we established a risk model based on hypoxia, immune, and EMT gene signatures, which was then used to divide patients into the high and low risk groups. We analyzed differences in enriched pathways between the two groups, following which we investigated whether the risk score was correlated with stemness scores, genes related to m6A, m5C, m1A and m7G modification, the immune microenvironment, immunotherapy response, and multiple anti-cancer drug sensitivity.

Results: Overall survival differed significantly between the high-risk and low-risk groups (HR = 4.26). The AUCs for predicting 1-, 3-, and 5-year survival were 0.763, 0.766, and 0.728, respectively. In the GSE68465 dataset, the HR was 2.03, while the AUCs for predicting 1-, 3-, and 5-year survival were 0.69, 0.651, and 0.618, respectively. The corresponding values in the GSE72094 dataset were an HR of 2.36 and AUCs of 0.653, 0.662, and 0.749, respectively. The risk score model could independently predict OS in patients with LUAD, and highly correlated with stemness scores and numerous m6A, m5C, m1A and m7G modification-related genes. Furthermore, the risk model was significantly correlated with multiple immune microenvironment characteristics. In the GSE135222 dataset, the HR was 4.26 and the AUC was 0.702. Evaluation of the GSE126044 and IMvigor210 cohorts indicated that PD-1/PD-LI inhibitor treatment may be indicated in patients with low risk scores, while anti-cancer therapy with various drugs may be indicated in patients with high risk scores.

Conclusion: Our novel risk model developed based on hypoxia, immune, and EMT gene signatures can aid in predicting clinical prognosis and guiding treatment in patients with LUAD.

Single-Cell RNA-Seq Reveals the Promoting Role of Ferroptosis Tendency During Lung Adenocarcinoma EMT Progression

Epithelial-mesenchymal transition (EMT) and ferroptosis are two important processes in biology. In tumor cells, they are intimately linked. We used single-cell RNA sequencing to investigate the regulatory connection between EMT and ferroptosis tendency in LUAD epithelial cells. We used Seurat to construct the expression matrix using the GEO dataset GSE131907 and extract epithelial cells. We found a positive correlation between the trends of EMT and ferroptosis tendency. Then we used SCENIC to analyze differentially activated transcription factors and constructed a molecular regulatory directed network by causal inference. Some ferroptosis markers (GPX4, SCP2, CAV1) were found to have strong regulatory effects on EMT. Cell communication networks were constructed by iTALK and implied that Ferro_High_EMT_High cells have a higher expression of SDC1, SDC4, and activation of LGALS9-HARVCR2 pathways. By deconvolution of bulk sequencing, the results of CIBERSORTx showed that the co-occurrence of ferroptosis tendency and EMT may lead to tumor metastasis and non-response to immunotherapy. Our findings showed there is a strong correlation between ferroptosis tendency and EMT. Ferroptosis may have a promotive effect on EMT. High propensities of ferroptosis and EMT may lead to poor prognosis and non-response to immunotherapy.