Pan-cancer proteogenomics characterization of tumor immunity.

Despite the successes of immunotherapy in cancer treatment over recent decades, less than <10%-20% cancer cases have demonstrated durable responses from immune checkpoint blockade. To enhance the efficacy of immunotherapies, combination therapies suppressing multiple immune evasion mechanisms are increasingly contemplated. To better understand immune cell surveillance and diverse immune evasion responses in tumor tissues, we comprehensively characterized the immune landscape of more than 1,000 tumors across ten different cancers using CPTAC pan-cancer proteogenomic data. We identified seven distinct immune subtypes based on integrative learning of cell type compositions and pathway activities. We then thoroughly categorized unique genomic, epigenetic, transcriptomic, and proteomic changes associated with each subtype. Further leveraging the deep phosphoproteomic data, we studied kinase activities in different immune subtypes, which revealed potential subtype-specific therapeutic targets. Insights from this work will facilitate the development of future immunotherapy strategies and enhance precision targeting with existing agents.

Proteomic analysis of breast cancer based on immune subtypes.

BACKGROUND: Immunotherapy is applied to breast cancer to resolve the limitations of survival gain in existing treatment modalities. With immunotherapy, a tumor can be classified into immune-inflamed, excluded and desert based on the distribution of immune cells. We assessed the clinicopathological features, each subtype's prognostic value and differentially expressed proteins between immune subtypes. METHODS: Immune subtyping and proteomic analysis were performed on 56 breast cancer cases with neoadjuvant chemotherapy. The immune subtyping was based on the level of tumor-infiltrating lymphocytes (TILs) and Klintrup criteria. If the level of TILs was ≥ 10%, it was classified as immune-inflamed type without consideration of the Klintrup criteria. In cases of 1-9% TIL, Klintrup criteria 1-3 were classified as the immune-excluded subtype and Klintrup criteria not available (NA) was classified as NA. Cases of 1% TILs and Klintrup 0 were classified as the immune-desert subtype. Mass spectrometry was used to identify differentially expressed proteins in formalin-fixed paraffin-embedded biopsy tissues. RESULTS: Of the 56 cases, 31 (55%) were immune-inflamed, 21 (38%) were immune-excluded, 2 (4%) were immune-desert and 2 (4%) were NA. Welch's t-test revealed two differentially expressed proteins between immune-inflamed and immune-excluded/desert subtypes. Coronin-1A was upregulated in immune-inflamed tumors (adjusted p = 0.008) and α-1-antitrypsin was upregulated in immune-excluded/desert tumors (adjusted p = 0.008). Titin was upregulated in pathologic complete response (pCR) than non-pCR among immune-inflamed tumors (adjusted p = 0.036). CONCLUSIONS: Coronin-1A and α-1-antitrypsin were upregulated in immune-inflamed and immune-excluded/desert subtypes, respectively. Titin's elevated expression in pCR within the immune-inflamed subtype may indicate a favorable prognosis. Further studies involving large representative cohorts are necessary to validate these findings.

Comprehensive analysis and identification of subtypes and hub genes of high immune response in lung adenocarcinoma.

BACKGROUND: The advent of immunotherapy targeting immune checkpoints has conferred significant clinical advantages to patients with lung adenocarcinoma (LUAD); However, only a limited subset of patients exhibit responsiveness to this treatment. Consequently, there is an imperative need to stratify LUAD patients based on their response to immunotherapy and enhance the therapeutic efficacy of these treatments. METHODS: The differentially co-expressed genes associated with CD8 + T cells were identified through weighted gene co-expression network analysis (WGCNA) and the Search Tool for the Retrieval of Interacting Genes (STRING) database. These gene signatures facilitated consensus clustering for TCGA-LUAD and GEO cohorts, categorizing them into distinct immune subtypes (C1, C2, C3, and C4). The Tumor Immune Dysfunction and Exclusion (TIDE) model and Immunophenoscore (IPS) analysis were employed to assess the immunotherapy response of these subtypes. Additionally, the impact of inhibitors targeting five hub genes on the interaction between CD8 + T cells and LUAD cells was evaluated using CCK8 and EDU assays. To ascertain the effects of these inhibitors on immune checkpoint genes and the cytotoxicity mediated by CD8 + T cells, flow cytometry, qPCR, and ELISA methods were utilized. RESULTS: Among the identified immune subtypes, subtypes C1 and C3 were characterized by an abundance of immune components and enhanced immunogenicity. Notably, both C1 and C3 exhibited higher T cell dysfunction scores and elevated expression of immune checkpoint genes. Multi-cohort analysis of Lung Adenocarcinoma (LUAD) suggested that these subtypes might elicit superior responses to immunotherapy and chemotherapy. In vitro experiments involved co-culturing LUAD cells with CD8 + T cells and implementing the inhibition of five pivotal genes to assess their function. The inhibition of these genes mitigated the immunosuppression on CD8 + T cells, reduced the levels of PD1 and PD-L1, and promoted the secretion of IFN-γ and IL-2. CONCLUSIONS: Collectively, this study delineated LUAD into four distinct subtypes and identified five hub genes correlated with CD8 + T cell activity. It lays the groundwork for refining personalized therapy and immunotherapy strategies for patients with LUAD.

Multi-omics analysis of immune-related microbiome and prognostic model in head and neck squamous cell carcinoma.

OBJECTIVES: The aim of our study is to explore the transcriptional and microbial characteristics of head and neck cancer's immune phenotypes using a multi-omics approach. MATERIALS AND METHODS: Employing TCGA data, we analyzed head and neck squamous cell carcinoma (HNSCC) immune cells with CIBERSORT and identified differentially expressed genes using DESeq2. Microbial profiles, obtained from the TCMA database, were analyzed using LEfSe algorithm to identify differential microbes in immune cell infiltration (ICI) subgroups. Random Forest algorithm and deep neural network (DNN) were employed to select microbial features and developed a prognosis model. RESULTS: We categorized HNSCC into three immune subtypes, finding ICI-2 with the worst prognosis and distinct microbial diversity. Our immune-related microbiome (IRM) model outperformed the TNM staging model in predicting survival, linking higher IRM model scores with poorer prognosis, and demonstrating clinical utility over TNM staging. Patients categorized as low-risk by the IRM model showed higher sensitivity to cisplatin and sorafenib treatments. CONCLUSIONS: This study offers a comprehensive exploration of the ICI landscape in HNSCC. We provide a detailed scenario of immune regulation in HNSCC and report a correlation between differing ICI patterns, intratumor microbiome, and prognosis. This research aids in identifying prime candidates for optimizing treatment strategies in HNSCC. CLINICAL RELEVANCE: This study revealed the microbial signatures associated with immunophenotyping of HNSCC and further found the microbial signatures associated with prognosis. The prognostic model based on IRM microbes is helpful for early prediction of patient prognosis and assisting clinical decision-making.

A comprehensive analysis of immune features and construction of an immune gene diagnostic model for sepsis.

Sepsis is a life-threatening syndrome resulting from immune system dysfunction that is caused by infection. It is of great importance to analyze the immune characteristics of sepsis, identify the key immune system related genes, and construct diagnostic models for sepsis. In this study, the sepsis transcriptome and expression profiling data were merged into an integrated dataset containing 277 sepsis samples and 117 non-sepsis control samples. Single-sample gene set enrichment analysis (ssGSEA) was used to assess the immune cell infiltration. Two sepsis immune subtypes were identified based on the 22 differential immune cells between the sepsis and the healthy control groups. Weighted gene co-expression network analysis (WCGNA) was used to identify the key module genes. Then, 36 differentially expressed immune-related genes were identified, based on which a robust diagnostic model was constructed with 11 diagnostic genes. The expression of 11 diagnostic genes was finally assessed in the training and validation datasets respectively. In this study, we provide comprehensive insight into the immune features of sepsis and establish a robust diagnostic model for sepsis. These findings may provide new strategies for the early diagnosis of sepsis in the future.

Comprehensive analysis of immune subtypes reveals the prognostic value of cytotoxicity and FAP+ fibroblasts in stomach adenocarcinoma.

BACKGROUND: The heterogeneity limits the effective application of immune checkpoint inhibitors for patients with stomach adenocarcinoma (STAD). Precise immunotyping can help select people who may benefit from immunotherapy and guide postoperative management by describing the characteristics of tumor microenvironment. METHODS: Gene expression profiles and clinical information of patients were collected from ACRG and TCGA-STAD datasets. The immune subtypes (ISs) were identified by consensus clustering analysis. The tumor immune microenvironments (TIME) of each IS were characterized using a series of immunogenomics methods and further confirmed by multiplex immunohistochemistry (mIHC) staining in clinical samples. Two online datasets and one in-house dataset were utilized to construct and validate a prognostic immune-related gene (IRG) signature. RESULTS: STAD patients were stratified into five reproducible ISs. IS1 (immune deserve subtype) had low immune infiltration and the highest degree of HER2 gene mutation. With abundant CD8+ T cells infiltration and activated cytotoxicity reaction, patients in the IS2 (immune-activated subtype) had the best overall survival (OS). IS3 and IS4 subtypes were both in the reactive stroma state and indicated the worst prognosis. However, IS3 (immune-inhibited subtype) was characterized by enrichment of FAP+ fibroblasts and upregulated TGF-ß signaling pathway, while IS4 (activated stroma subtype) was characterized by enrichment of ACTA2+ fibroblasts. In addition, mIHC staining confirmed that TGF-ß upregulated FAP+ fibroblasts were independent risk factor of OS. IS5 (chronic inflammation subtype) displayed moderate immune cells infiltration and had a relatively good survival. Lastly, we developed a nine-IRG signature model with a robust performance on overall survival prognostication. CONCLUSIONS: The immunotyping is indicative for characterize the TIME heterogeneity and the prediction of tumor prognosis for STADs, which may provide valuable stratification for the design of future immunotherapy.

Comprehensive molecular analysis of a four-pyroptosis-gene signature with prognosis and immune landscape in lung adenocarcinoma.

Pyroptosis plays an important role in tumor immunity. However, the biological behavior and prognostic significance of pyroptosis remain unclear. We identified 41 pyroptosis regulators differently expressed in lung adenocarcinoma (LUAD). All cases of LUAD can be classified into two molecular subtypes using unsupervised clustering algorithm. Using multiple analyses, a four-pyroptosis-gene signature was constructed, and all LUAD patients were categorized as low-risk or high-risk with a longer overall survival (OS) time in the low-risk group(P < 0.001). This signature had power prognosis and stratification which was validated by six independent datasets and clinical subtypes. Besides, this signature showed distinct clinical outcomes, immune landscapes in different risk groups. Moreover, the low-risk group had a higher response against immunotherapy with a lower TIDE score. Importantly, this signature surpassed other biomarkers (TIDE, TMB, PD-L1) in predicting prognosis. Overall, the current study might help with precise prognostic prediction and crucial treatment strategies, eventually promoting tailored therapy for LUAD patients.

The immune subtypes and landscape of sarcomas.

BACKGROUND: Considering the molecular heterogeneity of sarcomas and their immunologically quiet character, immunotherapy (e.g., immune checkpoint inhibitors) plays a viable role in only a subset of these tumors. This study aimed to determine the immune subtypes (IMSs) of sarcomas for selecting suitable patients from an extremely heterogeneous population. RESULTS: By performing consensus clustering analysis of the gene expression profiles of 538 patients with sarcomas in online databases, we stratified sarcomas into three IMSs characterized by different immune cell features, tumor mutational burdens (TMBs), gene mutations, and clinical outcomes. IMS1 showed an immune "hot" and immunosuppressive phenotype, the highest frequencies of CSMD3 mutation but the lowest frequencies of HMCN1 and LAMA2 mutations; these patients had the worst progression-free survival (PFS). IMS2 was defined by a high TMB and more gene mutations, but had the lowest frequency of MND1 mutations. IMS3 displayed the highest MDN1 expression level and an immune "cold" phenotype, these patients had the worst PFS. Each subtype was associated with different expression levels of immunogenic cell death modulators and immune checkpoints. Moreover, we applied graph learning-based dimensionality reduction to the immune landscape and identified significant intra-cluster heterogeneity within each IMS. Finally, we developed and validated an immune gene signature with good prognostic performance. CONCLUSIONS: Our results provide a conceptual framework for understanding the immunological heterogeneity of sarcomas. The identification of immune-related subtypes may facilitate optimal selection of sarcoma patients who will respond to appropriate therapeutic strategies.

Identification of cervical squamous cell carcinoma feature genes and construction of a prognostic model based on immune-related features.

As heterogeneity of cervical squamous cell carcinoma (CSCC), prognosis assessment for CSCC patients remain challenging. To develop novel prognostic strategies for CSCC patients, associated biomarkers are urgently needed. This study aimed to cluster CSCC samples from a molecular perspective. CSCC expression data sets were obtained from The Cancer Genome Atlas and based on the accessed expression profile, a co-expression network was constructed with weighted gene co-expression network analysis to form different gene modules. Tumor microenvironment was evaluated using ESTIMATE algorithm, observing that the brown module was highly associated with tumor immunity. CSCC samples were clustered into three subtypes by consensus clustering based on gene expression profiles in the module. Gene set variation analysis showed differences in immune-related pathways among the three subtypes. CIBERSORT and single-sample gene set enrichment analysis analyses showed the difference in immune cell infiltration among subtype groups. Also, Human leukocyte antigen protein expression varied considerably among subtypes. Subsequently, univariate, Lasso and multivariate Cox regression analyses were performed on the genes in the brown module and an 8-gene prognostic model was constructed. Kaplan-Meier analysis illuminated that the low-risk group manifested a favorable prognosis, and receiver operating characteristic curve showed that the model has good predictive performance. qRT-PCR was used to examine the expression status of the prognosis-associated genes. In conclusion, this study identified three types of CSCC from a molecular perspective and established an effective prognostic model for CSCC, which will provide guidance for clinical subtype identification of CSCC and treatment of patients.

A pan-cancer study of spalt-like transcription factors 1/2/3/4 as therapeutic targets.

Spalt-like transcription factors (SALLs) are evolutionarily conserved proteins that participate in embryonic development. Four members of the SALL family, SALL1, SALL2, SALL3, and SALL4, are involved in cellular apoptosis, angiogenesis, invasion, and metastasis of tumors. We used the TCGA pan-cancer data to conduct a comprehensive analysis of SALL genes. High heterogeneity in the expression of these genes was observed across various cancers, SALL1 and SALL2 were downregulated, whereas SALL4 was upregulated. Moreover, we verified that SALL4 was commonly associated with survival disadvantage, whereas others were linked to a better prognosis. In renal cancer, SALL1, SALL2, and SALL3 showed downregulation, suggesting that they acted as tumor suppressors. Furthermore, SALLs were associated with immune infiltrate subtypes, with a close association between different degrees of infiltration of stromal cells and immune cells. DNA and RNA analyses in different tumors suggested different degrees of negative or positive correlation with tumor stem cell-like features. Finally, we revealed that SALLs were related to cancer cell resistance. Our results highlight the necessity to further study each SALL gene as a separate entity in specific types of cancer. Although this article showed that SALLs could be promising targets for cancer therapy, it needs further studies to validate the findings.

Multi-omics analysis of an immune-based prognostic predictor in non-small cell lung cancer.

BACKGROUND: Inhibitors targeting immune checkpoints, such as PD-1/PD-L1 and CTLA-4, have prolonged survival in small groups of non-small cell lung cancer (NSCLC) patients, but biomarkers predictive of the response to the immune checkpoint inhibitors (ICIs) remain rare. METHODS: The nonnegative matrix factorization (NMF) was performed for TCGA-NSCLC tumor samples based on the LM22 immune signature to construct subgroups. Characterization of NMF subgroups involved the single sample gene set variation analysis (ssGSVA), and mutation/copy number alteration and methylation analyses. Construction of RNA interaction network was based on the identification of differentially expressed RNAs (DERs). The prognostic predictor was constructed by a LASSO-Cox regression model. Four GEO datasets were used for the validation analysis. RESULTS: Four immune based NMF subgroups among NSCLC patients were identified. Genetic and epigenetic analyses between subgroups revealed an important role of somatic copy number alterations in determining the immune checkpoint expression on specific immune cells. Seven hub genes were recognized in the regulatory network closely related to the immune phenotype, and a three-gene prognosis predictor was constructed. CONCLUSIONS: Our study established an immune-based prognosis predictor, which might have the potential to select subgroups benefiting from the ICI treatment, for NSCLC patients using publicly available databases.

Identification of glioblastoma immune subtypes and immune landscape based on a large cohort.

Glioblastomas (GBM) are the most common primary brain malignancy and also the most aggressive one. In addition, GBM have to date poor treatment options. Therefore, understanding the GBM microenvironment may help to design immunotherapy treatments and rational combination strategies. In this study, the gene expression profiles and clinical follow-up data were downloaded from TCGA-GBM, and the molecular subtypes were identified using ConsensusClusterPlus. Univariate and multivariate Cox regression were used to evaluate the prognostic value of immune subtypes. The Graph Structure Learning method was used for dimension reduction to reveal the internal structure of the immune system. A Weighted Correlation Network Analysis (WGCNA) was used to identify immune-related gene modules. Four immune subtypes (IS1, IS2, IS3, IS4) with significant prognosis differences were obtained. Interestingly, IS4 had the highest mutation rate. We also found significant differences in the distribution of the four subtypes at immune checkpoints, molecular markers, and immune characteristics. WGCNA identified 11 co-expressed module genes, and there were significant differences among the four subtypes. Finally, CD1A, CD1E, and IL23R genes with significant prognostic significance were selected as the final feature genes in the brown module. Overall, this study provided a conceptual framework for understanding the tumor immune microenvironment of GBM.

Identification of immune subtypes associated with neutrophils in tuberculosis infection based on weighted gene co-expression network analysis.

Tuberculosis (TB) is caused by Mycobacterium tuberculosis and is a major global health concern. Neutrophils play a significant role in TB infection and patient outcomes. This study aimed to identify gene modules associated with neutrophil infiltration in TB samples using WGCNA. Gene ontology and enrichment analyses were performed, and a random forest model was constructed to identify differentially expressed genes. K-means clustering was used to classify samples into subtypes, and immune-related scores, PD-L1 expression, HLA expression, and gene enrichment analysis were evaluated. The blue module showed significant correlation with neutrophils and enrichment in immune-related processes. The model exhibited good classification performance, and subtype 1 demonstrated higher immune-related scores, PD-L1 expression, HLA class I molecule expression, and immune-related pathway enrichment. These findings enhance our understanding of TB pathogenesis and provide potential targets for diagnosis and treatment strategies.

Prediction Model for Immunotherapy Efficacy in Hepatocellular Carcinoma Based on Alternative Splicing Sequencing Data.

Background: Integrating immune checkpoint inhibitors with multi-target tyrosine kinase inhibitors presents an innovative and hopeful strategy in liver cancer treatment. Nonetheless, a degree of resistance to this treatment is noticeable in certain patients. Alternative splicing (AS) represents a common biological process that controls the variety of life functions via isoforms. Purpose: Investigating how gene AS affects the effectiveness of combined immunotherapy in treating hepatocellular carcinoma (HCC). Methods: Our retrospective examination focused on AS's effect on immune therapy effectiveness, utilizing accessible tissue sequencing and clinical records for HCC. For corroborating our results, we gathered samples of drug-resistant HCC tissue, nearby tissues, HCC tissue with high drug responsiveness, and healthy liver tissue from clinical studies. Results: The study revealed a link between the frequency of AS occurrences, the expression levels of programmed cell death 1 ligand 1, and the resistance to tumor medications. Our study detailed the AS occurrences in HCC, leading to the creation of a risk-assessment function and a predictive model using AS data. The results of our study revealed that the risk score effectively distinguished between various immune subtypes and the effectiveness of immune therapy. Additional examination of the chosen AS occurrences uncovered their effects on both the immune microenvironment and cellular immunity. Our investigation also delved into the regulatory framework of AS, uncovering the role of stringently controlled splicing factors in the emergence of tumors and the modulation of the body's immune response. Conclusions: Increased AS in HCC diminishes the efficacy of immunotherapy; conversely, more AS in peritumoral tissue elevates the likelihood of tumor immune evasion.

A novel risk score system based on immune subtypes for identifying optimal mRNA vaccination population in hepatocellular carcinoma.

PURPOSE: Although mRNA vaccines have shown certain clinical benefits in multiple malignancies, their therapeutic efficacies against hepatocellular carcinoma (HCC) remains uncertain. This study focused on establishing a novel risk score system based on immune subtypes so as to identify optimal HCC mRNA vaccination population. METHODS: GEPIA, cBioPortal and TIMER databases were utilized to identify candidate genes for mRNA vaccination in HCC. Subsequently, immune subtypes were constructed based on the candidate genes. According to the differential expressed genes among various immune subtypes, a risk score system was established using machine learning algorithm. Besides, multi-color immunofluorescence of tumor tissues from 72 HCC patients were applied to validate the feasibility and efficiency of the risk score system. RESULTS: Twelve overexpressed and mutated genes associated with poor survival and APCs infiltration were identified as potential candidate targets for mRNA vaccination. Three immune subtypes (e.g. IS1, IS2 and IS3) with distinct clinicopathological and molecular profiles were constructed according to the 12 candidate genes. Based on the immune subtype, a risk score system was developed, and according to the risk score from low to high, HCC patients were classified into four subgroups on average (e.g. RS1, RS2, RS3 and RS4). RS4 mainly overlapped with IS3, RS1 with IS2, and RS2+RS3 with IS1. ROC analysis also suggested the significant capacity of the risk score to distinguish between the three immune subtypes. Higher risk score exhibited robustly predictive ability for worse survival, which was further independently proved by multi-color immunofluorescence of HCC samples. Notably, RS4 tumors exhibited an increased immunosuppressive phenotype, higher expression of the twelve potential candidate targets and increased genome altered fraction, and therefore might benefit more from vaccination. CONCLUSIONS: This novel risk score system based on immune subtypes enabled the identification of RS4 tumor that, due to its highly immunosuppressive microenvironment, may benefit from HCC mRNA vaccination.

A pan-cancer analysis unveiling the function of NR4A family genes in tumor immune microenvironment, prognosis, and drug response.

BACKGROUND: NR4A family genes play crucial roles in cancers. However, the role of NR4A family genes in cancers remains paradoxical as they promote or suppress tumorigenesis. OBJECTIVE: We aimed to conduct comprehensive analyses of the association between the expression of NR4A family genes and tumor microenvironment (TME) based on bioinformatics methods. METHODS: We collected RNA-seq data from 33 cancer types and 20 normal tissue sites from the TCGA and GTEx databases. Expression patterns of NR4A family genes and their associations with DNA methylation, miRNA, overall survival, drug responses, and tumor microenvironment were investigated. RESULTS: Significant downregulation of all NR4A family genes was observed in 15 cancer types. DNA promoter methylation and expression of NR4A family genes were negatively correlated in five cancers. The expression of 10 miRNAs targeting NR4A family genes was negatively correlated with the expression of NR4A family genes. High expression of all NR4A family genes was associated with poor prognosis in stomach adenocarcinoma and increased expressions of NR4A2 and NR4A3 were associated with poor prognosis in adrenocortical carcinoma. In addition, we found an elevated expression of NR4A2, which enhances the response to various chemotherapeutic drugs, whereas NR4A3 decreases drug sensitivity. Interestingly, in breast cancer, NR4A3 was significantly associated with C2 (IFN-γ dominant), C3 (inflammatory), and C6 (TGF-ß dominant) immune subtypes and infiltrated immune cell types, implying both oncogenic and tumor-suppressive functions of NR4A3 in breast cancer. CONCLUSION: The NR4A family genes have the potential to serve as a diagnostic, prognostic, and immunological marker of human cancers.

Circulating tumor DNA-based stratification strategy for chemotherapy plus PD-1 inhibitor in advanced non-small-cell lung cancer.

Stratification strategies for chemotherapy plus PD-1 inhibitors in advanced non-small-cell lung cancer (NSCLC) are critically demanded. We performed high-throughput panel-based deep next-generation sequencing and low-pass whole genome sequencing on prospectively collected circulating tumor DNA (ctDNA) specimens from 460 patients in the phase 3 CHOICE-01 study at different time points. We identified predictive markers for chemotherapy plus PD-1 inhibitor, including ctDNA status and genomic features such as blood-based tumor mutational burden, intratumor heterogeneity, and chromosomal instability. Furthermore, we established an integrated ctDNA-based stratification strategy, blood-based genomic immune subtypes (bGIS) scheme, to distinguish patients who benefit from the addition of PD-1 inhibitor to first-line chemotherapy. Moreover, we demonstrated potential applications for the dynamic monitoring of ctDNA. Overall, we proposed a potential therapeutic algorithm based on the ctDNA-based stratification strategy, shedding light on the individualized management of immune-chemotherapies for patients with advanced NSCLC.

Imaging phenotypes from MRI for the prediction of glioma immune subtypes from RNA sequencing: A multicenter study.

Tumor subtyping based on its immune landscape may guide precision immunotherapy. The aims of this study were to identify immune subtypes of adult diffuse gliomas with RNA sequencing data, and to noninvasively predict this subtype using a biologically interpretable radiomic signature from MRI. A subtype discovery dataset (n = 210) from a public database and two radiogenomic datasets (n = 130 and 55, respectively) from two local hospitals were included. Brain tumor microenvironment-specific signatures were constructed from RNA sequencing to identify the immune types. A radiomic signature was built from MRI to predict the identified immune subtypes. The pathways underlying the radiomic signature were identified to annotate their biological meanings. The reproducibility of the findings was verified externally in multicenter datasets. Three distinctive immune subtypes were identified, including an inflamed subtype marked by elevated hypoxia-induced immunosuppression, a "cold" subtype that exhibited scarce immune infiltration with downregulated antigen presentation, and an intermediate subtype that showed medium immune infiltration. A 10-feature radiomic signature was developed to predict immune subtypes, achieving an AUC of 0.924 in the validation dataset. The radiomic features correlated with biological functions underpinning immune suppression, which substantiated the hypothesis that molecular changes can be reflected by radiomic features. The immune subtypes, predictive radiomic signature, and radiomics-correlated biological pathways were validated externally. Our data suggest that adult-type diffuse gliomas harbor three distinctive immune subtypes that can be predicted by MRI radiomic features with clear biological significance. The immune subtypes, radiomic signature, and radiogenomic links can be replicated externally.

Corrigendum: Immunogenomic landscape in breast cancer reveals immunotherapeutically relevant gene signatures.

[This corrects the article DOI: 10.3389/fimmu.2022.805184.].

Metabolic-related gene signatures for survival prediction and immune cell subtypes associated with prognosis in intrahepatic cholangiocarcinoma.

OBJECTIVES: Our study aimed to reveal the metabolic-related gene signatures for survival prediction and immune cell subtypes associated with IHCC prognosis. METHODS: Differentially expressed metabolic genes were identified between survival group and dead group which were divided according to survival at discharge. Recursive feature elimination (RFE) and randomForest (RF) algorithms were applied to optimize the combination of feature metabolic genes, which were used to generate SVM classifier. Performance of SVM classifier was evaluated by receiver operating characteristic (ROC) curves. Gene set enrichment analysis (GSEA) was conducted to uncover the activated pathways in high risk group, and differences in immune cell distributions were revealed. RESULTS: There were 143 differentially expressed metabolic gens. RFE and RF identified 21 overlapping differentially expressed metabolic genes, and the constructed SVM classifier had excellent accuracy in training and validation dataset. RS survival prediction model was consisted of 10 metabolic genes. RS model had reliable predictive capability in the training and validation dataset. GSEA revealed 15 significant KEGG pathways that were relatively activated in the high risk group. High risk group had obviously lower counts of B cell naive and T cell CD4+ memory resting, while higher counts of B cell plasma and macrophage M2. CONCLUSION: Prognostic prediction model of 10 metabolic genes could accurately predict the prognosis of IHCC patients.