Identification of immune-enhanced molecular subtype associated with BRCA1 mutations, immune checkpoints and clinical outcome in ovarian carcinoma

A prognostic model for hepatocellular carcinoma patients based on polyunsaturated fatty acid-related genes

OBJECTIVE

Polyunsaturated fatty acids (PUFAs) have attracted increasing attention for their role in liver cancer development. The objective of this study is to develop a prognosis prediction model for patients with liver cancer based on PUFA-related metabolic gene characteristics.

METHOD

Transcriptome data and clinical data were obtained from public databases, while gene sets related to PUFAs were acquired from the gene set enrichment analysis (GSEA) database. Univariate Cox analysis was conducted on the training set, followed by LASSO logistic regression and multivariate Cox analysis on genes with p < .05. Subsequently, the stepwise Akaike information criterion method was employed to construct the model. The high- and low-risk groups were divided based on the median score, and the model's survival prediction ability, diagnostic efficiency, and risk score distribution of clinical features were validated. The above procedures were also validated in the validation set. Immune infiltration levels were evaluated using four algorithms, and the immunotherapeutic potential of different groups was explored. Significant enrichment pathways among different groups were selected based on the GSEA algorithm, and mutation analyses were conducted. Nomogram prognostic models were constructed by incorporating clinical factors and risk scores using univariate and multivariate Cox regression analysis, validated through calibration curves and clinical decision curves. Additionally, sensitivity analysis of drugs was performed to screen potential targeted drugs.

RESULTS

We constructed a prognostic model comprising eight genes (PLA2G12A, CYP2C8, ABCCI, CD74, CCR7, P2RY4, P2RY6, and YY1). Validation across multiple datasets indicated the model's favorable prognostic prediction ability and diagnostic efficiency, with poorer grading and staging observed in the high-risk group. Variations in mutation status and pathway enrichment were noted among different groups. Incorporating Stage, Grade, T.Stage, and RiskScore into the nomogram prognostic model demonstrated good accuracy and clinical decision benefits. Multiple immune analyses suggested greater benefits from immunotherapy in the low-risk group. We predicted multiple targeted drugs, providing a basis for drug development.

CONCLUSION

Our study's multifactorial prognostic model across multiple datasets demonstrates good applicability, offering a reliable tool for personalized therapy. Immunological and mutation-related analyses provide theoretical foundations for further research. Drug predictions offer important insights for future drug development and treatment strategies. Overall, this study provides comprehensive insights into tumor prognosis assessment and personalized treatment planning.

Innovative molecular subtypes of multiple signaling pathways in colon cancer and validation of FMOD as a prognostic-related marker

PURPOSE

Colon cancer is highly heterogeneous in terms of the immune and stromal microenvironment, genomic integrity, and oncogenic properties; therefore, molecular subtypes of the four characteristic dimensions are expected to provide novel clues for immunotherapy of colon cancer.

METHODS

According to the enrichment of four dimensions, we performed consensus cluster analysis and identified three robust molecular subtypes for colon cancer, namely immune enriched, immune deficiency, and stroma enriched. We characterized and validated the immune infiltration, gene mutations, copy number variants, methylation, protein expression, and clinical features in different datasets. Finally, we developed an 8-gene risk prognostic model and proposed the innovative RiskScore. In addition, a nomogram model was constructed combining clinical characteristics and RiskScore to validate its excellent clinical predictive power.

RESULTS

Combining clinical patient tissue samples and histochemical microarray data, we found that high FMOD expression in tumor epithelial cells was associated with poorer patient prognosis, but FMOD expression in the mesenchyme was not associated with prognosis. In pan-cancer, RiskScore, a prognostic model constructed based on characteristic pathway scores, was a poor prognostic factor for malignancy and was negatively associated with immunotherapy response.

CONCLUSION

The identification of molecular subtypes could provide innovative ideas for immunotherapy of colon cancer.

Prognostic signature construction and immunotherapy response analysis for Uterine Corpus Endometrial Carcinoma based on cuproptosis-related lncRNAs

BACKGROUND

As a general female malignant tumor, Uterine Corpus Endometrial Carcinoma (UCEC) has high mortality and relapses. Cuproptosis was found to play an essential role in tumor by more and more researches. However, it is still unclear of the prognostic value and function of cuproptosis related Long non-coding RNA (lncRNA) in UCEC.

METHODS

Sequencing data with the corresponding clinical data and cuproptosis-related genes (CRGs) data were obtained from the Cancer Gene Atlas (TCGA) database and cuproptosis related studies. Pearson test was applied to select cuproptosis-related lncRNAs (CRLs). Prognosis associated CRLs was identified by univariate Cox analysis and the predictors were determined by least absolute shrinkage and selection operator (Lasso)-Cox and multivariate Cox analyses to construct the cuproptosis-related lncRNA prognostic signature (CRLPS). The performance of the CRLPs was evaluated by consistency index (C-index) and Kaplan-Meier analysis. A nomogram model was constructed for survival prediction and the accuracy of the model was evaluated by calibration curve. Finally, immune related analyses were applied to predict immune responses and identify drugs with potential efficacy for the overall survival (OS).

RESULTS

A total of 734 CRLs were found and 29 of them were identified as prognosis related lncRNAs. 12 CRLs were finally determined to build the CRLPS which revealed good ability on prognosis predicting. Subsequently, risk score of the CRLPS and grade were assessed as independent prognosis factors for UCEC, based on which the prognostic model provided the highest prediction accuracy of 99.7%. The calibration curve suggested that the prediction results consisted well with the observation. Enrichment analysis showed the CRLPS was mainly associated with tumor development and immune response. Patients in low tumor mutation burden (TMB) group had poorer OS. Significant difference was found in tumor immune dysfunction and exclusion (TIDE) score between different risk score groups. Finally, based on the CRLPs, drug sensitivity analysis identified nine anticancer drugs with potential efficacy on prognosis.

CONCLUSION

Cuproptosis-related lncRNA prognostic signature was constructed for UCEC for the first time. Its high reliability and accuracy on predicting prognosis and immunotherapy response provided new perspective to explore the tumor mechanism and improve clinical prognosis. Nine discovered sensitive drugs provided important clues for personalized treatment of UCEC.

Early Stage Finding of an Immune-Enhanced Genetic Subtype of Nonsmall Cell Lung Cancer Related with T-Cell Depletion

Background

Molecular categorization of lung cancer in medical care is becoming increasingly important on a regular basis. One of the molecular classifications of NSCLC (early-stage NSCLC) supports that tumors of different biological varieties differ in terms of their genomes and clinical characteristics.

Methods

Based on published immune cell signatures and early-stage NSCLC gene expression data including cancer genome maps, we used consensus cluster analysis to identify immune molecular subtypes associated with early-stage NSCLC expression subtypes. These subtypes were correlated with early-stage NSCLC expression subtypes. Next, applying a wide range of statistical techniques, we evaluated the link between molecular subtypes and clinical features, immunological microenvironment, and immunotherapy reactivity. Molecular subtypes were defined as a classification of cancerous cells.

Results

Multiple RNAseq cross-platform datasets of immune genes were used to identify two molecular subtypes (C1 and C2) of NSCLC, with C1 showing a more favorable prognosis than C2. The results were validated on the CSE datasets. In terms of clinical characteristics, C2 subtype samples with a worse prognosis showed a more advanced tumor stage and higher mortality. C2 showed immuno-infiltrative characteristics but had depletion of T-cells. Biological functions such as EMT were enriched on C2. A low TIDE score in C1 indicated that C1 samples could benefit from taking immunotherapy. C2 were more susceptible to standard chemotherapeutic treatments such paclitaxel, cisplatin, sorafenib, crizotinib, and erlotinib.

Conclusion

According to our findings, early-stage NSCLC patients may benefit from receiving treatment with immune checkpoint blockade therapy.

Definition of immune molecular subtypes with distinct immune microenvironment, recurrence, and PANoptosis features to aid clinical therapeutic decision-making

Objective: Cervical cancer poses a remarkable health burden to females globally. Despite major advances in early detection and treatment modalities, some patients still relapse. The present study proposed a novel immune molecular classification that reflected distinct recurrent risk and therapeutic responses in cervical cancer.

Methods: We retrospectively collected two cervical cancer cohorts: TCGA and GSE44001. Consensus clustering approach was conducted based on expression profiling of recurrence- and immune-related genes. The abundance of immune cells was inferred via five algorithms. Immune functions and signatures were quantified through ssGSEA. Genetic mutations were analyzed by maftools package. Immunotherapeutic response was inferred via tumor mutation burden (TMB), Tumor Immune Dysfunction and Exclusion (TIDE), and Submap methods. Finally, we developed a LASSO model for recurrence prediction.

Results: Cervical cancer samples were categorized into two immune subtypes (IC1, and IC2). IC2 exhibited better disease free survival (DFS), increased immune cell infiltration within the immune microenvironment, higher expression of immune checkpoints, higher activity of immune-relevant pathways (APC co-inhibition and co-stimulation, inflammation-promoting, MHC class I, IFN response, leukocyte and stromal fractions, macrophage regulation, and TCR Shannon), and higher frequencies of genetic mutations. This molecular classification exhibited a remarkable difference with existing immune subtypes, with diverse PANoptosis (pyroptosis, apoptosis and necroptosis) features. Patients in IC2 were more likely to respond to immunotherapy and targeted, and chemotherapeutic agents. The immune subtype-relevant signature was quantified to predict patients’ recurrence risk.

Conclusion: Altogether, we developed an immune molecular classification, which can be utilized in clinical practice to aid decision-making on recurrence management.

Metabolic pathway-based molecular subtyping of colon cancer reveals clinical immunotherapy potential and prognosis

PURPOSE

Colon cancer presents challenges to clinical diagnosis and management due to its high heterogeneity. For more efficient and convenient diagnosis and treatment of colon cancer, we are committed to characterizing the molecular features of colon cancer by pioneering a classification system based on metabolic pathways.

METHODS

Based on the 113 metabolic pathways and genes collected in the previous stage, we scored and filtered the metabolic pathways of each sample in the training set by ssGSEA, and obtained 16 metabolic pathways related to colon cancer recurrence. In consistent clustering of training set samples with recurrence-related metabolic pathway scores, we identified two robust molecular subtypes of colon cancer (MC1 and MC2). Furthermore, we performed multi-angle analysis on the survival differences of subtypes, metabolic characteristics, clinical characteristics, functional enrichment, immune infiltration, differences with other subtypes, stemness indices, TIDE prediction, and drug sensitivity, and finally constructed colon cancer prognostic model.

RESULTS

The results showed that the MC1 subtype had a poor prognosis based on higher immune activity and immune checkpoint gene expression. The MC2 subtype is associated with high metabolic activity and low expression of immune checkpoint genes and a better prognosis. The MC2 subtype was more responsive to PD-L1 immunotherapy than the MC1 subclass. However, we did not observe significant differences in tumor mutational burden between the two.

CONCLUSION

Two molecular subtypes of colon cancer based on metabolic pathways have distinct immune signatures. Constructing prognostic models based on subtype differential genes provides valuable reference for personalized therapy targeting unique tumor metabolic signatures.

Molecular Typing of Gastric Cancer Based on Invasion-Related Genes and Prognosis-Related Features

Background

This study aimed to construct a prognostic stratification system for gastric cancer (GC) using tumour invasion-related genes to more accurately predict the clinical prognosis of GC.

Methodology

Tumour invasion-related genes were downloaded from CancerSEA, and their expression data in the TCGA-STAD dataset were used to cluster samples via non-negative matrix factorisation (NMF). Differentially expressed genes (DEGs) between subtypes were identified using the limma package. KEGG pathway and GO functional enrichment analyses were conducted using the WebGestaltR package (v0.4.2). The immune scores of molecular subtypes were evaluated using the R package ESTIMATE, MCPcounter and the ssGSEA function of the GSVA package. Univariate, multivariate and lasso regression analyses of DEGs were performed using the coxph function of the survival package and the glmnet package to construct a RiskScore model. The robustness of the model was validated using internal and external datasets, and a nomogram was constructed based on the model.

Results

Based on 97 tumour invasion-related genes, 353 GC samples from TCGA were categorised into two subtypes, thereby indicating the presence of inter-subtype differences in prognosis. A total of 569 DEGs were identified between the two subtypes; of which, four genes were selected to construct the risk model. This four-gene signature was robust and exhibited stable predictive performance in different platform datasets (GSE26942 and GSE66229), indicating that the established model performed better than other existing models.

Conclusion

A prognostic stratification system based on a four-gene signature was developed with a desirable area under the curve in the training and independent validation sets. Therefore, the use of this system as a molecular diagnostic test is recommended to assess the prognostic risk of patients with GC.

Developing an immune signature for triple-negative breast cancer to predict prognosis and immune checkpoint inhibitor response

Aim: We aimed to develop a new signature based on immune-related genes to predict prognosis and response to immune checkpoint inhibitors in patients with triple-negative breast cancer (TNBC). Materials & methods: Single-sample gene set enrichment was used to develop an immune-based prognostic signature (IPRS) for TNBC patients. We conducted multivariate Cox analysis to evaluate the prognosis value of the IPRS. Result: An IPRS based on 66 prognostic genes was developed. Multivariate Cox analysis indicated that the IPRS was an independent factor for prognosis. PD-1, PD-L1, PD-L2 and CTLA4 gene expression was higher in the low-risk group, suggesting IPRS could predict the response to immune checkpoint inhibitors. Conclusion: The IPRS might be a reliable signature to predict TNBC patients' prognosis and response to immune checkpoint inhibitors, but needs prospective validation.

Increased co-expression of MEST and BRCA1 is associated with worse prognosis and immune infiltration in ovarian cancer

OBJECTIVE

The crosstalk between tumor microenvironment (TME) and cancer cells plays a critical role in the occurrence and development of ovarian cancer. Imprinted gene MEST is a tumor-promoting factor that modulates several carcinogenic signaling pathways. This study aimed to investigate the expression pattern of MEST and its association with immune cell infiltration.

METHODS

The transcriptome data of The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) database was utilized, and the expression and immune characteristics of MEST were verified by immunohistochemistry of ovarian cancer specimens. Kaplan-Meier Plotter was used to assess the prognostic value in patients with ovarian cancer.

RESULTS

We found that high expression of MEST was associated with diminished immune cell infiltration and worse prognosis of ovarian cancer patients in independent cohorts. There was a positive correlation between MEST and BRCA1 expression. The MEST^highBRCA1^high ovarian cancer group was correlated with lower infiltration of CD4⁺ cells, CD57⁺ cells, CD68⁺ cells and MPO⁺ cells, had worse overall survival (OS) in TCGA (HR = 1.57, p = 0.0004) and GSE27651 (HR = 4.27, p = 0.0002) cohorts, and predicted poor progress free survival (PFS) in GSE9891 (HR = 1.76, p = 0.0098) and GSE15622 (HR = 4.80, p = 0.0121) cohorts. Moreover, the expression of PD-L1 predicted unfavorable OS (HR = 2.48, p = 0.0415) and PFS (HR = 2.36, p = 0.0215) in MEST^lowBRCA1^low ovarian cancer group in GSE9891 cohort.

CONCLUSION

These findings suggest that the co-expression of MEST and BRCA1 may be an ideal combination for predicting the prognosis and response to immunotherapy in patients with ovarian cancer.

Comprehensive analysis of prognostic immune-related genes and drug sensitivity in cervical cancer

Background

Cervical cancer (CC) is the leading cause of cancer-related death in women. A limited number of studies have investigated whether immune-prognostic features can be used to predict the prognosis of CC. This study aimed to develop an improved prognostic risk scoring model (PRSM) for CC based on immune-related genes (IRGs) to predict survival and determine the key prognostic IRGs.

Methods

We downloaded the gene expression profiles and clinical data of CC patients from the TCGA and GEO databases. The ESTIMATE algorithm was used to calculate the score for both immune and stromal cells. Differentially expressed genes (DEGs) in different subpopulations were analyzed by “Limma”. A weighted gene co-expression network analysis (WGCNA) was used to establish a DEG co-expression module related to the immune score. Immune-related gene pairs (IRGPs) were constructed, and univariate- and Lasso-Cox regression analyses were used to analyze prognosis and establish a PRSM. A log-rank test was used to verify the accuracy and consistency of the scoring model. Identification of the predicted key IRG was ensured by the application of functional enrichment, DisNor, protein–protein interactions (PPIs) and heatmap. Finally, we extracted the key prognostic immune-related genes from the gene expression data, validated the key genes by immunohistochemistry and analyzed the correlation between their expression and drug sensitivity.

Results

A new PRSM was developed based on 22 IRGPs. The prognosis of the low-risk group in the model group (P < 0.001) and validation group (P = 0.039) was significantly better than that in the high-risk group. Furthermore, M1 and M2 macrophages were highly expressed in the low-risk group. Retinoic acid-inducible gene-I (RIG-I)-like receptors (RLRs) and the Janus kinase-signal transducer and activator of transcription (JAK-STAT) signaling pathway were significantly enriched in the low-risk group. Three representative genes (CD80, CD28, and LCP2) were markers of CC prognosis. CD80 and CD28 may more prominent represent important indicators to improve patient prognosis. These key genes was positively correlated with drug sensitivity. Finally, we found that differences in the sensitivity to JNK inhibitors could be distinguished based on the use and risk grouping of this PRSM.

Conclusions

The prognostic model based on the IRGs and key genes have potential clinical significance for predicting the prognosis of CC patients, providing a foundation for clinical prognosis judgment and individualized treatment.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12935-021-02333-9.

Identification of metabolism-associated molecular subtype in ovarian cancer

Ovarian cancer (OC) is the most lethal gynaecological cancer with genomic complexity and extensive heterogeneity. This study aimed to characterize the molecular features of OC based on the gene expression profile of 2752 previously characterized metabolism-relevant genes and provide new strategies to improve the clinical status of patients with OC. Finally, three molecular subtypes (C1, C2 and C3) were identified. The C2 subtype displayed the worst prognosis, upregulated immune-cell infiltration status and expression level of immune checkpoint genes, lower burden of copy number gains and losses and suboptimal response to targeted drug bevacizumab. The C1 subtype showed downregulated immune-cell infiltration status and expression level of immune checkpoint genes, the lowest incidence of BRCA mutation and optimal response to targeted drug bevacizumab. The C3 subtype had an intermediate immune status, the highest incidence of BRCA mutation and a secondary optimal response to bevacizumab. Gene signatures of C1 and C2 subtypes with an opposite expression level were mainly enriched in proteolysis and immune-related biological process. The C3 subtype was mainly enriched in the T cell-related biological process. The prognostic and immune status of subtypes were validated in the Gene Expression Omnibus (GEO) dataset, which was predicted with a 45-gene classifier. These findings might improve the understanding of the diversity and therapeutic strategies for OC.

Identification of three immune molecular subtypes associated with immune profiles, immune checkpoints, and clinical outcome in multiple myeloma

Purpose

To identify the immune molecular subtype for MM to help achieve individualized and precise targeted therapy.

Methods

The GDC API was used to download the TCGA‐MM profile dataset, which contains 859 samples in total, all of which were anterior to the standard treatment after diagnosis. Moreover, 282, 298, and 258 samples were stage I, stage II, and stage III separately. We used the immune gene expression profile for consistent clustering; and used the R software package ConsensusClusterPlus to sort the immune molecular subtypes. Correlation between subtypes and clinical features, immunity, and prognosis was then analyzed.

Results

A total of 859 tumor samples were separated into these three subtypes, which were not meaningfully related to age or sex but showed a remarkable association with stage. The results suggested that obvious differences in immune metagene expression and expression of 10 immune checkpoint genes appeared among the three subtypes.

Conclusion

The three subtypes are distinctly different in terms of immune metagenes, immune checkpoint molecules, and clinical prognosis. The discovery of the immune microenvironment of MM could further reveal the strategy for immunotherapy in MM and provide a promising candidate prognostic tool for survival.

Significance of N6-Methyladenosine RNA Methylation Regulators in Immune Infiltrates of Ovarian Cancer

N6-methyladenosine (m6A) RNA methylation regulators play an important role in the occurrence and development of tumors. Here, we aimed to identify the potential roles of m6A RNA methylation regulators in immune infiltrates of ovarian cancer. We obtained two distinct m6A patterns (m6Acluster.A and m6Acluster.B) based on the expression levels of all 21 m6A RNA methylation regulators from The Cancer Genome Atlas (TCGA) database using a consensus clustering algorithm. Differential analysis of m6Acluster.A and m6Acluster.B identified 196 m6A-related genes. We further validated the m6A regulation mechanism based on the 196 m6A-related genes using another consensus clustering algorithm. Considering individual differences, principal component analysis algorithms were used to calculate an m6A score for each sample in order to quantify the m6A patterns. A low m6A score was associated with immune activation and enhanced response to immune checkpoint inhibitors, whereas a high m6A score was associated with tumor progression. Finally, we successfully verified the correlation between m6A regulators and immune microenvironment in OC using our microarray analysis data. In summary, m6A regulators play non-negligible roles in immune infiltrates of ovarian cancer. Our investigation of m6A patterns may help to guide future immunotherapy strategies for advanced ovarian cancer.

Integrative analysis of immune molecular subtypes and microenvironment characteristics of bladder cancer

The emergence of immunotherapy has provided an option of treatment methods for bladder cancer (BC). However, the beneficiaries of immunotherapy are still limited to small‐scale patients, and immunotherapy‐related adverse events often occur. It is a major challenge for clinical work to study the immune subtypes of BC and the molecular mechanism of immune escape, and identify the immune responders accurately. Here, we explore the immune molecular subtypes of bladder cancer and potential escape mechanisms. First, we screened the expression profiles of 303 differentially expressed immune‐related genes in BC patients from the Cancer Genome Atlas (TCGA) database, and successfully identified 4 molecular subtypes of BC. By comparing the clinical characteristics, immune cells infiltration, the expression of checkpoint genes, human leukocyte antigen (HLA) genes, and gene mutation status of different subtypes, we identified different clinical and immunological characteristics of 4 subtypes. Among 4 subtypes, Cluster 2 met the general characteristics of immunotherapy responders and responded well to immunotherapy, while Cluster 4 had the highest expression of immune characteristics, and is similar to the immune environment of normal bladder tissue. Then, the weighted gene co‐expression network analysis (WGCNA) of immune‐related genes revealed that brown module was positively correlated with subtypes. Pathway enrichment analysis explored the major pathways associated with subtypes, which are also associated with immune escape mechanisms. Moreover, the decision tree model, which was constructed by the principle of random forest screening factors, was also validated in internal validation set and external validation set from the Gene Expression Omnibus (GEO) cohort (GSE133624), and could achieve accurate subtypes prediction for BC patients with high‐throughput sequencing. Taken together, we explored the immune molecular subtypes and their mechanisms of BC, and these results may provide guidance for the development of new BC immunotherapy strategies.

Shaping Immune Responses in the Tumor Microenvironment of Ovarian Cancer

Reciprocal signaling between immune cells and ovarian cancer cells in the tumor microenvironment can alter immune responses and regulate disease progression. These signaling events are regulated by multiple factors, including genetic and epigenetic alterations in both the ovarian cancer cells and immune cells, as well as cytokine pathways. Multiple immune cell types are recruited to the ovarian cancer tumor microenvironment, and new insights about the complexity of their interactions have emerged in recent years. The growing understanding of immune cell function in the ovarian cancer tumor microenvironment has important implications for biomarker discovery and therapeutic development. This review aims to describe the factors that shape the phenotypes of immune cells in the tumor microenvironment of ovarian cancer and how these changes impact disease progression and therapy.

Bioinformatics Analysis of GFAP as a Potential Key Regulator in Different Immune Phenotypes of Prostate Cancer

Tumor immune escape plays an essential role in both cancer progression and immunotherapy responses. For prostate cancer (PC), however, the molecular mechanisms that drive its different immune phenotypes have yet to be fully elucidated. Patient gene expression data were analyzed from The Cancer Genome Atlas-prostate adenocarcinoma (TCGA-PRAD) and the International Cancer Genome Consortium (ICGC) databases. We used a Cell-type Identification by Estimating Relative Subsets of RNA Transcripts (CIBERSORT) analysis and an unsupervised clustering analysis to identify patient subgroups with distinct immune phenotypes. These distinct phenotypes were then explored for associations for differentially expressed genes (DEGs) and both epigenetic and genetic landscapes. Finally, we used a protein-protein interaction analysis to identify key hub genes. We identified two patient subgroups with independent immune phenotypes associated with the expression of Programmed death-ligand 1 (PD-L1). Patient samples in Cluster 1 (C1) had higher scores for immune-cell subsets compared to Cluster 2 (C2), and C2 samples had higher specific somatic mutations, MHC mutations, and genomic copy number variations compared to C1. We also found additional cluster phenotype differences for DNA methylation, microRNA (miRNA) expression, and long noncoding RNA (lncRNA) expression. Furthermore, we established a 4-gene model to distinguish between clusters by integrating analyses for DEGs, lncRNAs, miRNAs, and methylation. Notably, we found that glial fibrillary acidic protein (GFAP) might serve as a key hub gene within the genetic and epigenetic regulatory networks. These results improve our understanding of the molecular mechanisms underlying tumor immune phenotypes that are associated with tumor immune escape. In addition, GFAP may be a potential biomarker for both PC diagnosis and prognosis.

Development and Validation of a Prognostic Signature Based on Immune Genes in Cervical Cancer

BACKGROUND

Cervical cancer is one of the most common types of gynecological malignancies worldwide. This study aims to develop an immune signature to predict survival in cervical cancer.

METHOD

The gene expression data of 296 patients with cervical cancer from The Cancer Genome Atlas database (TCGA) and immune-related genes from the Immunology Database and Analysis Portal (ImmPort) database were included in this study. The immune signature was developed based on prognostic genes. The validation dataset was downloaded from the Gene Expression Omnibus (GEO) database.

RESULT

The immune signature namely immune-based prognostic score (IPRS) was developed with 229 genes. Multivariate analysis revealed that the IPRS was an independent prognostic factor for overall survival (OS) and progression-free survival (PFS) in patients with cervical cancer. Patients were stratified into high IPRS and low IPRS groups, and those in the high IPRS group were associated with better survival, which was validated in the validation set. A nomogram with IPRS and stage was constructed to predict mortality in cervical cancer.

CONCLUSIONS

We developed a robust prognostic signature IPRS that could be used to predict patients' survival outcome.

Reclassification of Kidney Clear Cell Carcinoma Based on Immune Cell Gene-Related DNA CpG Pairs

Background: A new method was developed based on the relative ranking of gene expression level, overcoming the flaw of the batch effect, and having reliable results in various studies. In the current study, we defined the two methylation sites as a pair. The methylation level in a specific sample was subject to pairwise comparison to calculate a score for each CpGs-pair. The score was defined as a CpGs-pair score. If the first immune-related CpG value was higher than the second one in a specific CpGs-pair, the output score of this immune-related CpGs-pair was 1; otherwise, the output score was 0. This study aimed to construct a new classification of Kidney Clear Cell Carcinoma (KIRC) based on DNA CpGs (methylation sites) pairs. Methods: In this study, the biomarkers of 28 kinds of immune infiltration cells and corresponding methylation sites were acquired. The methylation data were compared between KIRC and normal tissue samples, and differentially methylated sites (DMSs) were obtained. Then, DNA CpGs-pairs were obtained according to the pairs of DMSs. In total, 441 DNA CpGs-pairs were utilized to construct a classification using unsupervised clustering analysis. We also analyzed the potential mechanism and therapy of different subtypes, and validated them in a testing set. Results: The classification of KIRC contained three subgroups. The clinicopathological features were different across three subgroups. The distribution of immune cells, immune checkpoints and immune-related mechanisms were significantly different across the three clusters. The mutation and copy number variation (CNV) were also different. The clinicopathological features and potential mechanism in the testing dataset were consistent with those in the training set. Conclusions: Our findings provide a new accurate and stable classification for developing personalized treatments for the new specific subtypes.

Tumor Immune Microenvironment Related Gene-Based Model to Predict Prognosis and Response to Compounds in Ovarian Cancer

BACKGROUND

The tumor immune microenvironment (TIME) has been recognized to be an imperative factor facilitating the acquisition of many cancer-related hallmarks and is a critical target for targeted biological therapy. This research intended to construct a risk score model premised on TIME-associated genes for prediction of survival and identification of potential drugs for ovarian cancer (OC) patients.

METHODS AND RESULTS

The stromal and immune scores were computed utilizing the ESTIMATE algorithm in OC patient samples from The Cancer Genome Atlas (TCGA) database. Weighted gene co-expression network and differentially expressed genes analyses were utilized to detect stromal-and immune-related genes. The Least Absolute Shrinkage and Selection Operator (LASSO)-Cox regression was utilized for additional gene selection. The genes that were selected were utilized as the input for a stepwise regression to construct a TIME-related risk score (TIMErisk), which was then validated in Gene Expression Omnibus (GEO) database. For the evaluation of the protein expression levels of TIME regulators, the Human Protein Atlas (HPA) dataset was utilized, and for their biological functions, the TIMER and CIBERSORT algorithm, immunoreactivity, and Immune Cell Abundance Identifier (ImmuCellAI) were used. Possible OC medications were forecasted utilizing the Genomics of Drug Sensitivity in Cancer (GDSC) database and connectivity map (CMap). TIMErisk was developed based on ALPK2, CPA3, PTGER3, CTHRC1, PLA2G2D, CXCL11, and ZNF683. High TIMErisk was recognized as a poor factor for survival in the GEO and TCGA databases; subgroup analysis with FIGO stage, grade, lymphatic and venous invasion, debulking, and tumor site also indicated similar results. Functional immune cells corresponded to more incisive immune reactions, including secretion of chemokines and interleukins, natural killer cell cytotoxicity, TNF signaling pathway, and infiltration of activated NK cells, eosinophils, and neutrophils in patients with low TIMErisk. Several small molecular medications which may enhance the prognosis of patients in the TIMErisk subgroup were identified. Lastly, an enhanced predictive performance nomogram was constructed by compounding TIMErisk with the FIGO stage and debulking.

CONCLUSION

These findings may offer a valuable indicator for clinical stratification management and personalized therapeutic options for OC patients and may be a foundation for future mechanistic research of their association.

Classification of ovarian cancer associated with BRCA1 mutations, immune checkpoints, and tumor microenvironment based on immunogenomic profiling

Background

Ovarian cancer is a highly fatal gynecological malignancy and new, more effective treatments are needed. Immunotherapy is gaining attention from researchers worldwide, although it has not proven to be consistently effective in the treatment of ovarian cancer. We studied the immune landscape of ovarian cancer patients to improve the efficacy of immunotherapy as a treatment option.

Methods

We obtained expression profiles, somatic mutation data, and clinical information from The Cancer Genome Atlas. Ovarian cancer was classified based on 29 immune-associated gene sets, which represented different immune cell types, functions, and pathways. Single-sample gene set enrichment (ssGSEA) was used to quantify the activity or enrichment levels of the gene sets in ovarian cancer, and the unsupervised machine learning method was used sort the classifications. Our classifications were validated using Gene Expression Omnibus datasets.

Results

We divided ovarian cancer into three subtypes according to the ssGSEA score: subtype 1 (low immunity), subtype 2 (median immunity), and subtype 3 (high immunity). Most tumor-infiltrating immune cells and immune checkpoint molecules were upgraded in subtype 3 compared with those in the other subtypes. The tumor mutation burden (TMB) was not significantly different among the three subtypes. However, patients with BRCA1 mutations were consistently detected in subtype 3. Furthermore, most immune signature pathways were hyperactivated in subtype 3, including T and B cell receptor signaling pathways, PD-L1 expression and PD-1 checkpoint pathway the NF-κB signaling pathway, Th17 cell differentiation and interleukin-17 signaling pathways, and the TNF signaling pathway.

Conclusion

Ovarian cancer subtypes that are based on immune biosignatures may contribute to the development of novel therapeutic treatment strategies for ovarian cancer.

A Methylation-Based Reclassification of Bladder Cancer Based on Immune Cell Genes

Simple Summary

Bladder cancer (BC) development is highly related to immune cell infiltration. In this study, we aimed to construct a new classification of bladder cancer molecular subtypes based on immune-cell-associated CpG(Methylation) sites. The classification was accurate and stable. BC patients were successfully divided into three subtypes based on the immune-cell-associated CpG sites. The clinicopathologic features, distribution of immune cells, level of expression of checkpoints, stromal score, immune score, ESTIMATEScore, tumor purity, APC co_inhibition, APC co_stimulation, HLA, MHC class_I, Type I IFN_respons, Type II IFN response, and DNA stemness score (DNAss) presented significant differences among the three subgroups. The specific genomic alteration was also different across subgroups. High-level immune infiltration showed a correlation with high-level methylation. A lower RNA stemness score (RNAss) was associated with higher immune infiltration. Cluster 2 demonstrated a better response to chemotherapy. The anti-cancer targeted drug therapy results are different among the three subgroups.

Abstract

Background: Bladder cancer is highly related to immune cell infiltration. This study aimed to develop a new classification of BC molecular subtypes based on immune-cell-associated CpG sites. Methods: The genes of 28 types of immune cells were obtained from previous studies. Then, methylation sites corresponding to immune-cell-associated genes were acquired. Differentially methylated sites (DMSs) were identified between normal samples and bladder cancer samples. Unsupervised clustering analysis of differentially methylated sites was performed to divide the sites into several subtypes. Then, the potential mechanism of different subtypes was explored. Results: Bladder cancer patients were divided into three groups. The cluster 3 subtype had the best prognosis. Cluster 1 had the poorest prognosis. The distribution of immune cells, level of expression of checkpoints, stromal score, immune score, ESTIMATEScore, tumor purity, APC co_inhibition, APC co_stimulation, HLA, MHC class_I, Type I IFN Response, Type II IFN Response, and DNAss presented significant differences among the three subgroups. The distribution of genomic alterations was also different. Conclusions: The proposed classification was accurate and stable. BC patients could be divided into three subtypes based on the immune-cell-associated CpG sites. Specific biological signaling pathways, immune mechanisms, and genomic alterations were varied among the three subgroups. High-level immune infiltration was correlated with high-level methylation. The lower RNAss was associated with higher immune infiltration. The study of the intratumoral immune microenvironment may provide a new perspective for BC therapy.

Interactions between tumor mutation burden and immune infiltration in ovarian cancer

The aim of the study was to evaluate the relationship between tumor mutational burden (TMB) and immune infiltration in ovarian cancer. We extracted somatic mutational data and gene expression profiles of ovarian cancer from The Cancer Genome Atlas (TCGA). The samples were separated into low and high TMB groups. Correlations between TMB and cancer prognosis were analyzed and immune cell infiltration in the high and low TMB subgroups was calculated using the CIBERSORT package software. High TMB was significantly related to an improved survival rate. We identified 4 TMB-related core genes that were significantly associated with prognosis. Furthermore, mutations in the 4 genes were associated with immune cell infiltration. We also found a high proportion of naive B cells and activated NK cells in the high TMB group, while increased proportions of memory B cells and plasma cells were found in the low TMB group. Overall, our study indicated that patients with a higher TMB level experienced a favorable survival outcome and this may influence immune infiltration in ovarian cancer. Furthermore, the 4 TMB-related core genes were highly correlated with prognosis and the level of immune cell infiltration.

Identification of immune-enhanced molecular subtype associated with BRCA1 mutations, immune checkpoints and clinical outcome in ovarian carcinoma

Ovarian carcinoma has the highest mortality among the malignant tumours in gynaecology, and new treatment strategies are urgently needed to improve the clinical status of ovarian carcinoma patients. The Cancer Genome Atlas (TCGA) cohort were performed to explore the immune function of the internal environment of tumours and its clinical correlation with ovarian carcinoma. Finally, four molecular subtypes were obtained based on the global immune‐related genes. The correlation analysis and clinical characteristics showed that four subtypes were all significantly related to clinical stage; the immune scoring results indicated that most immune signatures were upregulated in C3 subtype, and the majority of tumour‐infiltrating immune cells were upregulated in both C3 and C4 subtypes. Compared with other subtypes, C3 subtype had a higher BRCA1 mutation, higher expression of immune checkpoints, and optimal survival prognosis. These findings of the immunological microenvironment in tumours may provide new ideas for developing immunotherapeutic strategies for ovarian carcinoma.