Construction of a Prognostic Immune Signature for Squamous-Cell Lung Cancer to Predict Survival

Development of a TGF-β signaling-related genes signature to predict clinical prognosis and immunotherapy responses in clear cell renal cell carcinoma

Background

Transforming growth factor (TGF)-β signaling is strongly related to the development and progression of tumor. We aimed to construct a prognostic gene signature based on TGF-β signaling-related genes for predicting clinical prognosis and immunotherapy responses of patients with clear cell renal cell carcinoma (ccRCC).

Methods

The gene expression profiles and corresponding clinical information of ccRCC were collected from the TCGA and the ArrayExpress (E-MTAB-1980) databases. LASSO, univariate and multivariate Cox regression analyses were conducted to construct a prognostic signature in the TCGA cohort. The E-MTAB-1980 cohort were used for validation. Kaplan-Meier (K-M) survival and time-dependent receiver operating characteristic (ROC) were conducted to assess effectiveness and reliability of the signature. The differences in gene enrichments, immune cell infiltration, and expression of immune checkpoints in ccRCC patients showing different risks were investigated.

Results

We constructed a seven gene (PML, CDKN2B, COL1A2, CHRDL1, HPGD, CGN and TGFBR3) signature, which divided the ccRCC patients into high risk group and low risk group. The K-M analysis indicated that patients in the high risk group had a significantly shorter overall survival (OS) time than that in the low risk group in the TCGA (p < 0.001) and E-MTAB-1980 (p = 0.012). The AUC of the signature reached 0.77 at 1 year, 0.7 at 3 years, and 0.71 at 5 years in the TCGA, respectively, and reached 0.69 at 1 year, 0.72 at 3 years, and 0.75 at 5 years in the E-MTAB-1980, respectively. Further analyses confirmed the risk score as an independent prognostic factor for ccRCC (p < 0.001). The results of ssGSEA that immune cell infiltration degree and the scores of immune-related functions were significantly increased in the high risk group. The CIBERSORT analysis indicated that the abundance of immune cell were significantly different between two risk groups. Furthermore, The risk score was positively related to the expression of PD-1, CTLA4 and LAG3.These results indicated that patients in the high risk group benefit more from immunotherapy.

Conclusion

We constructed a novel TGF-β signaling-related genes signature that could serve as an promising independent factor for predicting clinical prognosis and immunotherapy responses in ccRCC patients.

A Ferroptosis-Related LncRNA Signature Associated with Prognosis, Tumor Immune Environment, and Genome Instability in Hepatocellular Carcinoma

Background

Ferroptosis is an iron-dependent form of cell death. In this study, we identified ferroptosis-related long noncoding RNAs (FRlncRNAs) to investigate their association with hepatocellular carcinoma (HCC) in prognosis, tumor immune environment, and genome instability.

Methods

Transcriptome profile data of HCC were retrieved from a public database. FRlncRNAs were identified by co-expression analysis. Patients were randomly divided into training and test cohorts. Univariate Cox analysis and Least Absolute Shrinkage and Selection Operator (LASSO) Cox regression were performed to construct a risk model. Patients were divided into high- and low-risk groups based on the risk model. AUC and C index were used to assess the risk model. Survival analysis, immune status, and genome instability were compared between the two groups.

Results

Sixteen FRlncRNAs were identified and used to construct an FRlncRNA signature for the risk model. The Kaplan-Meier analysis revealed that patients in the high-risk group had poorer overall survival than patients in the low-risk group. The area under curve of the risk model was 0.879, 0.809, and 0.757 in the training cohort and 0.635, 0.688, and 0.739 in the test cohort at 1, 2, and 3 years, respectively. The risk model was an independent prognostic predictor and showed excellent prediction of prognosis compared with clinicopathological features. For the differentially expressed ferroptosis-related genes, many enriched metabolic pathways were identified in the functional enrichment analysis. Immune cells such as CD8+ T cells, macrophages M1, natural killer cells, and B cells, which may be associated with antitumor immune responses, differed between the high- and low-risk groups. Genome instability based on the risk model was also explored. A total of 61 genes were differently mutated between the two risk groups, and among them, TP53, HECW2, TRIM66, MCTP2, and KIAA1551 had the most significant mutation frequency differences.

Conclusion

The FRlncRNA signature is closely related with overall survival, tumor immune environment, and genome instability in HCC.

Development of a prognostic signature based on immune-related genes and the correlation with immune microenvironment in breast cancer

Breast cancer (BC) is an inflammatory tumor caused by a variety of pathological factors, and is still the most common malignant tumor in women. Immune-related genes (IRGs) play a prominent role in the oncogenesis and progression of BC, and are of tumor-specific expression patterns that would benefit the prognosis evaluation. However, there were no systematic studies concerning the possibilities of IRGs in BC prognosis. In this study, the Cancer Genome Atlas (TCGA) database was used to integrate the expression profiles of IRG with the overall survival (OS) rate of 1039 breast cancer patients. The Cox regression analysis was used to predict the survival-related IRGs in BC. Then, we successfully screened a total of 6 IRGs, including PSME2, ULBP2, IGHE, SCG2, SDC1, and SSTR1, and accordingly constructed a prognosis prediction model of BC. Based on the IRG-related model, the BC patients were divided into high- and low-risk groups, and the association between the prognostic model and tumor immune microenvironment (TME) was further explored. The prognostic model reflected the infiltration of various immune cells. Moreover, the low-risk group was found to be with higher immunophenoscore and distinct mutation signatures compared with the high-risk group. The histological validation showed that SDC1, as well as M2 macrophage biomarker CD206, were both of higher abundance in BC samples of high-risk patients, compared with those of low-risk patients. Our results identify the clinically significant IRGs and demonstrate the importance of the IRG-based immune prognostic model in BC monitoring, prognosis prediction, and therapy.

Comprehensive Analysis of Alteration Landscape and Its Clinical Significance of Mitochondrial Energy Metabolism Pathway-Related Genes in Lung Cancers

Background

Mitochondria are the energy factories of cells. The abnormality of mitochondrial energy metabolism pathways is closely related to the occurrence and development of lung cancer. The abnormal genes in mitochondrial energy metabolism pathways might be the novel targets and biomarkers to diagnose and treat lung cancers.

Method

Genes in major mitochondrial energy metabolism pathways were obtained from the KEGG database. The transcriptomic, mutation, and clinical data of lung cancers were obtained from The Cancer Genome Atlas (TCGA) database. Genes and clinical biomarkers were mined that affected lung cancer survival. Gene enrichment analysis was performed with ClusterProfiler and the gene set enrichment analysis (GSEA). STRING database and Cytoscape were used for protein-protein interaction (PPI) analysis. The diagnostic biomarker pattern of lung cancer was optimized, and its accuracy was verified with 10-fold cross-validation. The four genes screened by logistic regression model were verified by western blot in 5 pairs of lung cancer specimens collected in hospital.

Results

In total, 188 mitochondrial energy metabolism pathway-related genes (MMRGs) were included in this study. GSEA analysis found that MMRGs in the lung cancer group were mainly enriched in the metabolic pathway of oxidative phosphorylation and electron respiratory transport chain compared to the control group. Age did not affect the mutation frequency of MMRGs. Comparative analysis of these 188 MMRGs identified 43 differentially expressed MMRGs (24 upregulated and 19 downregulated) in the lung cancer group compared to the control group. The survival analysis of these 43 differentially expressed MMRGs found that the survival time was better in the low-expressed GAPDHS group than that in the high-expressed GAPDHS group of lung cancers. The advanced age, high expression of GAPDHS, low expressions of ACSBG1 and CYP4A11, and ACOX3 mutation were biomarkers of poor prognosis in lung cancers. PPI analysis showed that proteins such as GAPDH and GAPDHS interacted with many proteins in mitochondrial metabolic pathways. A four-MMRG-signature model (y = 0.0069∗ACADL - 0.001∗ALDH18A1 - 0.0405∗CPT1B + 0.0008∗PPARG - 1.625) was established to diagnose lung cancer with the accuracy up to 98.74%, AUC value up to 0.992, and a missed diagnosis rate of only 0.6%. Western blotting showed that ALDH18A1 and CPT1B proteins were significantly overexpressed in the lung cancer group (p < 0.05), and ACADL and PPARG proteins were slightly underexpressed in the lung cancer group (p < 0.05), which were consistent with the results of their corresponding mRNA expressions.

Conclusion

Mitochondrial energy metabolism pathway alterations are the important hallmarks of lung cancer. Age did not increase the risk of MMRG mutation. High expression of GAPDHS, low expression of ACSBG1, low expression of CYP4A11, mutated ACOX3, and old age predict a poor prognosis of lung cancer. Four differentially expressed MMRGs (ACADL, ALDH18A1, CPT1B, and PPARG) established a logistic regression model, which could effectively diagnose lung cancer. At the protein level, ALDH18A1 and CPT1B were significantly upregulated, and ACADL and PPARG were slightly underexpressed, in the lung cancer group compared to the control group, which were consistent with the results of their corresponding mRNA expressions.

Identification and Validation of Immune-Related Gene Signature for Predicting Lymph Node Metastasis and Prognosis in Lung Adenocarcinoma

Lung cancer is a serious malignancy, and lung adenocarcinoma (LUAD) is the most common pathological subtype. Immune-related factors play an important role in lymph node metastasis. In this study, we obtained gene expression profile data for LUAD and normal tissues from the TCGA database and analyzed their immune-related genes (IRGs), and observed that 459 IRGs were differentially expressed. Further analysis of the correlation between differentially expressed IRGs and lymph node metastasis revealed 18 lymph node metastasis-associated IRGs. In addition, we analyzed the mutations status, function and pathway enrichment of these IRGs, and regulatory networks established through TF genes. We then identified eight IRGs (IKBKB, LTBR, MIF, PPARD, PPIA, PSME3, S100A6, SEMA4B) as the best predictors by LASSO Logistic analysis and used these IRGs to construct a model to predict lymph node metastasis in patients with LUAD (AUC 0.75; 95% CI: 0.7064-0.7978), and survival analysis showed that the risk score independently affected patient survival. We validated the predictive effect of risk scores on lymph node metastasis and survival using the GEO database as a validation cohort and the results showed good agreement. In addition, the risk score was highly correlated with infiltration of immune cells (mast cells activated, macrophages M2, macrophages M0 and B cells naïve), immune and stromal scores, and immune checkpoint genes (LTBR, CD40LG, EDA2R, and TNFRSF19). We identified key IRGs associated with lymph node metastasis in LUAD and constructed a reliable risk score model, which may provide valuable biomarkers for LUAD patients and further reveal the mechanism of its occurrence.

Identification and Validation of Immune-Related Gene Signature for Predicting Lymph Node Metastasis and Prognosis in Lung Adenocarcinoma

Lung cancer is a serious malignancy, and lung adenocarcinoma (LUAD) is the most common pathological subtype. Immune-related factors play an important role in lymph node metastasis. In this study, we obtained gene expression profile data for LUAD and normal tissues from the TCGA database and analyzed their immune-related genes (IRGs), and observed that 459 IRGs were differentially expressed. Further analysis of the correlation between differentially expressed IRGs and lymph node metastasis revealed 18 lymph node metastasis-associated IRGs. In addition, we analyzed the mutations status, function and pathway enrichment of these IRGs, and regulatory networks established through TF genes. We then identified eight IRGs (IKBKB, LTBR, MIF, PPARD, PPIA, PSME3, S100A6, SEMA4B) as the best predictors by LASSO Logistic analysis and used these IRGs to construct a model to predict lymph node metastasis in patients with LUAD (AUC 0.75; 95% CI: 0.7064-0.7978), and survival analysis showed that the risk score independently affected patient survival. We validated the predictive effect of risk scores on lymph node metastasis and survival using the GEO database as a validation cohort and the results showed good agreement. In addition, the risk score was highly correlated with infiltration of immune cells (mast cells activated, macrophages M2, macrophages M0 and B cells naïve), immune and stromal scores, and immune checkpoint genes (LTBR, CD40LG, EDA2R, and TNFRSF19). We identified key IRGs associated with lymph node metastasis in LUAD and constructed a reliable risk score model, which may provide valuable biomarkers for LUAD patients and further reveal the mechanism of its occurrence.