A 14 immune-related gene signature predicts clinical outcomes of kidney renal clear cell carcinoma

A novel immune-related gene signature for predicting immunotherapy outcomes and survival in clear cell renal cell carcinoma

Clear cell renal carcinoma (ccRCC) is one of the most common cancers worldwide. In this study, a new model of immune-related genes was developed to predict the overall survival and immunotherapy efficacy in patients with ccRCC. Immune-related genes were obtained from the ImmPort database. Clinical data and transcriptomics of ccRCC samples were downloaded from GSE29609 and The Cancer Genome Atlas. An immune-related gene-based prognostic model (IRGPM) was developed using the least absolute shrinkage and selection operator regression algorithm and multivariate Cox regression. The reliability of the developed models was evaluated by Kaplan-Meier survival curves and time-dependent receiver operating characteristic curves. Furthermore, we constructed a nomogram based on the IRGPM and multiple clinicopathological factors, along with a calibration curve to examine the predictive power of the nomogram. Overall, this study investigated the association of IRGPM with immunotherapeutic efficacy, immune checkpoints, and immune cell infiltration. Eleven IRGs based on 528 ccRCC samples significantly associated with survival were used to construct the IRGPM. Remarkably, the IRGPM, which consists of 11 hub genes (SAA1, IL4, PLAUR, PLXNB3, ANGPTL3, AMH, KLRC2, NR3C2, KL, CSF2, and SEMA3G), was found to predict the survival of ccRCC patients accurately. The calibration curve revealed that the nomogram developed with the IRGPM showed high predictive performance for the survival probability of ccRCC patients. Moreover, the IRGPM subgroups showed different levels of immune checkpoints and immune cell infiltration in patients with ccRCC. IRGPM might be a promising biomarker of immunotherapeutic responses in patients with ccRCC. Overall, the established IRGPM was valuable for predicting survival, reflecting the immunotherapy response and immune microenvironment in patients with ccRCC.

Construction of a novel immune response prediction signature to predict the efficacy of immune checkpoint inhibitors in clear cell renal cell carcinoma patients

Background

Immune checkpoint inhibitor (ICI) treatment has enhanced survival outcomes in clear cell renal cell carcinoma (ccRCC) patients. Nevertheless, the effectiveness of immunotherapy in ccRCC patients is restricted and we intended to develop and characterize an immune response prediction signature (IRPS) to forecast the efficacy of immunotherapy.

Methods

RNA-seq expression profile and clinicopathologic characteristics of 539 kidney cancer and 72 patients with normal specimens, were downloaded from the Cancer Genome Atlas (TCGA) database, while the Gene Expression Omnibus (GEO) database was used as the validation set, which included 24 ccRCC samples. Utilization of the TCGA data and immune genes databases (ImmPort and the InnateDB), we explored through Weighted Gene Co-expression Network Analysis (WGCNA), along with Least Absolute Shrinkage and Selection Operator method (LASSO), and constructed an IRPS for kidney cancer patients. GSEA and CIBERSORT were performed to declare the molecular and immunologic mechanism underlying the predictive value of IRPS. The Human Protein Atlas (HPA) was deployed to verify the protein expressions of IRPS genes. Tumor immune dysfunction and exclusion (TIDE) score and immunophenoscore (IPS) were computed to determine the risk of immune escape and value the discrimination of IRPS. A ccRCC cohort with anti-PD-1 therapy was obtained as an external validation data set to verify the predictive value of IRPS.

Results

We constructed a 10 gene signature related to the prognosis and immune response of ccRCC patients. Considering the IRPS risk score, patients were split into high and low risk groups. Patients with high risk in the TCGA cohort tended towards advanced tumor stage and grade with poor prognosis (p < 0.001), which was validated in GEO database (p = 0.004). High-risk group tumors were related with lower PD-L1 expression, higher TMB, higher MSIsensor score, lower IPS, higher TIDE score, and enriched Treg cells, which might be the potential mechanism of immune dysfunction and exclusion. Patients in the IRPS low risk group had better PFS (HR:0.73; 95% CI: 0.54-1.0; P = 0.047).

Conclusion

A novel biomarker of IRPS was constructed to predict the benefit of immunotherapy, which might lead to more individualized prognoses and tailored therapy for kidney cancer patients.

A multivariate statistical test for differential expression analysis

Statistical tests of differential expression usually suffer from two problems. Firstly, their statistical power is often limited when applied to small and skewed data sets. Secondly, gene expression data are usually discretized by applying arbitrary criteria to limit the number of false positives. In this work, a new statistical test obtained from a convolution of multivariate hypergeometric distributions, the Hy-test, is proposed to address these issues. Hy-test has been carried out on transcriptomic data from breast and kidney cancer tissues, and it has been compared with other differential expression analysis methods. Hy-test allows implicit discretization of the expression profiles and is more selective in retrieving both differential expressed genes and terms of Gene Ontology. Hy-test can be adopted together with other tests to retrieve information that would remain hidden otherwise, e.g., terms of (1) cell cycle deregulation for breast cancer and (2) “programmed cell death” for kidney cancer.

A Novel Machine Learning 13-Gene Signature: Improving Risk Analysis and Survival Prediction for Clear Cell Renal Cell Carcinoma Patients

Simple Summary

Clear cell renal cell carcinoma is a type of kidney cancer which comprises the majority of all renal cell carcinomas. Many efforts have been made to identify biomarkers which could help healthcare professionals better treat this kind of cancer. With extensive public data available, we conducted a machine learning study to determine a gene signature that could indicate patient survival with high accuracy. Through the min-Redundancy and Max-Relevance algorithm we generated a signature of 13 genes highly correlated with patient outcomes. These findings reveal potential strategies for personalized medicine in the clinical practice.

Abstract

Patients with clear cell renal cell carcinoma (ccRCC) have poor survival outcomes, especially if it has metastasized. It is of paramount importance to identify biomarkers in genomic data that could help predict the aggressiveness of ccRCC and its resistance to drugs. Thus, we conducted a study with the aims of evaluating gene signatures and proposing a novel one with higher predictive power and generalization in comparison to the former signatures. Using ccRCC cohorts of the Cancer Genome Atlas (TCGA-KIRC) and International Cancer Genome Consortium (ICGC-RECA), we evaluated linear survival models of Cox regression with 14 signatures and six methods of feature selection, and performed functional analysis and differential gene expression approaches. In this study, we established a 13-gene signature (AR, AL353637.1, DPP6, FOXJ1, GNB3, HHLA2, IL4, LIMCH1, LINC01732, OTX1, SAA1, SEMA3G, ZIC2) whose expression levels are able to predict distinct outcomes of patients with ccRCC. Moreover, we performed a comparison between our signature and others from the literature. The best-performing gene signature was achieved using the ensemble method Min-Redundancy and Max-Relevance (mRMR). This signature comprises unique features in comparison to the others, such as generalization through different cohorts and being functionally enriched in significant pathways: Urothelial Carcinoma, Chronic Kidney disease, and Transitional cell carcinoma, Nephrolithiasis. From the 13 genes in our signature, eight are known to be correlated with ccRCC patient survival and four are immune-related. Our model showed a performance of 0.82 using the Receiver Operator Characteristic (ROC) Area Under Curve (AUC) metric and it generalized well between the cohorts. Our findings revealed two clusters of genes with high expression (SAA1, OTX1, ZIC2, LINC01732, GNB3 and IL4) and low expression (AL353637.1, AR, HHLA2, LIMCH1, SEMA3G, DPP6, and FOXJ1) which are both correlated with poor prognosis. This signature can potentially be used in clinical practice to support patient treatment care and follow-up.

Bioinformatics analysis and identification of hub genes and immune-related molecular mechanisms in chronic myeloid leukemia

BACKGROUND

Chronic myeloid leukemia (CML) is a malignant hyperplastic tumor of the bone marrow originating from pluripotent hematopoietic stem cells. The advent of tyrosine kinase inhibitors (TKIs) has greatly improved the survival rate of patients with CML. However, TKI-resistance leads to the disease recurrence and progression. This study aimed to identify immune-related genes (IRGs) associated with CML progression.

METHODS

We extracted the gene's expression profiles from the Gene Expression Omnibus (GEO). Bioinformatics analysis was used to determine the differentially expressed IRGs of CML and normal peripheral blood mononuclear cells (PBMCs). Functional enrichment and gene set enrichment analysis (GSEA) were used to explore its potential mechanism. Hub genes were identified using Molecular Complex Detection (MCODE) and the CytoHubba plugin. The hub genes' diagnostic value was evaluated using the receiver operating characteristic (ROC). The relative proportions of infiltrating immune cells in each CML sample were evaluated using CIBERSORT. Quantitative real-time PCR (RT-qPCR) was used to validate the hub gene expression in clinical samples.

RESULTS

A total of 31 differentially expressed IRGs were identified. GO analyses revealed that the modules were typically enriched in the receptor ligand activity, cytokine activity, and endopeptidase activity. KEGG enrichment analysis of IRGs revealed that CML involved Th17 cell differentiation, the NF-kappa B signaling pathway, and cytokine-cytokine receptor interaction. A total of 10 hub genes were selected using the PPI network. GSEA showed that these hub genes were related to the gamma-interferon immune response, inflammatory response, and allograft rejection. ROC curve analysis suggested that six hub genes may be potential biomarkers for CML diagnosis. Further analysis indicated that immune cells were associated with the pathogenesis of CML. The RT-qPCR results showed that proteinase 3 (PRTN3), cathepsin G (CTSG), matrix metalloproteinase 9 (MMP9), resistin (RETN), eosinophil derived neurotoxin (RNase2), eosinophil cationic protein (ECP, RNase3) were significantly elevated in CML patients' PBMCs compared with healthy controls.

CONCLUSION

These results improved our understanding of the functional characteristics and immune-related molecular mechanisms involved in CML progression and provided potential diagnostic biomarkers and therapeutic targets.

Characterization of Molecular Heterogeneity Associated With Tumor Microenvironment in Clear Cell Renal Cell Carcinoma to Aid Immunotherapy

Clear cell renal cell carcinoma (ccRCC) is the most common type of kidney cancer and has strong immunogenicity. A systematically investigation of the tumor microenvironment (TME) in ccRCC could contribute to help clinicians develop personalized treatment and facilitate clinical decision-making. In this study, we analyzed the immune-related subtype of ccRCC on the basis of immune-related gene expression data in The Cancer Genome Atlas (TCGA, N = 512) and E-MTAB-1980 (N = 101) dataset, respectively. As a result, two subtypes (C1 and C2) were identified by performing non-negative matrix factorization clustering. Subtype C1 was characterized by increased advance ccRCC cases and immune-related pathways. A higher immune score, stromal score, TMB value, Tumor Immune Dysfunction and Exclusion (TIDE) prediction score, and immune checkpoint genes expression level were also observed in C1. In addition, the C1 subtype might benefit from chemotherapy and immunotherapy. The patients in subtype C2 had more metabolism-related pathways, higher tumor purity, and a better prognosis. Moreover, some small molecular compounds for the treatment of ccRCC were identified between the two subtypes by using the Connectivity Map (CMap) database. Finally, we constructed and validated an immune-related (IR) score to evaluate immune modification individually. A high IR score corresponded to a favorable prognosis compared to a low IR score, while more advanced tumor stage and grade cases were enriched in the low IR score group. The two IR score groups also showed a distinct divergence among immune status, TME, and chemotherapy. The external validation dataset (E-MTAB-1980) and another immunotherapy cohort (IMvigor 210) demonstrated that patients in the high IR score group had a significantly prolonged survival time and clinical benefits compared to the low IR score group. Together, characterization of molecular heterogeneity and IR signature may help develop new insights into the TME of ccRCC and provide new strategies for personalized treatment.

A united risk model of 11 immune‑related gene pairs and clinical stage for prediction of overall survival in clear cell renal cell carcinoma patients

Clear cell renal cell carcinoma (ccRCC) is the most common subtype of renal cancer. Currently, we lack effective risk models for the prognosis of ccRCC patients. Given the significant role of cancer immunity in ccRCC, we aimed to establish a novel united risk model including clinical stage and immune-related gene pairs (IRGPs) to assess the prognosis. The gene expression profile and clinical data of ccRCC patients from The Cancer Genome Atlas and Arrayexpress were divided into training cohort (n = 381), validation cohort 1 (n = 156), and validation cohort 2 (n = 101). Through univariate Cox regression analysis and Least Absolute Shrinkage and Selection Operator analysis, 11 IRGPs were obtained. After further analysis, it was found that clinical stage could be an independent prognostic factor; hence, we used it to construct a united prognostic model with 11 IRGPs. Based on this model, patients were divided into high-risk and low-risk groups. In Kaplan–Meier analysis, a significant difference was observed in overall survival (OS) among all three cohorts (p < 0.001). The calibration curve revealed that the signature model is in high accordance with the observed values of each data cohort. The 1-year, 3-year, and 5-year receiver operating characteristic curves of each data cohort showed better performance than only IRGP signatures. The results of immune infiltration analysis revealed significantly (p < 0.05) higher abundance of macrophages M0, T follicular helper cells, and other tumor infiltrating cells. In summary, we successfully established a united prognostic risk model, which can effectively assess the OS of ccRCC patients.

A new thinking: extended application of genomic selection to screen multiomics data for development of novel hypoxia-immune biomarkers and target therapy of clear cell renal cell carcinoma

Increasing evidences show the clinical significance of the interaction between hypoxia and immune in clear cell renal cell carcinoma (ccRCC) microenvironment. However, reliable prognostic signatures based on a combination of hypoxia and immune have not been well established. Moreover, many studies have only used RNA-seq profiles to screen the prognosis feature of ccRCC. Presently, there is no comprehensive analysis of multiomics data to mine a better one. Thus, we try and get it. First, t-SNE and ssGSEA analysis were used to establish tumor subtypes related to hypoxia-immune, and we investigated the hypoxia-immune-related differences in three types of genetic or epigenetic characteristics (gene expression profiles, somatic mutation, and DNA methylation) by analyzing the multiomics data from The Cancer Genome Atlas (TCGA) portal. Additionally, a four-step strategy based on lasso regression and Cox regression was used to construct a satisfying prognostic model, with average 1-year, 3-year and 5-year areas under the curve (AUCs) equal to 0.806, 0.776 and 0.837. Comparing it with other nine known prognostic biomarkers and clinical prognostic scoring algorithms, the multiomics-based signature performs better. Then, we verified the gene expression differences in two external databases (ICGC and SYSU cohorts). Next, eight hub genes were singled out and seven hub genes were validated as prognostic genes in SYSU cohort. Furthermore, it was indicated high-risk patients have a better response for immunotherapy in immunophenoscore (IPS) analysis and TIDE algorithm. Meanwhile, estimated by GDSC and cMAP database, the high-risk patients showed sensitive responses to six chemotherapy drugs and six candidate small-molecule drugs. In summary, the signature can accurately predict the prognosis of ccRCC and may shed light on the development of novel hypoxia-immune biomarkers and target therapy of ccRCC.