Project description:BackgroundAlthough immune checkpoint inhibitors (ICI) are a promising therapeutic strategy for lung adenocarcinoma (LUAD), individual subgroups that might benefit from them are yet to be identified. As T cell-mediated tumor killing (TTK) is an underlying mechanism of ICI, we identified subtypes based on genes associated with TTK sensitivity and assessed their predictive significance for LUAD immunotherapies.MethodsUsing high-throughput screening techniques, genes regulating the sensitivity of T cell-mediated tumor killing (GSTTK) with differential expression and associations with prognosis were discovered in LUAD. Furthermore, patients with LUAD were divided into subgroups using unsupervised clustering based on GSTTK. Significant differences were observed in the tumor immune microenvironment (TIME), genetic mutation and immunotherapy response across subgroups. Finally, the prognostic significance of a scoring algorithm based on GSTTK was assessed.ResultsA total of 6 out of 641 GSTTK exhibited differential expression in LUAD and were associated with prognosis. Patients were grouped into two categories based on the expression of the six GSTTK, which represented different TTK immune microenvironments in LUAD. Immune cell infiltration, survival difference, somatic mutation, functional enrichment and immunotherapy responses also varied between the two categories. Additionally, a scoring algorithm accurately distinguished overall survival rates across populations.ConclusionsTTK had a crucial influence on the development of the varying TIME. Evaluation of the varied TTK modes of different tumors enhanced our understanding of TIME characteristics, wherein the changes in T cell activity in LUAD are reflected. Thus, this study guides the development of more effective therapeutic methods.
Project description:BackgroundIntra-tumor heterogeneity (ITH) plays a vital role in drug resistance and recurrence of lung cancer. We used a mutant-allele tumor heterogeneity (MATH) algorithm to assess ITH and investigated its association with clinical and molecular features in advanced lung adenocarcinoma.MethodsTissues from 63 patients with advanced lung adenocarcinoma were analyzed by next-generation sequencing (NGS) using a panel targeting 520 cancer-relevant genes. We calculated the MATH values from NGS data and further investigated their correlation with clinical and molecular characteristics.ResultsAmong the 63 patients with advanced lung adenocarcinoma, the median value of MATH was 33.06. Patients with EGFR mutation had higher level of MATH score than those with wild-type EGFR status (P = 0.008). Patients with stage IV disease showed a trend to have a higher MATH score than those with stage III (P = 0.052). MATH was higher in patients with disruptive TP53 mutations than in those with non-disruptive mutations (P = 0.036) or wild-type sequence (P = 0.023), but did not differ between tumors with non-disruptive mutations and wild-type TP53 (P = 0.867). High MATH is associated with mutations in mismatch repair (MMR) pathway (P = 0.026) and base excision repair (BER) pathway (P = 0.008). In addition, MATH was found to have a positive correlation with tumor mutational burden (TMB) (Spearman ρ = 0.354; P = 0.004). In 26 patients harboring EGFR mutation treated with first generation EGFR TKI as single-agent therapy, the objective response rate was higher in the Low-MATH group than in the High-MATH group (75% vs. 21%; P = 0.016) and Low-MATH group showed a significantly longer progression-free survival than High-MATH group (median PFS: 13.7 months vs. 10.1 months; P = 0.024).ConclusionsFor patients with advanced lung adenocarcinoma, MATH may serve as a clinically practical biomarker to assess intratumor heterogeneity.
Project description:BackgroundThe high heterogeneity of tumour and the complexity of tumour microenvironment (TME) greatly impacted the tumour development and the prognosis of cancer in the era of immunotherapy. In this study, we aimed to portray the single cell-characterised landscape of lung adenocarcinoma (LUAD), and develop an integrated signature incorporating both tumour heterogeneity and TME for prognosis stratification.MethodsSingle-cell tagged reverse transcription sequencing (STRT-seq) was performed on tumour tissues and matched normal tissues from 14 patients with LUAD for immune landscape depiction and candidate key genes selection for signature construction. Kaplan-Meier survival analyses and in-vitro cell experiments were conducted to confirm the gene functions. The transcriptomic profile of 1949 patients from 11 independent cohorts including nine public datasets and two in-house cohorts were obtained for validation.FindingsWe selected 11 key genes closely related to cell-to-cell interaction, tumour development, T cell phenotype transformation, and Ma/Mo cell distribution, including HLA-DPB1, FAM83A, ITGB4, OAS1, FHL2, S100P, FSCN1, SFTPD, SPP1, DBH-AS1, CST3, and established an integrated 11-gene signature, stratifying patients to High-Score or Low-Score group for better or worse prognosis. Moreover, the prognostically-predictive potency of the signature was validated by 11 independent cohorts, and the immunotherapeutic predictive potency was also validated by our in-house cohort treated by immunotherapy. Additionally, the in-vitro cell experiments and drug sensitivity prediction further confirmed the gene function and generalizability of this signature across the entire RNA profile spectrum.InterpretationThis single cell-characterised 11-gene signature might offer insights for prognosis stratification and potential guidance for treatment selection.FundingSupport for the study was provided by National key research and development project (2022YFC2505004, 2022YFC2505000 to Z.W. and J.W.), Beijing Natural Science Foundation (7242114 to J.X.), National Natural Science Foundation of China of China (82102886 to J.X., 81871889 and 82072586 to Z.W.), Beijing Nova Program (20220484119 to J.X.), NSFC general program (82272796 to J.W.), NSFC special program (82241229 to J.W.), CAMS Innovation Fund for Medical Sciences (2021-1-I2M-012, 2022-I2M-1-009 to Z.W. and J.W.), Beijing Natural Science Foundation (7212084 to Z.W.), CAMS Key lab of translational research on lung cancer (2018PT31035 to J.W.), Aiyou Foundation (KY201701 to J.W.). Medical Oncology Key Foundation of Cancer Hospital Chinese Academy of Medical Sciences (CICAMS-MOCP2022003 to J.X.).
Project description:BackgroundThe acquisition of invasive tumor cell behavior is considered to be the cornerstone of the metastasis cascade. Thus, genetic markers associated with invasiveness can be stratified according to patient prognosis. In this study, we aimed to identify an invasive genetic trait and study its biological relevance in lung adenocarcinoma.Methods250 TCGA patients with lung adenocarcinoma were used as the training set, and the remaining 250 TCGA patients, 500 ALL TCGA patients, 226 patients with GSE31210, 83 patients with GSE30219, and 127 patients with GSE50081 were used as the verification data sets. Subtype classification of all TCGA lung adenocarcinoma samples was based on invasion-associated genes using the R package ConsensusClusterPlus. Kaplan-Meier curves, LASSO (least absolute contraction and selection operator) method, and univariate and multivariate Cox analysis were used to develop a molecular model for predicting survival.ResultsAs a consequence, two molecular subtypes for LUAD were first identified from all TCGA all data sets which were significant on survival time. C1 subtype with poor prognosis has higher clinical characteristics of malignancy, higher mutation frequency of KRAS and TP53, and a lower expression of immune regulatory molecules. 2463 differentially expressed invasion genes between C1 and C2 subtypes were obtained, including 580 upregulation genes and 1883 downregulation genes. Functional enrichment analysis found that upregulated genes were associated with the development of tumor pathways, while downregulated genes were more associated with immunity. Furthermore, 5-invasion gene signature was constructed based on 2463 genes, which was validated in four data sets. This signature divided patients into high-risk and low-risk groups, and the LUDA survival rate of the high-risk group is significantly lower than that of the low-risk group. Multivariate Cox analysis revealed that this gene signature was an independent prognostic factor for LUDA. Compared with other existing models, our model has a higher AUC.ConclusionIn this study, two subtypes were identified. In addition, we developed a 5-gene signature prognostic risk model, which has a good AUC in the training set and independent validation set and is a model with independent clinical characteristics. Therefore, we recommend using this classifier as a molecular diagnostic test to assess the prognostic risk of patients with LUDA.
Project description:Background: Lung adenocarcinoma (LUAD) is a life-threatening malignant tumor, contributing for the largest cancer burden worldwide. Tumor microenvironment (TME) is composed of various immune cells, stromal cells and tumor cells, which is highly associated with the cancer prognosis and the response to immunotherapy, in which macrophages in TME have been revealing a potential target for cancer treatment. In this study, we sought to further explore the role of macrophages in LUAD progression and establish a risk model related to macrophages for LUAD. Methods: We explored immune-related pathways that might be affected by counting positively associated genes in macrophages. Molecular typing was also constructed by mining macrophage-associated genes with prognostic value through COX regression and other analyses. RiskScore prognostic models were constructed using lasso regression and stepwise multifactorial regression analysis. The differences on clinical characteristics among three subtypes (C1, C2, and C3) and RiskScore subtypes were analyzed in TCGA dataset. Immunological algorithms such as TIMER, ssGSEA, MCP-Counter, ESTIMATE, and TIDE were used to calculate the level of difference in immune infiltration between the different subtypes. The TCGA mutation dataset processed by mutect2 was used to demonstrate the frequency of mutations between different molecular subtypes. Finally, nomograms, calibration curves, and decision curves were created to assess the predictive accuracy and reliability of the model. Results: The C1 subtype demonstrated the best prognostic outcome, accompanied by higher levels of immune infiltration and lower mutation frequency, while the majority of patients in the C1 subtype were women under 65 years of age. Myeloid-derived suppressor cell (MDSC) scores were higher in the C3 subtype, suggesting a more severe immune escape, which may have contributed to the tumor evading the immune system resulting in a poorer prognosis for patients. In addition, our RiskScore prognostic model had good predictive accuracy and reliability. Conclusion: This paper provides a study of macrophage-related pathways, immunosuppression, and their mechanisms of action in lung cancer, along with targets for future treatment to guide the optimal treatment of lung cancer.
Project description:Lung adenocarcinoma (LUAD) exhibits high heterogeneity and is well known for its high genetic variation. Recently, the understanding of non-genetic variation provides a new perspective to study the heterogeneity of LUAD. Little is known about whether super-enhancers (SEs) may be primarily responsible for the inter-tumor heterogeneity of LUAD. We used super-enhancer RNA (seRNA) levels of a large-scale clinical well-annotated LUAD cohort to stratify patients into three clusters with different prognosis and other malignant characteristics. Mechanistically, estrogen-related receptor alpha (ERRα) in cluster 3-like cell lines acts as a cofactor of BRD4 to assist SE-promoter loops to activate glycolysis-related target gene expression, thereby promoting glycolysis and malignant progression, which confers a therapeutic vulnerability to glycolytic inhibitors. Our study identified three groups of patients according to seRNA levels, among which patients in cluster 3 have the worst prognosis and vulnerability of glycolysis dependency. We also proposed a 3-TF index model to stratify patients with glycolysis-addicted tumors according to tumor SE stratification.
Project description:BackgroundAlthough checkpoint blockade is a promising approach for the treatment of hepatocellular carcinoma (HCC), subsets of patients expected to show a response have not been established. As T cell-mediated tumor killing (TTK) is the fundamental principle of immune checkpoint inhibitor therapy, we established subtypes based on genes related to the sensitivity to TKK and evaluated their prognostic value for HCC immunotherapies.MethodsGenes regulating the sensitivity of tumor cells to T cell-mediated killing (referred to as GSTTKs) showing differential expression in HCC and correlations with prognosis were identified by high-throughput screening assays. Unsupervised clustering was applied to classify patients with HCC into subtypes based on the GSTTKs. The tumor microenvironment, metabolic properties, and genetic variation were compared among the subgroups. A scoring algorithm based on the prognostic GSTTKs, referred to as the TCscore, was developed, and its clinical and predictive value for the response to immunotherapy were evaluated.ResultsIn total, 18 out of 641 GSTTKs simultaneously showed differential expression in HCC and were correlated with prognosis. Based on the 18 GSTTKs, patients were clustered into two subgroups, which reflected distinct TTK patterns in HCC. Tumor-infiltrating immune cells, immune-related gene expression, glycolipid metabolism, somatic mutations, and signaling pathways differed between the two subgroups. The TCscore effectively distinguished between populations with different responses to chemotherapeutics or immunotherapy and overall survival.ConclusionsTTK patterns played a nonnegligible role in formation of TME diversity and metabolic complexity. Evaluating the TTK patterns of individual tumor will contribute to enhancing our cognition of TME characterization, reflects differences in the functionality of T cells in HCC and guiding more effective therapy strategies.
Project description:Recently, several molecular subtypes with different prognosis have been found in lung adenocarcinoma (LUAD). However, the characteristics of the ferroptosis molecular subtypes and the associated tumor microenvironment (TME) cell infiltration have not been fully studied in LUAD. Using 1160 lung adenocarcinoma samples, we explored the molecular subtypes mediated by ferroptosis-related genes, along with the associated TME cell infiltration. The ferroptosis score was constructed using the least absolute shrinkage and selection operator regression (LASSO) method to quantify the ferroptosis characteristics of a single tumor. Three different molecular subtypes related to ferroptosis, with different prognoses, were identified in LUAD. Analysis of TME cell infiltration revealed immune heterogeneity among the three subtypes. Cluster A was characterized by immunosuppression and was associated with stromal activation. Cluster C was characterized by a large number of immune cells infiltrating the TME, promoting tumor immune response, and it was significantly enriched in immune activation-related signaling pathways. Relatively less infiltration of immune cells was a feature of cluster B. The ferroptosis score can predict tumor subtype, immunity and prognosis. A low ferroptosis score was characterized by immune activation and good prognosis, as seen in the cluster C subtype. Relative immunosuppression and poor prognosis were the characteristics of a high ferroptosis score, as seen in cluster A and B subtypes. At the same time, the anti-PD-1/L1 immunotherapy cohort demonstrated that a low ferroptosis score was associated with higher efficacy of immunotherapy. The ferroptosis score is a promising biomarker that could be of great significance to determine the prognosis, molecular subtypes, TME cell infiltration characteristics and immunotherapy effects in patients with LUAD.
Project description:Intra-tumor heterogeneity (ITH) of human tumors is important for tumor progression, treatment response, and drug resistance. However, the spatial distribution of ITH remains incompletely understood. Here, we present spatial analysis of ITH in lung adenocarcinomas from 147 patients using multi-region mass spectrometry of >5,000 regions, single-cell copy number sequencing of ~2,000 single cells, and cyclic immunofluorescence of >10 million cells. We identified two distinct spatial patterns among tumors, termed clustered and random geographic diversification (GD). These patterns were observed in the same samples using both proteomic and genomic data. The random proteomic GD pattern, which is characterized by decreased cell adhesion and lower levels of tumor-interacting endothelial cells, was significantly associated with increased risk of recurrence or death in two independent patient cohorts. Our study presents comprehensive spatial mapping of ITH in lung adenocarcinoma and provides insights into the mechanisms and clinical consequences of GD.
Project description:Lung cancer is one of the main cancer types due to its persistently high incidence and mortality, yet a simple and effective prognostic model is still lacking. This study aimed to identify independent prognostic genes related to the heterogeneity of lung adenocarcinoma (LUAD), generate a prognostic risk score model, and construct a nomogram in combination with other pathological characteristics to predict patients' overall survival (OS). A significant amount of data pertaining to single-cell RNA sequencing (scRNA-seq), RNA sequencing (RNA-seq), and somatic mutation were used for data mining. After statistical analyses, a risk scoring model was established based on eight independent prognostic genes, and the OS of high-risk patients was significantly lower than that of low-risk patients. Interestingly, high-risk patients were more sensitive and effective to immune checkpoint blocking therapy. In addition, it was noteworthy that CCL20 not only affected prognosis and differentiation of LUAD but also led to poor histologic grade of tumor cells. Ultimately, combining risk score, clinicopathological information, and CCL20 mutation status, a nomogram with good predictive performance and high accuracy was established. In short, our research established a prognostic model that could be used to guide clinical practice based on the constantly updated big multi-omics data. Finally, this analysis revealed that CCL20 may become a potential therapeutic target for LUAD.