A 5-Gene Signature Is Closely Related to Tumor Immune Microenvironment and Predicts the Prognosis of Patients with Non-Small Cell Lung Cancer

A new immune signature for survival prediction and immune checkpoint molecules in non-small cell lung cancer

Background

Immune checkpoint blockade (ICB) therapy has brought remarkable clinical benefits to patients with advanced non-small cell lung carcinoma (NSCLC). However, the prognosis remains largely variable.

Methods

The profiles of immune-related genes for patients with NSCLC were extracted from TCGA database, ImmPort dataset, and IMGT/GENE-DB database. Coexpression modules were constructed using WGCNA and 4 modules were identified. The hub genes of the module with the highest correlations with tumor samples were identified. Then integrative bioinformatics analyses were performed to unveil the hub genes participating in tumor progression and cancer-associated immunology of NSCLC. Cox regression and Lasso regression analyses were conducted to screen prognostic signature and to develop a risk model.

Results

Functional analysis showed that immune-related hub genes were involved in the migration, activation, response, and cytokine-cytokine receptor interaction of immune cells. Most of the hub genes had a high frequency of gene amplifications. MASP1 and SEMA5A presented the highest mutation rate. The ratio of M2 macrophages and naïve B cells revealed a strong negative association while the ratio of CD8 T cells and activated CD4 memory T cells showed a strong positive association. Resting mast cells predicted superior overall survival. Interactions including protein-protein, lncRNA and transcription factor interactions were analyzed and 9 genes were selected by LASSO regression analysis to construct and verify a prognostic signature. Unsupervised hub genes clustering resulted in 2 distinct NSCLC subgroups. The TIDE score and the drug sensitivity of gemcitabine, cisplatin, docetaxel, erlotinib and paclitaxel were significantly different between the 2 immune-related hub gene subgroups.

Conclusions

These findings suggested that our immune-related genes can provide clinical guidance for the diagnosis and prognosis of different immunophenotypes and facilitate the management of immunotherapy in NSCLC.

Integrating Radiomics with Genomics for Non-Small Cell Lung Cancer Survival Analysis

Purpose

The objectives of our study were to assess the association of radiological imaging and gene expression with patient outcomes in non-small cell lung cancer (NSCLC) and construct a nomogram by combining selected radiomic, genomic, and clinical risk factors to improve the performance of the risk model.

Methods

A total of 116 cases of NSCLC with CT images, gene expression, and clinical factors were studied, wherein 87 patients were used as the training cohort, and 29 patients were used as an independent testing cohort. Handcrafted radiomic features and deep-learning genomic features were extracted and selected from CT images and gene expression analysis, respectively. Two risk scores were calculated through Cox regression models for each patient based on radiomic features and genomic features to predict overall survival (OS). Finally, a fusion survival model was constructed by incorporating these two risk scores and clinical factors.

Results

The fusion model that combined CT images, gene expression data, and clinical factors effectively stratified patients into low- and high-risk groups. The C-indexes for OS prediction were 0.85 and 0.736 in the training and testing cohorts, respectively, which was better than that based on unimodal data.

Conclusions

Combining radiomics and genomics can effectively improve OS prediction for NSCLC patients.

Predicting Differences in Treatment Response and Survival Time of Lung Adenocarcinoma Patients Based on a Prognostic Risk Model of Glycolysis-Related Genes

Background: Multiple factors influence the survival of patients with lung adenocarcinoma (LUAD). Specifically, the therapeutic outcomes of treatments and the probability of recurrence of the disease differ among patients with the same stage of LUAD. Therefore, effective prognostic predictors need to be identified.

Methods: Based on the tumor mutation burden (TMB) data obtained from The Cancer Genome Atlas (TCGA) database, LUAD patients were divided into high and low TMB groups, and differentially expressed glycolysis-related genes between the two groups were screened. The least absolute shrinkage and selection operator (LASSO) and Cox regression were used to obtain a prognostic model. A receiver operating characteristic (ROC) curve and a calibration curve were generated to evaluate the nomogram that was constructed based on clinicopathological characteristics and the risk score. Two data sets (GSE68465 and GSE11969) from the Gene Expression Omnibus (GEO) were used to verify the prognostic performance of the gene. Furthermore, differences in immune cell distribution, immune-related molecules, and drug susceptibility were assessed for their relationship with the risk score.

Results: We constructed a 5-gene signature (FKBP4, HMMR, B4GALT1, SLC2A1, STC1) capable of dividing patients into two risk groups. There was a significant difference in overall survival (OS) times between the high-risk group and the low-risk group (p < 0.001), with the low-risk group having a better survival outcome. Through multivariate Cox analysis, the risk score was confirmed to be an independent prognostic factor (HR = 2.709, 95% CI = 1.981–3.705, p < 0.001), and the ROC curve and nomogram exhibited accurate prediction performance. Validation of the data obtained in the GEO database yielded similar results. Furthermore, there were significant differences in sensitivity to immunotherapy, cisplatin, paclitaxel, gemcitabine, docetaxel, gefitinib, and erlotinib between the low-risk and high-risk groups.

Conclusion: Our results reveal that glycolysis-related genes are feasible predictors of survival and the treatment response of patients with LUAD.

Natural Killer Cells and Dendritic Cells: Expanding Clinical Relevance in the Non-Small Cell Lung Cancer (NSCLC) Tumor Microenvironment

Non-small cell lung cancer (NSCLC) is a major subtype of lung cancer that accounts for almost 85% of lung cancer cases worldwide. Although recent advances in chemotherapy, radiotherapy, and immunotherapy have helped in the clinical management of these patients, the survival rate in advanced stages remains dismal. Furthermore, there is a critical lack of accurate prognostic and stratification markers for emerging immunotherapies. To harness immune response modalities for therapeutic benefits, a detailed understanding of the immune cells in the complex tumor microenvironment (TME) is required. Among the diverse immune cells, natural killer (NK cells) and dendritic cells (DCs) have generated tremendous interest in the scientific community. NK cells play a critical role in tumor immunosurveillance by directly killing malignant cells. DCs link innate and adaptive immune systems by cross-presenting the antigens to T cells. The presence of an immunosuppressive milieu in tumors can lead to inactivation and poor functioning of NK cells and DCs, which results in an adverse outcome for many cancer patients, including those with NSCLC. Recently, clinical intervention using modified NK cells and DCs have shown encouraging response in advanced NSCLC patients. Herein, we will discuss prognostic and predictive aspects of NK cells and DC cells with an emphasis on NSCLC. Additionally, the discussion will extend to potential strategies that seek to enhance the anti-tumor functionality of NK cells and DCs.

Systematic profiling of invasion-related gene signature predicts prognostic features of lung adenocarcinoma

Due to the high heterogeneity of lung adenocarcinoma (LUAD), molecular subtype based on gene expression profiles is of great significance for diagnosis and prognosis prediction in patients with LUAD. Invasion-related genes were obtained from the CancerSEA database, and LUAD expression profiles were downloaded from The Cancer Genome Atlas. The ConsensusClusterPlus was used to obtain molecular subtypes based on invasion-related genes. The limma software package was used to identify differentially expressed genes (DEGs). A multi-gene risk model was constructed by Lasso-Cox analysis. A nomogram was also constructed based on risk scores and meaningful clinical features. 3 subtypes (C1, C2 and C3) based on the expression of 97 invasion-related genes were obtained. C3 had the worst prognosis. A total of 669 DEGs were identified among the subtypes. Pathway enrichment analysis results showed that the DEGs were mainly enriched in the cell cycle, DNA replication, the p53 signalling pathway and other tumour-related pathways. A 5-gene signature (KRT6A, MELTF, IRX5, MS4A1 and CRTAC1) was identified by using Lasso-Cox analysis. The training, validation and external independent cohorts proved that the model was robust and had better prediction ability than other lung cancer models. The gene expression results showed that the expression levels of MS4A1 and KRT6A in tumour tissues were higher than in normal tissues, while CRTAC1 expression in tumour tissues was lower than in normal tissues. The 5-gene signature prognostic stratification system based on invasion-related genes could be used to assess prognostic risk in patients with LUAD.

A bioinformatic analysis of the inhibin-betaglycan-endoglin/CD105 network reveals prognostic value in multiple solid tumors

Inhibins and activins are dimeric ligands belonging to the TGFβ superfamily with emergent roles in cancer. Inhibins contain an α-subunit (INHA) and a β-subunit (either INHBA or INHBB), while activins are mainly homodimers of either βA (INHBA) or βB (INHBB) subunits. Inhibins are biomarkers in a subset of cancers and utilize the coreceptors betaglycan (TGFBR3) and endoglin (ENG) for physiological or pathological outcomes. Given the array of prior reports on inhibin, activin and the coreceptors in cancer, this study aims to provide a comprehensive analysis, assessing their functional prognostic potential in cancer using a bioinformatics approach. We identify cancer cell lines and cancer types most dependent and impacted, which included p53 mutated breast and ovarian cancers and lung adenocarcinomas. Moreover, INHA itself was dependent on TGFBR3 and ENG/CD105 in multiple cancer types. INHA, INHBA, TGFBR3, and ENG also predicted patients' response to anthracycline and taxane therapy in luminal A breast cancers. We also obtained a gene signature model that could accurately classify 96.7% of the cases based on outcomes. Lastly, we cross-compared gene correlations revealing INHA dependency to TGFBR3 or ENG influencing different pathways themselves. These results suggest that inhibins are particularly important in a subset of cancers depending on the coreceptor TGFBR3 and ENG and are of substantial prognostic value, thereby warranting further investigation.

Homology-Based Image Processing for Automatic Classification of Histopathological Images of Lung Tissue

Simple Summary

The purpose of this study was to develop a computer-aided diagnosis (CAD) system for automatic classification of histopathological images of lung tissues. Homology-based image processing (HI) was proposed for CAD. For developing and validating CAD with HI, two datasets of histopathological images of lung tissues were used. The private dataset consists of 94 histopathological images that were obtained for the following five categories: normal, emphysema, atypical adenomatous hyperplasia, lepidic pattern of adenocarcinoma, and invasive adenocarcinoma. The public dataset consists of 15,000 histopathological images that were obtained for the following three categories: lung adenocarcinoma, lung squamous cell carcinoma, and benign lung tissue. For the two datasets, our results show that HI was more useful than conventional texture analysis for the CAD system.

Abstract

The purpose of this study was to develop a computer-aided diagnosis (CAD) system for automatic classification of histopathological images of lung tissues. Two datasets (private and public datasets) were obtained and used for developing and validating CAD. The private dataset consists of 94 histopathological images that were obtained for the following five categories: normal, emphysema, atypical adenomatous hyperplasia, lepidic pattern of adenocarcinoma, and invasive adenocarcinoma. The public dataset consists of 15,000 histopathological images that were obtained for the following three categories: lung adenocarcinoma, lung squamous cell carcinoma, and benign lung tissue. These images were automatically classified using machine learning and two types of image feature extraction: conventional texture analysis (TA) and homology-based image processing (HI). Multiscale analysis was used in the image feature extraction, after which automatic classification was performed using the image features and eight machine learning algorithms. The multicategory accuracy of our CAD system was evaluated in the two datasets. In both the public and private datasets, the CAD system with HI was better than that with TA. It was possible to build an accurate CAD system for lung tissues. HI was more useful for the CAD systems than TA.

Development and validation of an immune-related prognostic signature in lung adenocarcinoma

BACKGROUND

Lung adenocarcinomas (LUAD) is the most common histological subtype of lung cancers. Tumor immune microenvironment (TIME) is involved in tumorigeneses, progressions, and metastases. This study is aimed to develop a robust immune-related signature of LUAD.

METHODS

A total of 1774 LUAD cases sourced from public databases were included in this study. Immune scores were calculated through ESTIMATE algorithm and weighted gene co-expression network analysis (WGCNA) was applied to identify immune-related genes. Stability selections and Lasso COX regressions were implemented to construct prognostic signatures. Validations and comparisons with other immune-related signatures were conducted in independent Gene Expression Omnibus (GEO) cohorts. Abundant infiltrated immune cells and pathway enrichment analyses were carried out, respectively, through ImmuCellAI and gene set enrichment analysis (GSEA).

RESULTS

In Cancer Genome Atlas (TCGA) LUAD cohorts, immune scores of higher levels were significantly associated with better prognoses (P < .05). Yellow (n = 270) and Blue (n = 764) colored genes were selected as immune-related genes, and after univariate Cox regression analysis (P < .005), a total of 133 genes were screened out for subsequent model constructions. A four-gene signature (ARNTL2, ECT2, PPIA, and TUBA4A) named IPSLUAD was developed through stability selection and Lasso COX regression. It was suggested by multivariate and subgroup analyses that IPSLUAD was an independent prognostic factor. It was suggested by Kaplan-Meier survival analysis that eight out of nine patients in high-risk groups had significantly worse prognoses in validation data sets (P < .05). IPSLUAD outperformed other signatures in two independent cohorts.

CONCLUSIONS

A robust immune-related prognostic signature with great performances in multiple LUAD cohorts was developed in this study.