Construction of immune-related risk signature for renal papillary cell carcinoma

Prediction of clear cell renal cell carcinoma prognosis based on an immunogenomic landscape analysis

Immune cell infiltration and tumor-related immune molecules play key roles in tumorigenesis and tumor progression. The influence of immune interactions on the molecular characteristics and prognosis of clear cell renal cell carcinoma (ccRCC) remains unclear. A machine learning algorithm was applied to the transcriptome data from The Cancer Genome Atlas database to determine the immunophenotypic and immunological characteristics of ccRCC patients. These algorithms included single-sample gene set enrichment analyses and cell type identification. Using bioinformatics techniques, we examined the prognostic potential and regulatory networks of immune-related genes (IRGs) involved in ccRCC immune interactions. Fifteen IRGs (CCL7, CHGA, CMA1, CRABP2, IFNE, ISG15, NPR3, PDIA2, PGLYRP2, PLA2G2A, SAA1, TEK, TGFA, TNFSF14, and UCN2) were identified as prognostic IRGs associated with overall survival and were used to construct a prognostic model. The area under the receiver operating characteristic curve at 1 year was 0.927; 3 years, 0.822; and 5 years, 0.717, indicating good predictive accuracy. Molecular regulatory networks were found to govern immune interactions in ccRCC. Additionally, we developed a nomogram containing the model and clinical characteristics with high prognostic potential. By systematically examining the sophisticated regulatory mechanisms, molecular characteristics, and prognostic potential of ccRCC immune interactions, we provided an important framework for understanding the molecular mechanisms of ccRCC and identifying new prognostic markers and therapeutic targets for future research.

Identification of prognostic stemness-related genes in kidney renal papillary cell carcinoma

BACKGROUND

Kidney renal papillary cell carcinoma (KIRP) is the second most prevalent malignant cancer originating from the renal epithelium. Nowadays, cancer stem cells and stemness-related genes (SRGs) are revealed to play important roles in the carcinogenesis and metastasis of various tumors. Consequently, we aim to investigate the underlying mechanisms of SRGs in KIRP.

METHODS

RNA-seq profiles of 141 KIRP samples were downloaded from the TCGA database, based on which we calculated the mRNA expression-based stemness index (mRNAsi). Next, we selected the differentially expressed genes (DEGs) between low- and high-mRNAsi groups. Then, we utilized weighted gene correlation network analysis (WGCNA) and univariate Cox analysis to identify prognostic SRGs. Afterwards, SRGs were included in the multivariate Cox regression analysis to establish a prognostic model. In addition, a regulatory network was constructed by Pearson correlation analysis, incorporating key genes, upstream transcription factors (TFs), and downstream signaling pathways. Finally, we used Connectivity map analysis to identify the potential inhibitors.

RESULTS

In total, 1124 genes were characterized as DEGs between low- and high-RNAsi groups. Based on six prognostic SRGs (CCKBR, GPR50, GDNF, SPOCK3, KC877982.1, and MYO15A), a prediction model was established with an area under curve of 0.861. Furthermore, among the TFs, genes, and signaling pathways that had significant correlations, the CBX2-ASPH-Notch signaling pathway was the most significantly correlated. Finally, resveratrol might be a potential inhibitor for KIRP.

CONCLUSIONS

We suggested that CBX2 could regulate ASPH through activation of the Notch signaling pathway, which might be correlated with the carcinogenesis, development, and unfavorable prognosis of KIRP.

Based on the prognosis model of immunogenes, the prognosis model was constructed to predict the invasion of immune genes and immune cells related to primary liver cancer and its experimental validation

Background

Primary liver cancer (PLC) is a prevalent malignancy of the digestive system characterized by insidious symptom onset and a generally poor prognosis. Recent studies have highlighted a significant correlation between the initiation and prognosis of liver cancer and the immune function of PLC patients.

Purpose

Revealing the expression of PLC-related immune genes and the characteristics of immune cell infiltration provides assistance for the analysis of clinical pathological parameters and prognosis of PLC patients.

Methods

PLC-related differentially expressed genes (DEGs) with a median absolute deviation (MAD > 0.5) were identified from TCGA and GEO databases. These DEGs were intersected with immune-related genes (IRGs) from the ImmPort database to obtain PLC-related IRGs. The method of constructing a prognostic model through immune-related gene pairs (IRGPs) is used to obtain IRGPs and conduct the selection of central immune genes. The central immune genes obtained from the selection of IRGPs are validated in PLC. Subsequently, the relative proportions of 22 types of immune cells in different immune risk groups are evaluated, and the differential characteristics of PLC-related immune cells are verified through animal experiments.

Results

Through database screening and the construction of an IRGP prognosis model, 84 pairs of IRGPs (P < 0.001) were ultimately obtained. Analysis of these 84 IRGPs revealed 11 central immune genes related to PLC, showing differential expression in liver cancer tissues compared to normal liver tissues. Results from the CiberSort platform indicate differential expression of immune cells such as naive B cells, macrophages, and neutrophils in different immune risk groups. Animal experiments demonstrated altered immune cell proportions in H22 tumor-bearing mice, validating findings from peripheral blood and spleen homogenate analyses.

Conclusion

Our study successfully predicted and validated PLC-related IRGs and immune cells, suggesting their potential as prognostic indicators and therapeutic targets for PLC.

Network diffusion-based approach for survival prediction and identification of biomarkers using multi-omics data of papillary renal cell carcinoma

Identification of cancer subtypes based on molecular knowledge is crucial for improving the patient diagnosis, prognosis, and treatment. In this work, we integrated copy number variations (CNVs) and transcriptomic data of Kidney Papillary Renal Cell Carcinoma (KIRP) using a network diffusion strategy to stratify cancers into clinically and biologically relevant subtypes. We constructed GeneNet, a KIRP specific gene expression network from RNA-seq data. The copy number variation data was projected onto GeneNet and propagated on the network for clustering. We identified robust subtypes that are biologically informative and significantly associated with patient survival, tumor stage and clinical subtypes of KIRP. We performed a Singular Value Decomposition (SVD) analysis of KIRP subtypes, which revealed the genes/silent players related to poor survival. A differential gene expression analysis between subtypes showed that genes related to immune, extracellular matrix organization, and genomic instability are upregulated in the poor survival group. Overall, the network-based approach revealed the molecular subtypes of KIRP and captured the relationship between gene expression and CNVs. This framework can be further expanded to integrate other omics data.

Seven Immune-Related Genes' Prognostic Value and Correlation with Treatment Outcome in Head and Neck Squamous Cell Carcinoma

Background

Head and neck squamous cell carcinoma (HNSCC) is a growing concern worldwide, due to its poor prognosis, low responsiveness to treatment, and drug resistance. Since immunotherapy effectively improves HNSCC patients' survival status, it is important to continuously explore new immune-related predictive factors to accurately predict the immune landscape and clinical outcomes of individuals suffering from HNSCC.

Methods

The HNSCC transcriptome profiling of RNA-sequencing data was retrieved from TCGA database, and the microarray of GSE27020 was obtained from the GEO database for validation. The differentially expressed genes (DEGs) between HNSCC and normal samples were identified by multiple test corrections in TCGA database. The univariate and multivariate Cox analyses were performed to identify proper immune-related genes (IRGs) to construct a risk model. The Cox regression coefficient was employed for calculation of the risk score (RS) of IRG signature. The median value of RS was utilized as a basis to classify individuals with HNSCC into high- and low-risk groups. The Kaplan-Meier (K-M) survival analysis and receiver operating characteristic (ROC) curves were employed for the identification of the prognostic significance and precision of the IRG signature. The signature was also evaluated based on clinical variables, predictive nomogram, mutation analysis, infiltrating immune cells, immune-related pathways, and chemotherapeutic efficacy. The protein-protein interaction (PPI) network and functional enrichment pathway investigations were utilized to explore possible potential molecular mechanisms. Finally, the hub gene's differential mRNA expression levels were evaluated by means of the Gene Expression Profiling Interactive Analysis (GEPIA), and the Human Protein Atlas (HPA) was utilized for the validation of their translational levels.

Results

Collectively, 1593 DEGs between HNSCC and normal samples were identified, of which 136 IRGs were differentially expressed. Then, the 136 immune-related DEGs were mostly enriched in the cytokine-related signaling pathways by GO and KEGG analyses. After that, a valuable signature based on seven genes (DKK1, GAST, IGHM, IL12RB2, SLURP1, STC2, and TNFRSF4) was designed. The HNSCC patients into the low-risk group and the high-risk group were divided by using the median RS; the HNSCC patients in the high-risk group had a worse survival than those in the low-risk group. The risk signature was verified to be an independent predictive marker for HNSCC patients. Meanwhile, the RS had the largest contribution to survival of these patients based on the predictive nomogram. In addition, the low-risk HNSCC patients exhibited significantly enriched immune cells, along with an association with high chemosensitivity.

Conclusion

The constructed gene signature can independently function as a predictive indicator for the clinical features of HNSCC patients. The low-risk HNSCC subjects might benefit from immunotherapy and chemotherapy.

Phosphogluconate dehydrogenase is a predictive biomarker for immunotherapy in hepatocellular carcinoma

Background

Phosphogluconate dehydrogenase (PGD) is involved in the regulation of various tumors. However, its role in hepatocellular carcinoma (HCC) is poorly understood. This study tried to determine the prognostic efficacy of PGD and its value for immunotherapy in HCC.

Methods

The data from the TCGA database was used to explore the predictive power of PGD expression and methylation on the overall survival (OS) of HCC through Cox regression and the Kaplan-Meier analysis. Then, we used the GEO and ICGC database to further verify the predictive power. Finally, the relationship between PGD and immune cells and the relationship between PGD and the efficacy of immunotherapy were explored through bioinformatics analysis in HCC.

Results

PGD is highly expressed in HCC tissues, which is negatively regulated by PGD methylation. Low PGD expression and PGD hypermethylation predict better OS in HCC patients. Besides, a meta-analysis based on the TCGA, GSE14520, and ICGC databases further confirms that low PGD expression is closely related to favorable OS. Then, we find significant differences of immune cell infiltrations between high and low PGD expression groups. Expressions of immune checkpoints, most HLA members and tumor mutation burden (TMB) are higher in the high PGD expression group, which indicates beneficial efficacy of immunotherapy in this group. And the potential mechanisms of PGD are exhibited.

Conclusion

PGD is an independent prognostic factor of HCC patients and plays an important role in immune cell infiltration and immunotherapy, which indicates that PGD can be used as a predictive biomarker for HCC immunotherapy.

Immune-related 3-lncRNA signature with prognostic connotation in a multi-cancer setting

Background

Advances in our understanding of the tumor microenvironment have radically changed the cancer field, highlighting the emerging need for biomarkers of an active, favorable tumor immune phenotype to aid treatment stratification and clinical prognostication. Numerous immune-related gene signatures have been defined; however, their prognostic value is often limited to one or few cancer types. Moreover, the area of non-coding RNA as biomarkers remains largely unexplored although their number and biological roles are rapidly expanding.

Methods

We developed a multi-step process to identify immune-related long non-coding RNA signatures with prognostic connotation in multiple TCGA solid cancer datasets.

Results

Using the breast cancer dataset as a discovery cohort we found 2988 differentially expressed lncRNAs between immune favorable and unfavorable tumors, as defined by the immunologic constant of rejection (ICR) gene signature. Mapping of the lncRNAs to a coding-non-coding network identified 127 proxy protein-coding genes that are enriched in immune-related diseases and functions. Next, we defined two distinct 20-lncRNA prognostic signatures that show a stronger effect on overall survival than the ICR signature in multiple solid cancers. Furthermore, we found a 3 lncRNA signature that demonstrated prognostic significance across 5 solid cancer types with a stronger association with clinical outcome than ICR. Moreover, this 3 lncRNA signature showed additional prognostic significance in uterine corpus endometrial carcinoma and cervical squamous cell carcinoma and endocervical adenocarcinoma as compared to ICR.

Conclusion

We identified an immune-related 3-lncRNA signature with prognostic connotation in multiple solid cancer types which performed equally well and in some cases better than the 20-gene ICR signature, indicating that it could be used as a minimal informative signature for clinical implementation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12967-022-03654-7.

A prognostic Risk Score model for oral squamous cell carcinoma constructed by 6 glycolysis-immune-related genes

Background

Oral squamous cell carcinoma (OSCC) is the most frequent tumor of the head and neck. The glycolysis-related genes and immune-related genes have been proven prognostic values in various cancers. Our study aimed to test the prognostic value of glycolysis-immune-related genes in OSCC.

Methods

Data of OSCC patients were obtained from the Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases. Enrichment analysis was applied to the glycolysis- and immune-related genes screened by differential expression analysis. Univariate Cox and LASSO Cox analyses were used to filtrate the genes related to the prognosis of OSCC and to construct Risk Score model.

Results

A Risk Score model was constructed by six glycolysis-immune-related genes (including ALDOC, VEGFA, HRG, PADI3, IGSF11 and MIPOL1). High risk OSCC patients (Risk Score >−0.3075) had significantly worse overall survival than that of low risk patients (Risk Score <−0.3075).

Conclusions

The Risk Score model constructed basing on 6 glycolysis-immune-related genes was reliable in stratifying OSCC patients with different prognosis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12903-022-02358-0.

Upregulation of Centromere Proteins as Potential Biomarkers for Esophageal Squamous Cell Carcinoma Diagnosis and Prognosis

Esophageal squamous cell carcinoma (ESCC) has a high incidence and low survival rate, necessitating the identification of novel specific biomarkers. Centromere-associated proteins (CENPs) have been reported to be biomarkers for many cancers, but their roles in ESCC have seldom been investigated. Here, the potential clinical roles of CENPs in ESCC patients were demonstrated by a systematic bioinformatics analysis. Most CENP-encoding genes were differentially expressed between tumor and normal tissues. CENPA, CENPE, CENPF, CENPI, CENPM, CENPN, CENPQ, and CENPR were upregulated universally in the three datasets. Survival analysis demonstrated that high expression of CENPE and CENPQ was positively correlated with the outcomes of ESCC patients. The CENPE-based forecast model was more accurate than the tumor-node-metastasis (TNM) staging-based model, which was classified as stage I/II vs. III/IV. More importantly, the forecast model based on the commonly upregulated CENPs exhibited a much higher area under the curve (AUC) value (0.855) than the currently known TTL, ZNF750, AC016205.1, and BOLA3 biomarkers. The nomogram model integrating the CENPs, TNM stage, and sex was highly accurate in the prognosis of ESCC patients ( $\text{AUC}=0.906$ ). Besides, gene set enrichment analysis (GSEA) demonstrated that CENPE expression is significantly correlated with cell cycle, G2/M checkpoint, mitotic spindle, p53, etc. Finally, in validation experiments, we also found that CENPE and CENPQ were significantly overexpressed in esophageal cancer cells. Taken together, these results clearly suggest that CENPs are clinically promising diagnostic and prognostic biomarkers for ESCC patients.

Functions, Roles, and Biological Processes of Ferroptosis-Related Genes in Renal Cancer: A Pan-Renal Cancer Analysis

Ferroptosis is a cell death process discovered in recent years, highly related to cancer, acute kidney injury, and other diseases. In this study, a pan-renal cancer analysis of ferroptosis-associated genes in renal cancer was performed to construct a multigene joint signature for predicting prognosis in renal cancer patients. First, gene expression profiles were downloaded from the TCGA and GTEx databases to search for genes significantly associated with renal cancer prognosis through differential gene expression analysis, weighted gene co-expression network analysis (WGCNA), and survival analysis. Thereafter, the gene-set enrichment analysis (GSEA) was used to identify the biological processes in which ferroptosis-associated genes might be involved. Weighted gene co-expression network analysis resulted in 4,434 differentially expressed genes (DEGs) and 42 co-expression modules, among which ferroptosis-related genes were distributed in 11 gene modules. The survival analysis screening resulted in three DEGs associated with renal cancer prognosis, namely SLC7A11, HMOX1, and MT1G. Specifically, SLC7A11 and HMOX1 were upregulated in renal cancer tissues, while MT1G was downregulated. Receiver operating characteristic (ROC) curves, combined with Kaplan–Meier and Cox regression analysis, revealed that high expression of SLC7A11 was a prognostic risk factor for four different renal cancers, that low expression of HMOX1 was a poor prognostic marker for patients, and that increased expression of MT1G increased the prognostic risk for three additional classes of renal cancer patients, except for renal papillary cell carcinoma. The GSEA results showed that the ferroptosis-related genes from these screens were mainly associated with signaling pathways related to tumor progression and tumor immunity. This study provides potential biological markers for prognosis prediction in renal cancer patients with different subtypes, and these results imply that ferroptosis is highly associated with renal carcinogenesis progression.

The Identification of a Tumor Infiltration CD8+ T-Cell Gene Signature That Can Potentially Improve the Prognosis and Prediction of Immunization Responses in Papillary Renal Cell Carcinoma

Background

CD8+ T cells, vital effectors pertaining to adaptive immunity, display close relationships to the immunization responses to kill tumor cells. Understanding the effect exerted by tumor infiltration CD8+ T cells in papillary renal cell carcinoma (papRCC) is critical for assessing the prognosis process and responses to immunization therapy in cases with this disease.

Materials and Approaches

The single-cell transcriptome data of papRCC were used for screening CD8+ T-cell-correlated differentially expressed genes to achieve the following investigations. On that basis, a prognosis gene signature associated with tumor infiltration CD8+ T cell was built and verified with The Cancer Genome Atlas data set. Risk scores were determined for papRCC cases and categorized as high- or low-risk groups. The prognosis significance for risk scores was assessed with multiple-variate Cox investigation and Kaplan–Meier survival curves. In addition, the possible capability exhibited by the genetic profiles of cases to assess the response to immunization therapy was further explored.

Results

Six hundred twenty-one cell death-inhibiting RNA genes were screened using single-cell RNA sequencing. A gene signature consisting of seven genes (LYAR, YBX1, PNRC1, TCF25, MYL12B, MINOS1, and LINC01420) was then identified, and this collective was considered to be an independent prognosis indicator that could strongly assess overall survival in papRCC. In addition, the data allowed papRCC cases to fall to cohorts at high and low risks, exhibiting a wide range of clinically related features as well as different CD8+ T-cell immunization infiltration and immunization therapy responses.

Conclusions

Our work provides a possible explanation for the limited response of current immunization checkpoint-inhibiting elements for combating papRCC. Furthermore, the researchers built a novel genetic signature that was able to assess the prognosis and immunotherapeutic response of cases. This may also be considered as a promising therapeutic target for the disease.

SKA3 Serves as a Biomarker for Poor Prognosis in Kidney Renal Papillary Cell Carcinoma

Background

There is a surprising paucity of studies investigating the potential mechanism of SKA3 in the progression and prognosis of kidney renal papillary cell carcinoma (KIRP).

Methods

We used TCGA and other databases to analyze the expression, clinical value, and potential mechanisms of SKA3 in KIRP patients. We also explored therapeutic agents for KIRP through GSCALite.

Results

SKA3 mRNA expression was significantly upregulated and the area under the curve was 0.792 (95% CI 0.727–0.856). Increased SKA3 expression was related to shorter overall survival, disease-specific survival and progression-free survival. Hub genes in protein–protein interactions were CDK1, CDC20, CCNB1, CCNA2, BUB1, AURKB, BUB1B, PLK1, CCNB2, and MAD2L1, which were differentially expressed and also associated with KIRP prognosis. Gene-set enrichment analysis indicated that E2F targets, epithelial–mesenchymal transition, glycolysis, the WNT signaling pathway, and other pathways were highly enriched upon SKA3 upregulation. Gene-set variation analysis of SKA3 and its ten hub genes showed that the significant correlation of cancer-related pathways included the cell cycle, DNA damage, hormone androgen receptor, hormone estrogen receptor, PI3K/Akt, and Ras/MAPK. In addition, we found that MEK inhibitors, ie, trametinib, selumetinib, PD0325901, and RDEA119, may be feasible targeting agents for KIRP patients.

Conclusion

SKA3 might contribute to poor prognosis of KIRP through cell cycle, DNA damage, hormone androgen receptor, hormone estrogen receptor, PI3K/Akt, and RAS/MAPK. SKA3 potentially serves as a prognostic biomarker and target for KIRP.

A novel prognostic cancer-related lncRNA signature in papillary renal cell carcinoma

Background

Papillary renal cell carcinoma (pRCC) ranks second in renal cell carcinoma and the prognosis of pRCC remains poor. Here, we aimed to screen and identify a novel prognostic cancer-related lncRNA signature in pRCC.

Methods

The RNA-seq profile and clinical feature of pRCC cases were downloaded from TCGA database. Significant cancer-related lncRNAs were obtained from the Immlnc database. Differentially expressed cancer-related lncRNAs (DECRLs) in pRCC were screened for further analysis. Cox regression report was implemented to identify prognostic cancer-related lncRNAs and establish a prognostic risk model, and ROC curve analysis was used to evaluate its precision. The correlation between RP11-63A11.1 and clinical characteristics was further analyzed. Finally, the expression level and role of RP11-63A11.1 were studied in vitro.

Results

A total of 367 DECRLs were finally screened and 26 prognostic cancer-related lncRNAs were identified. Among them, ten lncRNAs (RP11-573D15.8, LINC01317, RNF144A-AS1, TFAP2A-AS1, LINC00702, GAS6-AS1, RP11-400K9.4, LUCAT1, RP11-63A11.1, and RP11-156L14.1) were independently associated with prognosis of pRCC. These ten lncRNAs were incorporated into a prognostic risk model. In accordance with the median value of the riskscore, pRCC cases were separated into high and low risk groups. Survival analysis indicated that there was a significant difference on overall survival (OS) rate between the two groups. The area under curve (AUC) in different years indicated that the model was of high efficiency in prognosis prediction. RP11-63A11.1 was mainly expressed in renal tissues and it correlated with the tumor stage, T, M, N classifications, OS, PFS, and DSS of pRCC patients. Consistent with the expression in pRCC tissue samples, RP11-63A11.1 was also down-regulated in pRCC cells. More importantly, up-regulation of RP11-63A11.1 attenuated cell survival and induced apoptosis.

Conclusions

Ten cancer-related lncRNAs were incorporated into a powerful model for prognosis evaluation. RP11-63A11.1 functioned as a cancer suppressor in pRCC and it might be a potential therapeutic target for treating pRCC.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12935-021-02247-6.

Identification of immunization-related new prognostic biomarkers for papillary renal cell carcinoma by integrated bioinformatics analysis

Background

Despite papillary renal cell carcinoma (pRCC) being the second most common type of kidney cancer, the underlying molecular mechanism remains unclear. Targeted therapies in the past have not been successful because of the lack of a clear understanding of the molecular mechanism. Hence, exploring the underlying mechanisms and seeking novel biomarkers for pursuing a precise prognostic biomarker and appropriate therapies are critical.

Material and methods

In our research, the differentially expressed genes (DEGs) were screened from the TCGA and GEO databases, and a total of 149 upregulated and 285 downregulated genes were sorted. This was followed by construction of functional enrichment and protein–protein interaction (PPI) network, and then the top 15 DEGs were selected for further analysis. The P4HB gene was chosen as our target gene by repetitively validating multiple datasets, and higher levels of P4HB expression predicted lower overall survival (OS) in patients with pRCC.

Results

We found that P4HB not only connects with immune cell infiltration and co-expression with PD-1, PD-L2, and CTLA-4, but also has a strong connection with the newly discovered hot gene, TOX.

Conclusion

We speculate that P4HB is a novel gene involved in the progression of pRCC through immunomodulation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12920-021-01092-w.

Construction of a Signature Composed of 14 Immune Genes to Judge the Prognosis and Immune Infiltration of Colon Cancer

Background: Colon cancer (CC) is an immunogenic tumor and immune-targeting disease. In this study, we analyzed differentially expressed genes (DEGs) from the expression profile data in CC of The Cancer Genome Atlas. Methods and Results: Using univariate and multivariate Cox regression analysis, an immune gene-risk model containing 14 immune genes was established. Four hundred seventeen CC samples were divided into high-risk and low-risk groups, and Kaplan-Meier analysis revealed that high-risk score predicted poor survival. Meanwhile, we found the model was an independent prognostic factor for CC. Weighted gene coexpression network analysis was used to identify key gene modules between high- and low-risk groups. The methods of CIBERSORT and single-sample Gene Set Enrichment Analysis were used to evaluate the correlation between immune cells and our model. Conclusion: Taken together, our study suggested that the immune gene-related risk model may be developed as a potential tool in the prognostic assessment of CC.

A novel immune signature to predict the prognosis of patients with hepatocellular carcinoma

Aberrant immunity has been associated with the initiation and progression of cancers such as hepatocellular carcinoma (HCC). Here, we aim to develop a signature based on immune-related genes (IRGs) to predict the prognosis of HCC patients. The gene expression profiles of 891 HCC samples were derived from 4 publicly accessible datasets. A total of 1534 IRGs from Immunology Database and Analysis Portal website were obtained as candidate genes for prognostic assessment. Using least absolute shrinkage and selection operator (LASSO) regression analysis, 12 IRGs were selected as prognostic biomarkers and were then aggregated to generate an IRG score for each HCC sample. In the training dataset (n = 365), patients with high IRG scores showed a remarkably poorer overall survival than those with low IRG scores (log-rank P < .001). Similar results were documented in 3 independent testing datasets (n = 226, 221, 79, respectively). Multivariate Cox regression and stratified analyses indicated that the IRG score was an independent and robust signature to predict the overall survival in HCC patients. Patients with high IRG scores tended to be in advanced TNM stages, with increased risks of tumor recurrence and metastasis. More importantly, the IRG score was strongly associated with certain immune cell counts, gene expression of immune checkpoints, estimated immune score, and mutation of critical genes in HCC. In conclusion, the proposed IRG score can predict the prognosis and reflect the tumor immune microenvironment of HCC patients, which may facilitate the individualized treatment and provide potential immunotherapeutic targets.

Development of a Predictive Immune-Related Gene Signature Associated With Hepatocellular Carcinoma Patient Prognosis

Background:

Hepatocellular carcinoma (HCC) remains the third leader cancer-associated cause of death globally, but the etiological basis for this complex disease remains poorly clarified. The present study was thus conceptualized to define a prognostic immune-related gene (IRG) signature capable of predicting immunotherapy responsiveness and overall survival (OS) in patients with HCC.

Methods:

Five differentially expressed IRG associated with HCC were established the immune-related risk model through univariate Cox regression and least absolute shrinkage and selection operator (LASSO) regression analyses. Patients were separated at random into training and testing cohorts, after which the association between the identified IRG signature and OS was evaluated using the “survival” R package. In addition, maftools was leveraged to assess mutational data, with tumor mutation burden (TMB) scores being calculated as follows: (total mutations/total bases) × 10⁶. Immune-related risk term abundance was quantified via “ssGSEA” algorithm using the “gsva” R package.

Results:

HCC patients were successfully stratified into low-risk and high-risk groups based upon a signature composed of 5 differentially expressed IRGs, with overall survival being significantly different between these 2 groups in training cohort, testing cohort and overall patient cohort (P = 1.745e-06, P = 1.888e-02, P = 4.281e-07). No association was observed between TMB and this IRG risk score in the overall patient cohort (P = 0.461). Notably, 19 out of 29 immune-related risk terms differed substantially in the overall patient dataset. These risk terms mainly included checkpoints, human leukocyte antigens, natural killer cells, dendritic cells, and major histocompatibility complex class I.

Conclusion:

In summary, an immune-related prognostic gene signature was successfully developed and used to predict survival outcomes and immune system status in patients with HCC. This signature has the potential to help guide immunotherapeutic treatment planning for patients affected by this deadly cancer.

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.

Bioinformatics Profiling of Five Immune-Related lncRNAs for a Prognostic Model of Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC), one of the most common tumors worldwide, has the fifth highest mortality rate, which is increasing every year. At present, many studies have revealed that immunotherapy has an important effect on many malignant tumors. The main purpose of our research was to verify and establish a new immune-related lncRNA model and to explore the potential immune mechanisms. We analysed the pathways and mechanisms of immune-related lncRNAs by bioinformatics analysis, screened key lncRNAs based on Cox regression analysis, and determined the characteristics of the immune-related lncRNAs. On this basis, a predictive model was established. Through a comparison of specificity and sensitivity, we found that the constructed model was superior to the known markers of HCC. Then, the cell types were identified by the relative subgroup (CIBERSORT) algorithm for RNA transcripts. A signature model was eventually constructed, and we proved that it was a survival factor for HCC. Moreover, five kinds of immune cells were significantly positively correlated with the signature. The results indicated that these five kinds of lncRNAs may be related to the immune infiltration of hepatocellular carcinoma. To verify these findings, we selected the top coexpressed lncRNA, AC099850.3, for further study. We found that AC099850.3 could promote the migration and proliferation of hepatocellular carcinoma cells in vitro. RT-PCR experiments found that AC099850.3 could promote the expression of the cell cycle molecules BUB1, CDK1, PLK1, and TTK, and western blotting to prove that the expression of the molecules CD155 and PD-L1 was inhibited in the interference group. In conclusion, we used five kinds of immune-related lncRNAs to construct prognostic signatures to explore the mechanism, which provides a new way to study therapies for HCC.

Multi-Omics Data Analyses Construct a Six Immune-Related Genes Prognostic Model for Cervical Cancer in Tumor Microenvironment

The cross-talk between tumor cells and the tumor microenvironment (TME) is an important factor in determining the tumorigenesis and progression of cervical cancer (CC). However, clarifying the potential mechanisms which trigger the above biological processes remains a challenge. The present study focused on immune-relevant differences at the transcriptome and somatic mutation levels through an integrative multi-omics analysis based on The Cancer Genome Atlas database. The objective of the study was to recognize the specific immune-related prognostic factors predicting the survival and response to immunotherapy of patients with CC. Firstly, eight hub immune-related prognostic genes were ultimately identified through construction of a protein–protein interaction network and Cox regression analysis. Secondly, 32 differentially mutated genes were simultaneously identified based on the different levels of immune infiltration. As a result, an immune gene-related prognostic model (IGRPM), including six factors (chemokine receptor 7 [CCR7], CD3d molecule [CD3D], CD3e molecule [CD3E], and integrin subunit beta 2 [ITGB2], family with sequence similarity 133 member A [FAM133A], and tumor protein p53 [TP53]), was finally constructed to forecast clinical outcomes of CC. Its predictive capability was further assessed and validated using the Gene Expression Omnibus validation set. In conclusion, IGRPM may be a promising prognostic signature to predict the prognoses and responses to immunotherapy of patients with CC. Moreover, the multi-omics study showed that IGRPM could be a novel therapeutic target for CC, which is a promising biomarker for indicating the immune-dominant status of the TME and revealing the potential mechanisms responsible for the tumorigenesis and progression of CC.

Construction and Validation of an Immune Infiltration-Related Gene Signature for the Prediction of Prognosis and Therapeutic Response in Breast Cancer

Breast cancer patients show significant heterogeneity in overall survival. Current assessment models are insufficient to accurately predict patient prognosis, and models for predicting treatment response are lacking. We evaluated the relationship between various immune cells and breast cancer and confirmed the association between immune infiltration and breast cancer progression. Different bioinformatics and statistical approaches were combined to construct a robust immune infiltration-related gene signature for predicting patient prognosis and responses to immunotherapy and chemotherapy. Our research found that a higher immune infiltration-related risk score (IRS) indicates that the patient has a worse prognosis and is not very sensitive to immunotherapy. In addition, a new nomogram was constructed based on the gene signature and clinicopathological features to improve the risk stratification and quantify the risk assessment of individual patients. Our study might contribute to the optimization of the risk stratification for survival and the personalized management of breast cancer.

Construction of a prognostic value model in papillary renal cell carcinoma by immune-related genes

Papillary renal cell carcinoma (PRCC) is the second most common type of renal carcinoma following clear cell renal cell carcinoma, and the role of immune-related genes (IRGs) in tumorigenesis and metastasis is evident; its prognostic value in PRCC remains unclear. In this study, we downloaded the gene expression profiles and clinical data of patients with PRCC from The Cancer Genome Atlas (TCGA) database and obtained IRGs from the ImmPort database. A total of 371 differentially expressed IRGs (DEIRGs) were discovered between PRCC and normal kidney tissues. Prognostic DEIRGs (PDEIRGs) were identified by univariate Cox regression analysis. Then, we screened the four most representative PDEIRGs (IL13RA2, CCL19, BIRC5, and INHBE) and used them to construct a risk model to predict the prognosis of patients with PRCC. This model precisely stratified survival outcome and accurately identified mutation burden in PRCC. Thus, our results suggest that these four PDEIRGs are available prognostic predictors for PRCC. They could be used to assess the prognosis and to guide individualized treatments for patients with PRCC.

Identification of an Immune-Related Risk Signature Correlates With Immunophenotype and Predicts Anti-PD-L1 Efficacy of Urothelial Cancer

Immune checkpoint inhibitor (ICI) treatment has been used to treat advanced urothelial cancer. Molecular markers might improve risk stratification and prediction of ICI benefit for urothelial cancer patients. We analyzed 406 cases of bladder urothelial cancer from The Cancer Genome Atlas (TCGA) data set and identified 161 messenger RNAs (mRNAs) as differentially expressed immunity genes (DEIGs). Using the LASSO Cox regression model, an eight-mRNA-based risk signature was built. We validated the prognostic and predictive accuracy of this immune-related risk signature in 348 metastatic urothelial cancer (mUC) samples treated with anti-PD-L1 (atezolizumab) from IMvigor210. We built an immune-related risk signature based on the eight mRNAs: ANXA1, IL22, IL9R, KLRK1, LRP1, NRG3, SEMA6D, and STAP2. The eight-mRNA-based risk signature successfully categorizes patients into high-risk and low-risk groups. Overall survival was significantly different between these groups, regardless if the initial TCGA training set, the internal TCGA testing set, all TCGA set, or the ICI treatment set. The hazard ratio (HR) of the high-risk group to the low-risk group was 3.65 (p < 0.0001), 2.56 (p < 0.0001), 3.36 (p < 0.0001), and 2.42 (p = 0.0009). The risk signature was an independent prognostic factor for prediction survival. Moreover, the risk signature was related to immunity characteristics. In different tumor mutational burden (TMB) subgroups, it successfully categorizes patients into high-risk and low-risk groups, with significant differences of clinical outcome. Our eight-mRNA-based risk signature is a stable biomarker for urothelial cancer and might be able to predict which patients benefit from ICI treatment. It might play a role in precision individualized immunotherapy.

Expression patterns and prognostic value of m6A RNA methylation regulators in adrenocortical carcinoma

Adrenocortical carcinoma (ACC) is considered a rare cancer with poor prognosis. We used public datasets from The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) databases to assess the relationships between N6-methyladenosine (m6A)-related genes and ACC.We used the Wilcoxon signed-rank test to compare m6A-related gene expression in ACC tissues with that in normal tissues. Then, ACC patients were grouped based on a cluster analysis of m6A-related gene expression. m6A-related genes that were significantly associated with survival were incorporated into a risk signature, and 2 groups were divided according to median risk score. Fisher exact tests were utilized to analyze differences in clinical variables between groups. We compared the overall survival (OS) rates of the groups by means of Kaplan-Meier curves and Cox regression analyses.We found that RBM15, ZC3H3, YTDHF1, YTDHF2, and ALBH5 were overexpressed in ACC and that KIAA1429, YTHDC1, HNRNPC, WTAP, METTL3, and FTO were down regulated in ACC. In addition, membership in cluster 2 or the high-risk group was associated with advanced clinical factors and poor prognosis. The univariable and multivariable Cox regression analyses showed that risk score can be considered an independent prognostic factor for ACC.We found that the expression of m6A-related genes could be used as an independent prognostic factor in ACC. However, the current study has some limitations, and further studies of m6A-related genes in ACC are needed.

Development and verification of an immune-related gene pairs prognostic signature in ovarian cancer

Ovarian cancer (OV) is the most common gynaecological cancer worldwide. Immunotherapy has recently been proven to be an effective treatment strategy. The work here attempts to produce a prognostic immune-related gene pair (IRGP) signature to estimate OV patient survival. The Gene Expression Omnibus (GEO) and Cancer Genome Atlas (TCGA) databases provided the genetic expression profiles and clinical data of OV patients. Based on the InnateDB database and the least absolute shrinkage and selection operator (LASSO) regression model, we first identified a 17-IRGP signature associated with survival. The average area under the curve (AUC) values of the training, validation, and all TCGA sets were 0.869, 0.712, and 0.778, respectively. The 17-IRGP signature noticeably split patients into high- and low-risk groups with different prognostic outcomes. As suggested by a functional study, some biological pathways, including the Toll-like receptor and chemokine signalling pathways, were significantly negatively correlated with risk scores; however, pathways such as the p53 and apoptosis signalling pathways had a positive correlation. Moreover, tumour stage III, IV, grade G1/G2, and G3/G4 samples had significant differences in risk scores. In conclusion, an effective 17-IRGP signature was produced to predict prognostic outcomes in OV, providing new insights into immunological biomarkers.

The role of an immune signature for prognosis and immunotherapy response in endometrial cancer

Immunotherapy is a practical and promising treatment for advanced and recurrent endometrial cancer (EC). In this study, we identified an immune-related gene (IRG) signature to predict the overall survival (OS) and response to immune checkpoints inhibitors (ICIs) in patients with EC. The RNA expression profiles of EC were obtained from The Cancer Genome Atlas database and then were filtered for IRGs based on the Immport database. Using the conjoint Cox regression model, an immune signature consisting of seven risk IRGs (CBLC, PLA2G2A, TNF, NR3C1, APOD, TNFRSF18, and LTB) was developed. The immune signature was independent of other clinical factors and was superior to the traditional staging method for OS prediction in EC. Immunohistochemistry staining from the Human Protein Atlas database and quantitative real-time PCR analysis of EC samples were also performed to validate the expression levels of risk IRGs. By further analyzing the tumor microenvironment in EC, patients in the low-risk subgroup showed a higher immune cell infiltration status, which was associated with a better prognosis. Moreover, the tumor mutational burden and immunophenoscore analysis demonstrated that the low-risk subgroup was more sensitive to ICI-based immunotherapy. These findings might shed light on the development of targeted treatment and novel biomarkers for patients with EC.

Integrated Analysis of a Competing Endogenous RNA Network Reveals a Prognostic Signature in Kidney Renal Papillary Cell Carcinoma

The kidney renal papillary cell carcinoma (KIRP) is a relatively rare type of renal cell carcinoma (RCC). Currently, most kidney cancer studies primarily focus on RCC, and there has been no investigation to find a robust signature to predict the survival outcome of KIRP patients. In this study, we constructed a competing endogenous RNA (ceRNA) network, including 1,251 lncRNA-miRNA-mRNA interactions. Eight differentially expressed genes (IGF2BP3, PLK1, LINC00200, NCAPG, CENPF, miR-217, GAS6-As1, and LRRC4) based on the TCGA database were selected. The prognostic signature was established by combining the univariate Cox regression method and a stepwise regression method, with its predictive value validated by time-dependent receiver operating characteristic (ROC) curves. In conclusion, we identified eight prognostic signatures with using ceRNA networks. Our study provided a global view and a systematic dissection on KIRP prognosis biomarkers, and the eight identified genes might be used as new and important prognostic factors involved in KIRP pathogenesis.

Identification of 6 Hub Proteins and Protein Risk Signature of Colorectal Cancer

Background

Colorectal cancer (CRC) is the second most common cause of cancer death in the United States and the third most common cancer globally. The incidence of CRC tends to be younger, and we urgently need a reliable prognostic assessment strategy.

Methods

Protein expression profile and clinical information of 390 CRC patients/samples were downloaded from the TCPA and TCGA database, respectively. The Kaplan-Meier, Cox regression, and Pearson correlation analysis were applied in this study.

Results

Based on the TCPA and TCGA database, we screened 6 hub proteins and first constructed protein risk signature, all of which were significantly associated with CRC patients' overall survival (OS). The risk score was an independent prognostic factor and significantly related with the size of the tumor in situ (T). 6 hub proteins were differentially expressed in cancer and normal tissues and in different CRC stages, which were validated at the ONCOMINE database. Next, 40 coexpressed proteins of 6 hub proteins were extracted from the TCPA database. In the protein-protein interaction (PPI) network, HER1, HER2, and CTNNB1 were at the center. Function enrichment analysis illustrated that 46 proteins were mainly involved in the EGFR (HER1) tyrosine kinase inhibitor resistance pathway.

Conclusion

Studies indicated that 6 hub proteins might be considered as new targets for CRC therapies, and the protein risk signature can be used to predict the OS of CRC patients.

Survival prediction and response to immune checkpoint inhibitors: A prognostic immune signature for hepatocellular carcinoma

•

As is known to us, this is the first immune-gene-related signature of HCC which was validated in various ways.

•

This signature was linked to the mutation status, and we found the difference of mutation between two risk groups.

•

We estimated the association among the signature, TMB and the survival and found the result contrary to previous studies.

Hepatocellular carcinoma (HCC) is one of the most common cancers all over the world. Several studies have explored if immune-related genes and tumor immune microenvironment could play roles in HCC prognoses. This study is aimed at developing a prognostic signature of HCC based on immune-related genes or tumor immune microenvironment to predict survival and response to immune checkpoint inhibitors (ICIs). We constructed a prognostic signature using bioinformatics method and validated its predictive capability. The mechanisms of the signature prediction were explored with The Cancer Immunome Atlas (TCIA) and mutation analysis. We also explored the association between the signature and immunophenoscore (IPS), which is the marker of ICIs response. A 6 immune-related-gene (6-IRG) signature was developed. It was revealed in a multivariate analysis that the 6-IRG signature was an independent prognostic factor of overall survival and progression-free interval among HCC patients. In the high-risk group of 6-IRG signature score, macrophage M0 cells and regulatory T cells, which are observed associated with poor overall survival in our study, were higher. The low-risk group had a higher IPS, which meant a better response to ICIs. Taken together, we constructed a reliable 6-IRG signature for prediction of survival and response to ICIs. The signature needs further testing for clinical application.

Development and Validation of an Individualized Immune-Related Gene Pairs Prognostic Signature in Papillary Renal Cell Carcinoma

Papillary renal carcinoma (PRCC) is one of the important subtypes of kidney cancer, with a high degree of heterogeneity. At present, there is still a lack of robust and accurate biomarkers for the diagnosis, prognosis and treatment selection of PRCC. Considering the important role of tumor immunity in PRCC, we aim to construct a signature based on immune-related gene pairs (IRGPs) to estimate the prognostic of patients with PRCC. We obtained gene expression profiling and clinical information of patients with PRCC from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO), which were divided into discovery (n = 287) and validation (n = 28) cohorts, respectively. By univariate analysis, multivariate Cox analysis, and least absolute shrinkage and selection operator (Lasso) analysis, we selected 14 IRGPs with a panel of 22 unique genes to construct the prognostic signature. According to the signature, we stratified patients into high-risk group and low-risk group. In both discovery and validation cohorts, the results of Kaplan-Meier analysis showed that there were significant differences in OS between the two groups (p < 0.001). Combined with multiple clinical and pathological factors, the results of multivariate analyses confirmed that this signature was an independent predictor of OS (HR, 3.548; 95%CI, 2.096-6.006; p < 0.001). The results of immune infiltration analysis demonstrated that the abundance of multiple tumor-infiltration lymphocytes such as CD8 + T cells, Tregs, and T follicular cell helper were significantly higher in the high-risk group. Functional analysis showed that multiple immune-related signaling pathways were enriched in the high-risk group. In conclusion, we successfully established an individualized prognostic IRGPs signature, which can accurately assess and predict the OS of patients with PRCC.

Comprehensive profiling of immune-related genes in soft tissue sarcoma patients

Background

Immune-related genes (IRGs) have been confirmed to have an important role in tumorigenesis and tumor microenvironment formation. Nevertheless, a systematic analysis of IRGs and their clinical significance in soft tissue sarcoma (STS) patients is lacking.

Methods

Gene expression files from The Cancer Genome Atlas (TCGA) database and Genotype-Tissue Expression (GTEx) were used to select differentially expressed genes (DEGs). Differentially expressed immune-related genes (DEIRGs) were determined by matching the DEG and ImmPort gene sets, which were evaluated by functional enrichment analysis. Unsupervised clustering of the identified DEIRGs was conducted, and associations with prognosis, the tumor microenvironment (TME), immune checkpoints, and immune cells were analyzed simultaneously. Two prognostic signatures, one for overall survival (OS) and one for progression free survival (PFS), were established and validated in an independent set. Finally, two transcription factor (TF)-IRG regulatory networks were constructed, and a crucial regulatory axis was validated.

Results

In total, 364 DEIRGs and four clusters were identified. OS, TME scores, five immune checkpoints, and 12 types of immune cells were found to be significantly different among the four clusters. The two prognostic signatures incorporating 20 DEIRGs showed favorable discrimination and were successfully validated. Two nomograms combining signature and clinical variables were generated. The C-indexes were 0.879 (95%CI 0.832 ~ 0.926) and 0.825 (95%CI 0.776 ~ 0.874) for the OS and PFS signatures, respectively. Finally, TF-IRG regulatory networks were established, and the MYH11-ADM regulatory axis was verified in three independent datasets.

Conclusion

This comprehensive analysis of the IRG landscape in soft tissue sarcoma revealed novel IRGs related to carcinogenesis and the immune microenvironment. These findings have implications for prognosis and therapeutic responses, which reveal novel potential prognostic biomarkers, promote precision medicine, and provide potential novel targets for immunotherapy.

Development of an Immune-Related Risk Signature for Predicting Prognosis in Lung Squamous Cell Carcinoma

Lung squamous cell carcinoma (LSCC) is the most common subtype of non-small cell lung cancer. Immunotherapy has become an effective treatment in recent years, while patients showed different responses to the current treatment. It is vital to identify the potential immunogenomic signatures to predict patient' prognosis. The expression profiles of LSCC patients with the clinical information were downloaded from TCGA database. Differentially expressed immune-related genes (IRGs) were extracted using edgeR algorithm, and functional enrichment analysis showed that these IRGs were primarily enriched in inflammatory- and immune-related processes. "Cytokine-cytokine receptor interaction" and "PI3K-AKT signaling pathway" were the most enriched KEGG pathways. 27 differentially expressed IRGs were significantly correlated with the overall survival (OS) of patients using univariate Cox regression analysis. A prognostic risk signature that comprises seven IRGs (GCCR, FGF8, CLEC4M, PTH, SLC10A2, NPPC, and FGF4) was developed with effective predictive performance by multivariable Cox stepwise regression analysis. Most importantly, the signature could be an independent prognostic predictor after adjusting for clinicopathological parameters, and also validated in two independent LSCC cohorts (GSE4573 and GSE17710). Potential molecular mechanisms and tumor immune landscape of these IRGs were investigated through computational biology. Analysis of tumor infiltrating lymphocytes and immune checkpoint molecules revealed distinct immune landscape in high- and low-risk group. The study was the first time to construct IRG-based immune signature in the recognition of disease progression and prognosis of LSCC patients.

Seven immune-related genes prognostic power and correlation with tumor-infiltrating immune cells in hepatocellular carcinoma

Background

Given poor prognosis and the lack of efficient therapy for advanced hepatocellular carcinoma, immunotherapy has emerged as an increasingly important role. However, there were few reports on the correlation between immune‐related genes and HCC. The purpose of this study is to construct a novel immune‐related gene‐based prognostic signature for HCC and to explore the potential mechanisms.

Methods

We organized expression data of 374 HCC samples and 50 nontumor samples from TCGA database. A robust signature was constructed by Cox regression analysis based on the immune‐related genes, which were filtered by differential genes analysis and Cox regression analysis. Then, the correlation analysis between the signature and clinical characteristics was conducted. And the signature was validated in ICGC database. Furthermore, the relationships between immune cell infiltration and the signature were explored by bioinformatics analysis.

Results

Seven genes‐based model (Risk score = BIRC5 * 0.0238 + FOS * 0.0055 + DKK1 * 0.0085 + FGF13 * 0.3432 + IL11 * 0.0135 + IL17D * 0.0878 + SPP1 * 0.0003) was constructed eventually and it was proved to be an independent prognostic factor for HCC patients. The signature‐calculated risk scores were shown to be positively correlated with the infiltration of these five immune cells, including macrophages, neutrophils, CD8+T, dendritic, and B cells. And the results suggested that high amplication of BIRC5, FGF13, IL11, IL17D, and SPP1 were more likely correlated with immune cell infiltration. Finally, PPI network, TFs‐based regulatory network and gene enrichment plots were performed to show potential molecular mechanisms.

Conclusion

We construct a robust immune‐related gene‐based prognostic signature with seven genes and explore potential mechanisms about it, which may contribute to the immunotherapy research for HCC.

Identification of core miRNA prognostic markers in patients with laryngeal cancer using bioinformatics analysis

PURPOSE

Lots of studies indicated that many microRNAs (miRNAs) are associated with the prognosis of patients with laryngeal cancer (LC). The objective of our study is to identify potential core miRNAs associated with the pathogenesis and prognosis of LC.

METHODS

Using the Cancer Genome Atlast database, we identified 70 differentially expressed miRNAs between LC tumor specimens and non-tumor specimens. Then Cox regression analyses and the least absolute shrinkage and selection operator regression signature were performed to detect miRNA prognostic markers. A nomogram integrating miRNA prognostic markers was constructed to predict overall survival (OS) for LC patients. The potential target genes of the key miRNA were predicted by miRTarBase and miRDB databases. Subsequently, their potential functions were revealed by gene ontology annotation and kyoto encyclopedia of genes and genomes pathway enrichment analysis. Related biological pathways of the key target gene involved in LC were detected through gene set enrichment analysis (GSEA).

RESULTS

A prognostic miRNA signature was constructed. The up-regulated miR-105-1 was related to a worse OS (p = 0.043), which suggested that miR-105-1 may likely be the key miRNA prognostic marker. Survival analyses and paired expression analyses of target genes indicated that ENDOU may be the key target gene. Finally, we conducted GSEA to elucidate the pathways enriched between low- and high-ENDOU expression datasets.

CONCLUSION

Our findings might bring some new light on the pathogenesis of LC. Then, it might facilitate doctors to predict the prognosis and improve treatment outcomes for LC patients. However, the behaviors of LC are relatively heterogeneous, and the TCGA database cannot provide detailed information about the subsites and treatment modalities of LC. Further molecular biological experiments and clinical investigations would be required to confirm this conclusion.

A novel prognostic signature of immune-related genes for patients with colorectal cancer

Colorectal cancer (CRC) is one of the most commonly diagnosed cancers with an estimated 1.8 million new cases worldwide and associated with high mortality rates of 881 000 CRC-related deaths in 2018. Screening programs and new therapies have only marginally improved the survival of CRC patients. Immune-related genes (IRGs) have attracted attention in recent years as therapeutic targets. The aim of this study was to identify an immune-related prognostic signature for CRC. To this end, we combined gene expression and clinical data from the CRC data sets of The Cancer Genome Atlas (TCGA) into an integrated immune landscape profile. We identified a total of 476 IRGs that were differentially expressed in CRC vs normal tissues, of which 18 were survival related according to univariate Cox analysis. Stepwise multivariate Cox proportional hazards analysis established an immune-related prognostic signature consisting of SLC10A2, FGF2, CCL28, NDRG1, ESM1, UCN, UTS2 and TRDC. The predictive ability of this signature for 3- and 5-year overall survival was determined using receiver operating characteristics (ROC), and the respective areas under the curve (AUC) were 79.2% and 76.6%. The signature showed moderate predictive accuracy in the validation and GSE38832 data sets as well. Furthermore, the 8-IRG signature correlated significantly with tumour stage, invasion, lymph node metastasis and distant metastasis by univariate Cox analysis, and was established an independent prognostic factor by multivariate Cox regression analysis for CRC. Gene set enrichment analysis (GSEA) revealed a relationship between the IRG prognostic signature and various biological pathways. Focal adhesions and ECM-receptor interactions were positively correlated with the risk scores, while cytosolic DNA sensing and metabolism-related pathways were negatively correlated. Finally, the bioinformatics results were validated by real-time RT-qPCR. In conclusion, we identified and validated a novel, immune-related prognostic signature for patients with CRC, and this signature reflects the dysregulated tumour immune microenvironment and has a potential for better CRC patient management.

Identification of Immune-Related Prognostic Genes and LncRNAs Biomarkers Associated With Osteosarcoma Microenvironment

Osteosarcoma (OS) is the most common malignancy of the bone that occurs majorly in young people and adolescents. Although the survival of OS patients markedly improved by complete surgical resection and chemotherapy, the outcome is still poor in patients with recurrent and/or metastasized OS. Thus, identifying prognostic biomarkers that reflect the biological heterogeneity of OS could lead to better interventions for OS patients. Increasing studies have indicated the association between immune-related genes (IRGs) and cancer prognosis. In the present study, based on the data concerning OS obtained from TARGET (Therapeutically Applicable Research to Generate Effective Treatments) database, we constructed a classifier containing 12 immune-related (IR) long non-coding RNAs (lncRNAs) and 3 IRGs for predicting the prognosis of OS by using the least absolute shrinkage and selection operation Cox regression. Besides, based on the risk score calculated by the classifier, the samples were divided into high- and low-risk groups. We further investigated the tumor microenvironment of the OS samples by ESTIMATE and CIBERSORT algorithms between the two groups. Finally, we identified three small molecular drugs with potential therapeutic value for OS patients with high-risk score. Our results suggest that the IRGs and IR-lncRNAs–based classifier could be used as a reliable prognostic predictor for OS survival.

Development of an immune-related prognostic model for pediatric acute lymphoblastic leukemia patients

Background

Acute lymphoblastic leukemia (ALL) is the most common hematological malignancy in pediatrics, and immune‐related genes (IRGs) play crucial role in its development. Our study aimed to identify prognostic immune biomarkers of pediatric ALL and construct a risk assessment model.

Methods

Pediatric ALL patients’ gene expression data were downloaded from Therapeutically Applicable Research to Generate Effective Treatments (TARGET) database. We screened differentially expressed IRGs (DEIRGs) between the relapse and non‐relapse groups. Cox regression analysis was used to identify optimal prognostic genes, then, a risk model was constructed, and its accuracy was verified in different cohorts.

Results

We screened 130 DEIRGs from 251 pediatric ALL samples. The top three pathways that DEIRGs may influence tumor progression are NABA matrisome‐associated, chemotaxis, and antimicrobial humoral response. A set of 84 prognostic DEIRGs was identified by using univariate Cox analysis. Then, Lasso regression and multivariate Cox regression analysis screened four optimal genes (PRDX2, S100A10, RORB, and SDC1), which were used to construct the prognostic risk model. The risk score was calculated and the survival analysis results showed that high‐risk score was associated with poor overall survival (OS) (p = 3.195 × 10⁻⁷). The time‐dependent survival receiver operating characteristic curves showed good prediction accuracy (Area Under Curves for 3‐year, 5‐year OS were 0.892 and 0.89, respectively). And the predictive performance of our risk model was successfully verified in testing cohort and entire cohort.

Conclusions

Our prognostic risk model can effectively divide pediatric ALL patients into high‐risk and low‐risk groups, which may help predict clinical prognosis and optimize individualized treatment.

Development of an Immune Infiltration-Related Prognostic Scoring System Based on the Genomic Landscape Analysis of Glioblastoma Multiforme

Introduction: Glioblastoma multiforme (GBM) is the most common deadly brain malignancy and lacks effective therapies. Immunotherapy acts as a promising novel strategy, but not for all GBM patients. Therefore, classifying these patients into different prognostic groups is urgent for better personalized management.

Materials and Methods: The Cell type Identification by Estimating Relative Subsets of RNA Transcripts (CIBERSORT) algorithm was used to estimate the fraction of 22 types of immune-infiltrating cells, and least absolute shrinkage and selection operator (LASSO) Cox regression analysis was performed to construct an immune infiltration-related prognostic scoring system (IIRPSS). Additionally, a quantitative predicting survival nomogram was also established based on the immune risk score (IRS) derived from the IIRPSS. Moreover, we also preliminarily explored the differences in the immune microenvironment between different prognostic groups.

Results: There was a total of 310 appropriate GBM samples (239 from TCGA and 71 from CGGA) included in further analyses after CIBERSORT filtering and data processing. The IIRPSS consisting of 17 types of immune cell fractions was constructed in TCGA cohort, the patients were successfully classified into different prognostic groups based on their immune risk score (p = 1e-10). What's more, the prognostic performance of the IIRPSS was validated in CGGA cohort (p = 0.005). The nomogram also showed a superior predicting value. (The predicting AUC for 1-, 2-, and 3-year were 0.754, 0.813, and 0.871, respectively). The immune microenvironment analyses reflected a significant immune response and a higher immune checkpoint expression in high-risk immune group.

Conclusion: Our study constructed an IIRPSS, which maybe valuable to help clinicians select candidates most likely to benefit from immunological checkpoint inhibitors (ICIs) and laid the foundation for further improving personalized immunotherapy in patients with GBM.

Identification of an immune-related risk signature for predicting prognosis in clear cell renal cell carcinoma

Immune status affects the initiation and progression of clear cell renal cell carcinoma (ccRCC), the most common subtype of renal cell carcinoma. In this study, we identified an immune-related, five-gene signature that improves survival prediction in ccRCC. Patients were classified as high- and low-risk based on the signature risk score. Survival analysis showed differential prognosis, while principal component analysis revealed distinctly different immune phenotypes between the two risk groups. High-risk patients tended to have advanced stage, higher grade disease, and poorer prognoses. Functional enrichment analysis showed that the signature genes were mainly involved in the cytokine-cytokine receptor interaction pathway. Moreover, we found that tumors from high-risk patients had higher relative abundance of T follicular helper cells, regulatory T cells, and M0 macrophages, and higher expression of PD-1, CTLA-4, LAG3, and CD47 than low-risk patients. This suggests our gene signature may not only serve as an indicator of tumor immune status, but may be a promising tool to select high-risk patients who may benefit from immune checkpoint inhibitor therapy. Multivariate Cox regression analysis showed that the signature remained an independent prognostic factor after adjusting for clinicopathological variables, while prognostic accuracy was further improved after integrating clinical parameters into the analysis.

Prognostic value of immune-related genes in clear cell renal cell carcinoma

Clear cell renal cell carcinoma (ccRCC) is the most common pathological subtype of renal cell carcinoma, and immune-related genes (IRGs) are key contributors to its development. In this study, the gene expression profiles and clinical data of ccRCC patients were downloaded from The Cancer Genome Atlas database and the cBioPortal database, respectively. IRGs were obtained from the ImmPort database. We analyzed the expression of IRGs in ccRCC, and discovered 681 that were differentially expressed between ccRCC and normal kidney tissues. Univariate Cox regression analysis was used to identify prognostic differentially expressed IRGs (PDEIRGs). Using Lasso regression and multivariate Cox regression analyses, we detected seven optimal PDEIRGs (PLAU, ISG15, IRF9, ARG2, RNASE2, SEMA3G and UCN) and used them to construct a risk model to predict the prognosis of ccRCC patients. This model accurately stratified patients with different survival outcomes and precisely identified patients with different mutation burdens. Our findings suggest the seven PDEIRGs identified in this study are valuable prognostic predictors in ccRCC patients. These genes could be used to investigate the developmental mechanisms of ccRCC and to design individualized treatments for ccRCC patients.

The clinical significance of collagen family gene expression in esophageal squamous cell carcinoma

BACKGROUND

Esophageal squamous cell carcinoma (ESCC) is a subtype of esophageal cancer with high incidence and mortality. Due to the poor 5-year survival rates of patients with ESCC, exploring novel diagnostic markers for early ESCC is emergent. Collagen, the abundant constituent of extracellular matrix, plays a critical role in tumor growth and epithelial-mesenchymal transition. However, the clinical significance of collagen genes in ESCC has been rarely studied. In this work, we systematically analyzed the gene expression of whole collagen family in ESCC, aiming to search for ideal biomarkers.

METHODS

Clinical data and gene expression profiles of ESCC patients were collected from The Cancer Genome Atlas and the gene expression omnibus databases. Bioinformatics methods, including differential expression analysis, survival analysis, gene sets enrichment analysis (GSEA) and co-expression network analysis, were performed to investigate the correlation between the expression patterns of 44 collagen family genes and the development of ESCC.

RESULTS

A total of 22 genes of collagen family were identified as differentially expressed genes in both the two datasets. Among them, COL1A1, COL10A1 and COL11A1 were particularly up-regulated in ESCC tissues compared to normal controls, while COL4A4, COL6A5 and COL14A1 were notably down-regulated. Besides, patients with low COL6A5 expression or high COL18A1 expression showed poor survival. In addition, a 7-gene prediction model was established based on collagen gene expression to predict patient survival, which had better predictive accuracy than the tumor-node-metastasis staging based model. Finally, GSEA results suggested that collagen genes might be tightly associated with PI3K/Akt/mTOR pathway, p53 pathway, apoptosis, cell cycle, etc.

CONCLUSION

Several collagen genes could be potential diagnostic and prognostic biomarkers for ESCC. Moreover, a novel 7-gene prediction model is probably useful for predicting survival outcomes of ESCC patients. These findings may facilitate early detection of ESCC and help improves prognosis of the patients.

Identification of a prognostic immune signature for cervical cancer to predict survival and response to immune checkpoint inhibitors

Cervical cancer (CC) is a leading cause of cancer-related death in women. Limited studies have investigated whether immune-related genes (IRGs) or tumor immune microenvironment (TIME) could be indicators for CC prognoses. The aim of this study was to develop an improved prognostic signature for CC based on IRGs or TIME to predict survival and response to immune checkpoint inhibitors (ICIs). A prognostic signature was constructed using bioinformatics method and its predictive capability was validated. The mechanisms underlying the signature's predictive capability were explored with CIBERSORT algorithm and mutation analysis. Immunophenoscore (IPS) is validated for ICIs response, and was therefore explored in relation to the signature. A prognostic signature based on 11 IRGs was developed. A multivariate analysis revealed that the 11-IRG signature was an independent prognostic factor for overall survival (OS) and progression-free interval in CC patients. In the 11-IRG signature high-risk group, CD8 T cells and resting mast cells, which are found to associate with better OS in our study, were lower; activated mast cells, associated with poorer OS, were higher, compared with the low-risk group. An IPS analysis suggested that the 11-IRG signature low-risk group, which possessed a higher IPS, represented a more immunogenic phenotype that was more inclined to respond to ICIs. In short, an 11-IRG prognostic signature for predicting CC patients' survival and response to ICIs was firmly established. The predictive capability of this model in CC requires further testing with the goal of better prognostic stratification and treatment management.

Identification of novel prognostic alternative splicing signature in papillary renal cell carcinoma

Papillary renal cell carcinoma (pRCC) is a heterogeneous disease containing multifocal or solitary tumors with an aggressive phenotype. Increasing evidence has indicated the involvement of aberrant splicing variants in renal cell cancer, while systematic profiling of aberrant alternative splicing (AS) in pRCC was lacking and largely unknown. In the current study, comprehensive profiling of AS events were performed based on the integration of pRCC cohort from the Cancer Genome Atlas database and SpliceSeq software. With rigorous screening and univariate Cox analysis, a total of 2077 prognoses AS events from 1642 parent genes were identified. Then, stepwise least absolute shrinkage and selection operator method-penalized Cox regression analyses with 10-fold cross-validation followed by multivariate Cox regression were used to construct the prognostic AS signatures within each AS type. And a final 21 AS event-based signature was proposed which showed potent prognostic capability in stratifying patients into low- and high-risk subgroups (P < .0001). Furthermore, time-dependent receiver operating characteristics curves confirmed that the final AS signature was effective and robust in predicting overall survival for pRCC patients with the area under the curve above 0.9 from 1 to 5 years. In addition, splicing correlation network was built to uncover the potential regulatory pattern among prognostic splicing factors and candidate AS events. Besides, gene set enrichment analysis revealed the involvement of these candidates AS events in tumor-related pathways including extracellular matrix organization, oxidative phosphorylation, and P53 signaling pathways. Taken together, our results could contribute to elucidating the underlying mechanism of AS in the oncogenesis process and broaden the novel field of prognostic and clinical application of molecule-targeted approaches in pRCC.

Construction of immune-related risk signature for renal papillary cell carcinoma

The kidney renal papillary cell carcinoma (KIRP) is a relatively rare type of kidney cancer. There has been no investigation to find a robust signature to predict the survival outcome of KIRP patients in the aspect of tumor immunology. In this study, 285 KIRP samples from The Cancer Genome Atlas (TCGA) were randomly divided into training and testing set. A total of 1534 immune‐related genes from The Immunology Database and Analysis Portal (ImmPort) were used as candidates to construct the signature. Using univariate Cox analysis, we evaluated the relationship between overall survival and immune‐related genes expression and found 272 immune‐related genes with predicting prognostic ability. In order to construct an efficient predictive model, the Cox proportional hazards model with an elastic‐net penalty was used and identified 23 groups after 1000 iterations. As a result, 15‐genes model showing more stable than other gene groups was chosen to construct our immune‐related risk signature. In line with our expectations, the signature can independently predict the survival outcome of KIRP patients. Patients with high‐immune risk were found correlated with advanced stage. We also found that the high‐immune risk patients with higher PBRM1 and SETD2 mutations, increasing chromosomal instability, together with the gene set enrichment analysis (GSEA) results showing intensive connection of our signature with immune pathways. In conclusion, our study constructs a robust 15‐gene signature for predicting KIRP patients’ survival outcome on the basis of tumor immune environment and may provide possible relationship between prognosis and immune‐related biological function.