Analysis of expression differences of immune genes in non-small cell lung cancer based on TCGA and ImmPort data sets and the application of a prognostic model

Abstract

Background: There has been little investigation carried out into the activity of immune-related genes in the prognosis of non-small cell lung cancer (NSCLC). Our study set out to analyze the correlation between the differential expression of immune genes and NSCLC prognosis by screening the differential expression of immune genes. Based on the immune genes identified, we aimed to construct a prognostic risk model and explore some novel molecules which have predictive potential for therapeutic effect and prognosis in lung cancer.Methods: Immune gene transcriptome data and clinical data of NSCLC samples were extracted from TCGA database, and transcription factors in the ImmPort dataset were obtained. The data were divided into two groups: normal tissues and tumor tissues. The expression levels of immune genes were compared using the edgeR algorithm, and then differential expression analysis was performed. The survival analysis was carried out by combining differential immune genes with clinical survival time, so that the immune genes influencing the prognosis of NSCLC could be determined. A risk score was calculated based on the expression levels of the immune genes related to the prognosis of NSCLC and their corresponding coefficients to construct a prognostic risk model. This model was used to calculate patient risk scores and perform clinical correlation analysis. The selected molecules were further verified by clinical samples.Results: By comparing NSCLC tissues with normal tissues, a total of 6,778 differentially expressed genes were found (P<0.05), of which 490 were differential immune-related genes. Survival analysis determined 28 differential immune genes to be associated with prognosis (P<0.05). We calculated the patient risk value based on the immune gene prognosis model. The survival curve was drawn according to the patient risk score and showed that the survival prognosis was significantly different for the high-risk and the low-risk groups (P<0.05). The area under the receiver operating characteristic (ROC) curve (AUC) was 0.723, which represented a relatively high true-positive rate. All of the results proved the reliability of our immune gene risk prognostic model. After drawing the risk curve, S100A16, IGKV4, S100P, ANGPTL4, SEMA4B, andLGR4 were found to be the high-risk immune genes in NSCLC. Clinical correlation analysis of survivalrelated differential immune genes revealed that in patients with lymph node metastasis, ANGPTL4 was positively correlated with T stage, S100a16 and SEMA4B were upregulated, and VIPR1 was downregulated.Further analysis revealed that VIPR1 was decreased in metastatic lung cancer compared to non-metastatic lung cancer. Furthermore, the real-time PCR detection of the clinical samples showed that S100A16 expression in lung cancer was increased, while VIPR1 expression in lung cancer was downregulated, which was consistent with the results of our bioinformatics analysis.Conclusions: Based on big data from the TCGA and ImmPort databases, our study analyzed the

Identification of a three-miRNA signature as a novel prognostic model for papillary renal cell carcinoma

“…Although surgery and some replacement therapies including immunotherapy have been widely applied for pRCC patients, the limited response rate and unsatisfied outcomes motivated us to further explore the more appropriate methods to improve the therapeutic efficiency and achieve more personalized treatment from the perspective of predicting prognosis [ 21 – 23 ]. Because of the significantly individual variation, a class of genetic markers including coding and non-coding genes attracted increasing attentions in the past few years [ 24 – 26 ]. Liu et al indicated lncRNA KTN1-AS1 played the remarkable roles on predicting the poor prognosis of non-small cell lung cancer and facilitate tumor progression through regulating miR-23b/DEPDC1 axis [ 27 ].…”

“…The overlaps represented the numbers of genes predicted by more than one database (a miR-34a-5p; b miR-410-3p; c miR-6720-3p) efficiency and achieve more personalized treatment from the perspective of predicting prognosis[21][22][23]. Because of the significantly individual variation, a class of genetic markers including coding and non-coding genes attracted increasing attentions in the past few years[24][25][26]. Liu et al indicated lncRNA KTN1-AS1 played the remarkable roles on predicting the poor prognosis of non-small cell lung cancer and facilitate tumor progression through regulating miR-23b/DEPDC1 axis[27].…”

Identification and validation of glycolysis-related signatures for prognosis in patients with breast cancer

“…Notably, the prediction of tumor prognosis by single mRNA is limited, so it is of great significance to construct comprehensive predict model. For example, the risk model based on systematic immune-related genes shows a superior prognostic values in non-small cell lung cancer [37] . In our work, we constructed a riskscore model according to glycolysis-related mRNAs in breast cancer.…”

Identification of Robust Immune‐Associated lncRNAs Signature for Immune Checkpoint Blockade and Prognosis in Breast Cancer

“…Brie y, genes were recognized as protein-coding genes or non-coding genes based on their Ensembl IDs or Refseq IDs, and only the long non-coding genes in NetAffx Annotation les were retained. We downloaded the immune gene data from the ImmPort data project, and 2,483 immune-associated genes were gained [20,21]. We employed the Pearson correlation to analyze the correlation between immunerelated genes and lncRNAs.…”