Exome sequencing of primary breast cancers with paired metastatic lesions reveals metastasis-enriched mutations in the A-kinase anchoring protein family (AKAPs)

TM-Score predicts immunotherapy efficacy and improves the performance of the machine learning prognostic model in gastric cancer

Immunotherapy is becoming increasingly important, but the overall response rate is relatively low in the treatment of gastric cancer (GC). The application of tumor mutational burden (TMB) in predicting immunotherapy efficacy in GC patients is limited and controversial, emphasizing the importance of optimizing TMB-based patient selection. By combining TMB and major histocompatibility complex (MHC) related hub genes, we established a novel TM-Score. This score showed superior performance for immunotherapeutic selection (AUC = 0.808) compared to TMB, MSI status, and EBV status. Additionally, it predicted the prognosis of GC patients. Subsequently, a machine learning model adjusted by the TM-Score further improved the accuracy of survival prediction (AUC > 0.8). Meanwhile, we found that GC patients with low TM-Score had a higher mutation frequency, higher expression of HLA genes and immune checkpoint genes, and higher infiltration of CD8+ T cells, CD4+ helper T cells, and M1 macrophages. This suggests that TM-Score is significantly associated with tumor immunogenicity and tumor immune environment. Notably, based on the RNA-seq and scRNA-seq, it was found that AKAP5, a key component gene of TM-Score, is involved in anti-tumor immunity by promoting the infiltration of CD4+ T cells, NK cells, and myeloid cells. Additionally, siAKAP5 significantly reduced MHC-II mRNA expression in the GC cell line. In addition, our immunohistochemistry assays confirmed a positive correlation between AKAP5 and human leukocyte antigen (HLA) expression. Furthermore, AKAP5 levels were higher in patients with longer survival and those who responded to immunotherapy in GC, indicating its potential value in predicting prognosis and immunotherapy outcomes. In conclusion, TM-Score, as an optimization of TMB, is a more precise biomarker for predicting the immunotherapy efficacy of the GC population. Additionally, AKAP5 shows promise as a therapeutic target for GC.

Divergent regulation of KCNQ1/E1 by targeted recruitment of protein kinase A to distinct sites on the channel complex

The slow delayed rectifier potassium current, IKs, conducted through pore-forming Q1 and auxiliary E1 ion channel complexes is important for human cardiac action potential repolarization. During exercise or fright, IKs is up-regulated by protein kinase A (PKA)-mediated Q1 phosphorylation to maintain heart rhythm and optimum cardiac performance. Sympathetic up-regulation of IKs requires recruitment of PKA holoenzyme (two regulatory - RI or RII - and two catalytic Cα subunits) to Q1 C-terminus by an A kinase anchoring protein (AKAP9). Mutations in Q1 or AKAP9 that abolish their functional interaction result in long QT syndrome type 1 and 11, respectively, which increases the risk of sudden cardiac death during exercise. Here, we investigated the utility of a targeted protein phosphorylation (TPP) approach to reconstitute PKA regulation of IKs in the absence of AKAP9. Targeted recruitment of endogenous Cα to E1-YFP using a GFP/YFP nanobody (nano) fused to RIIα enabled acute cAMP-mediated enhancement of IKs, reconstituting physiological regulation of the channel complex. By contrast, nano-mediated tethering of RIIα or Cα to Q1-YFP constitutively inhibited IKs by retaining the channel intracellularly in the endoplasmic reticulum and Golgi. Proteomic analysis revealed that distinct phosphorylation sites are modified by Cα targeted to Q1-YFP compared to free Cα. Thus, functional outcomes of synthetically recruited PKA on IKs regulation is critically dependent on the site of recruitment within the channel complex. The results reveal insights into divergent regulation of IKs by phosphorylation across different spatial and time scales, and suggest a TPP approach to develop new drugs to prevent exercise-induced sudden cardiac death.

Using the Proteomics Toolbox to Resolve Topology and Dynamics of Compartmentalized cAMP Signaling

cAMP is a second messenger that regulates a myriad of cellular functions in response to multiple extracellular stimuli. New developments in the field have provided exciting insights into how cAMP utilizes compartmentalization to ensure specificity when the message conveyed to the cell by an extracellular stimulus is translated into the appropriate functional outcome. cAMP compartmentalization relies on the formation of local signaling domains where the subset of cAMP signaling effectors, regulators and targets involved in a specific cellular response cluster together. These domains are dynamic in nature and underpin the exacting spatiotemporal regulation of cAMP signaling. In this review, we focus on how the proteomics toolbox can be utilized to identify the molecular components of these domains and to define the dynamic cellular cAMP signaling landscape. From a therapeutic perspective, compiling data on compartmentalized cAMP signaling in physiological and pathological conditions will help define the signaling events underlying disease and may reveal domain-specific targets for the development of precision medicine interventions.

Genetic variation as a long-distance modulator of RAD21 expression in humans

Somatic mutations and changes in expression of RAD21 are common in many types of cancer. Moreover, sub-optimal levels of RAD21 expression in early development can result in cohesinopathies. Altered RAD21 levels can result directly from mutations in the RAD21 gene. However, whether DNA variants outside of the RAD21 gene could control its expression and thereby contribute to cancer and developmental disease is unknown. In this study, we searched for genomic variants that modify RAD21expression to determine their potential to contribute to development or cancer by RAD21 dysregulation. We searched 42,953,834 genomic variants for a spatial-eQTL association with the transcription of RAD21. We identified 123 significant associations (FDR < 0.05), which are local (cis) or long-distance (trans) regulators of RAD21 expression. The 123 variants co-regulate a further seven genes (AARD, AKAP11, GRID1, KCNIP4, RCN1, TRIOBP, and USP32), enriched for having Sp2 transcription factor binding sites in their promoter regions. The Sp2 transcription factor and six of the seven genes had previously been associated with cancer onset, progression, and metastasis. Our results suggest that genome-wide variation in non-coding regions impacts on RAD21 transcript levels in addition to other genes, which then could impact on oncogenesis and the process of ubiquitination. This identification of distant co-regulation of oncogenes represents a strategy for discovery of novel genetic regions influencing cancer onset and a potential for diagnostics.

Targeted amplicon sequencing for primary tumors and matched lymph node metastases in patients with extrahepatic cholangiocarcinoma

BACKGROUND

Lymph node metastasis (LNM) is one of the most adverse prognostic factors in extrahepatic cholangiocarcinoma (EHCC) cases. As next-generation sequencing technology has become more widely available, the genomic profile of biliary tract carcinoma has been clarified. However, whether LNMs have additional genomic alterations in patients with EHCC has not been investigated. Here, we aimed to compare the genomic alterations between primary tumors and matched LNMs in patients with EHCC.

METHODS

Sixteen patients with node-positive EHCCs were included. Genomic DNA was extracted from tissue samples of primary tumors and matched LNMs. Targeted amplicon sequencing of 160 cancer-related genes was performed.

RESULTS

Among the 32 tumor samples from 16 patients, 91 genomic mutations were identified. Genomic mutations were noted in 31 genes, including TP53, MAP3K1, SMAD4, APC, and ARID1A. TP53 mutations were most frequently observed (12/32; 37.5%). Genomic mutation profiles were highly concordant between primary tumors and matched LNMs (13/16; 81.3%), and an additional genomic mutation of CDK12 was observed in only one patient.

CONCLUSION

Genomic mutations were highly concordant between primary tumors and matched LNMs, suggesting that genotyping of archived primary tumor samples may help predict genomic mutations of metastatic tumors in patients with EHCC.

Heterogeneity and tumor evolution reflected in liquid biopsy in metastatic breast cancer patients: a review

Breast cancer is a spatially and temporally dynamic disease in which differently evolving genetic clones are responsible for progression and clinical outcome. We review tumor heterogeneity and clonal evolution from studies comparing primary tumors and metastasis and discuss plasma circulating tumor DNA as a powerful real-time approach for monitoring the clonal landscape of breast cancer during treatment and recurrence. We found only a few early studies exploring clonal evolution and heterogeneity through analysis of multiregional tissue biopsies of different progression steps in comparison with circulating tumor DNA (ctDNA) from blood plasma. The model of linear progression seemed to be more often reported than the model of parallel progression. The results show complex routes to metastasis, however, and plasma most often reflected metastasis more than primary tumor. The described patterns of evolution and the polyclonal nature of breast cancer have clinical consequences and should be considered during patient diagnosis and treatment selection. Current studies focusing on the relevance of clonal evolution in the clinical setting illustrate the role of liquid biopsy as a noninvasive biomarker for monitoring clonal progression and response to treatment. In the clinical setting, circulating tumor DNA may be an ideal support for tumor biopsies to characterize the genetic landscape of the metastatic disease and to improve longitudinal monitoring of disease dynamics and treatment effectiveness through detection of residual tumor after resection, relapse, or metastasis within a particular patient.

Identification of key gene signatures for the overall survival of ovarian cancer

BACKGROUND

The five-year overall survival (OS) of advanced-stage ovarian cancer remains nearly 25-35%, although several treatment strategies have evolved to get better outcomes. A considerable amount of heterogeneity and complexity has been seen in ovarian cancer. This study aimed to establish gene signatures that can be used in better prognosis through risk prediction outcome for the survival of ovarian cancer patients. Different studies' heterogeneity into a single platform is presented to explore the penetrating genes for poor or better survival. The integrative analysis of multiple data sets was done to determine the genes that influence poor or better survival. A total of 6 independent data sets was considered. The Cox Proportional Hazard model was used to obtain significant genes that had an impact on ovarian cancer patients. The gene signatures were prepared by splitting the over-expressed and under-expressed genes parallelly by the variable selection technique. The data visualisation techniques were prepared to predict the overall survival, and it could support the therapeutic regime.

RESULTS

We preferred to select 20 genes in each data set as upregulated and downregulated. Irrespective of the selection of multiple genes, not even a single gene was found common among data sets for the survival of ovarian cancer patients. However, the same analytical approach adopted. The chord plot was presented to make a comprehensive understanding of the outcome.

CONCLUSIONS

This study helps us to understand the results obtained from different studies. It shows the impact of the heterogeneity from one study to another. It shows the requirement of integrated studies to make a holistic view of the gene signature for ovarian cancer survival.

Clinical, Pathological, and Molecular Features of Breast Carcinoma Cutaneous Metastasis

Cutaneous metastases (CMs) account for 2% of all skin malignancies, and nearly 70% of CMs in women originate from breast cancer (BC). CMs are usually associated with poor prognosis, are difficult to treat, and can pose diagnostic problems, such as in histopathological diagnosis when occurring long after development of the primary tumor. In addition, the molecular differences between the primary tumors and their CMs, and between CMs and metastases in other organs, are not well defined. Here, we review the main clinical, pathological, and molecular characteristics of breast cancer CMs. Identifying molecular markers in primary BC that predict CM and can be used to determine the molecular differences between primary tumors and their metastases is of great interest for the design of new therapeutic approaches.

Somatic Copy Number Alterations in Human Cancers: An Analysis of Publicly Available Data From The Cancer Genome Atlas

Somatic copy number alterations (SCNAs) are a pervasive trait of human cancers that contributes to tumorigenesis by affecting the dosage of multiple genes at the same time. In the past decade, The Cancer Genome Atlas (TCGA) and the International Cancer Genome Consortium (ICGC) initiatives have generated and made publicly available SCNA genomic profiles from thousands of tumor samples across multiple cancer types. Here, we present a comprehensive analysis of 853,218 SCNAs across 10,729 tumor samples belonging to 32 cancer types using TCGA data. We then discuss current models for how SCNAs likely arise during carcinogenesis and how genomic SCNA profiles can inform clinical practice. Lastly, we highlight open questions in the field of cancer-associated SCNAs.

Subcellular Organization of the cAMP Signaling Pathway

The field of cAMP signaling is witnessing exciting developments with the recognition that cAMP is compartmentalized and that spatial regulation of cAMP is critical for faithful signal coding. This realization has changed our understanding of cAMP signaling from a model in which cAMP connects a receptor at the plasma membrane to an intracellular effector in a linear pathway to a model in which cAMP signals propagate within a complex network of alternative branches and the specific functional outcome strictly depends on local regulation of cAMP levels and on selective activation of a limited number of branches within the network. In this review, we cover some of the early studies and summarize more recent evidence supporting the model of compartmentalized cAMP signaling, and we discuss how this knowledge is starting to provide original mechanistic insight into cell physiology and a novel framework for the identification of disease mechanisms that potentially opens new avenues for therapeutic interventions. SIGNIFICANCE STATEMENT: cAMP mediates the intracellular response to multiple hormones and neurotransmitters. Signal fidelity and accurate coordination of a plethora of different cellular functions is achieved via organization of multiprotein signalosomes and cAMP compartmentalization in subcellular nanodomains. Defining the organization and regulation of subcellular cAMP nanocompartments is necessary if we want to understand the complex functional ramifications of pharmacological treatments that target G protein-coupled receptors and for generating a blueprint that can be used to develop precision medicine interventions.

Oestrogen-regulated protein SLC39A6: a biomarker of good prognosis in luminal breast cancer

PURPOSE

The outcome of the luminal oestrogen receptor-positive (ER +) subtype of breast cancer (BC) is highly variable and patient stratification needs to be refined. We assessed the prognostic significance of oestrogen-regulated solute carrier family 39 member 6 (SLC39A6) in BC, with emphasis on ER + tumours.

MATERIALS AND METHODS

SLC39A6 mRNA expression and copy number alterations were assessed using the METABRIC cohort (n = 1980). SLC39A6 protein expression was evaluated in a large (n = 670) and annotated series of early-stage (I-III) operable BC using tissue microarrays and immunohistochemistry. The associations between SLC39A6 expression and clinicopathological parameters, patient outcomes and other ER-related markers were evaluated using Chi-square tests and Kaplan-Meier curves.

RESULTS

High SLC39A6 mRNA and protein expression was associated with features characteristic of less aggressive tumours in the entire BC cohort and ER + subgroup. SLC39A6 protein expression was detected in the cytoplasm and nuclei of the tumour cells. High SLC39A6 nuclear expression and mRNA levels were positively associated with ER + tumours and expression of ER-related markers, including the progesterone receptor, forkhead box protein A1 and GATA binding protein 3. In the ER + luminal BC, high SLC39A6 expression was independently associated with longer BC-specific survival (BCSS) (P = 0.015, HR 0.678, 95% CI 0.472‒0.972) even in those who did not receive endocrine therapy (P = 0.001, HR 0.701, 95% CI 0.463‒1.062).

CONCLUSION

SLC39A6 may be prognostic for a better outcome in ER + luminal BC. Further functional studies to investigate the role of SLC39A6 in ER + luminal BC are warranted.

Drugs That Regulate Local Cell Signaling: AKAP Targeting as a Therapeutic Option

Cells respond to environmental cues by mobilizing signal transduction cascades that engage protein kinases and phosphoprotein phosphatases. Correct organization of these enzymes in space and time enables the efficient and precise transmission of chemical signals. The cyclic AMP-dependent protein kinase A is compartmentalized through its association with A-kinase anchoring proteins (AKAPs). AKAPs are a family of multivalent scaffolds that constrain signaling enzymes and effectors at subcellular locations to drive essential physiological events. More recently, it has been recognized that defective signaling in certain endocrine disorders and cancers proceeds through pathological AKAP complexes. Consequently, pharmacologically targeting these macromolecular complexes unlocks new therapeutic opportunities for a growing number of clinical indications. This review highlights recent findings on AKAP signaling in disease, particularly in certain cancers, and offers an overview of peptides and small molecules that locally regulate AKAP-binding partners.

Low expression of A-kinase anchor protein 5 predicts poor prognosis in non-mucin producing stomach adenocarcinoma based on TCGA data

BACKGROUND

In the past, there were not a lot of studies on how A-kinase anchor protein 5 (AKAP5) involving in the pathogenesis and prognosis of non-mucin producing stomach adenocarcinoma (NMSA). Therefore, we studied the relationship between AKAP5 and the prognosis of NMSA and its possible mechanisms using publicly available data from The Cancer Genome Atlas (TCGA).

METHODS

RNA high-throughput sequencing and clinicopathologic data of NMSA were downloaded from the TCGA. Clinical pathologic features associated with AKAP5 expression were analyzed using the chi-square and Fisher exact tests. The relationship between the overall survival (OS) and AKAP5 expression was analyzed by the Kaplan-Meier method and the Cox regression analysis. GSEA analysis was performed using the TCGA dataset.

RESULTS

Our results indicated that the AKAP5 expression was increased in NMSA (all tumor vs. adjacent mucosa). Also, histologic grade, clinical stage, N classification, and survival status were significantly correlated with AKAP5 expression. Kaplan-Meier curves showed that low AKAP5 expression was associated with a poor OS among the NMSA patients (P=5.003e-05), and in the clinical stage III and IV (P=4.646e-05), TNM stage T3 (P=0.016), T4 (P=0.001), N2 (P=0.012), N3 (P=0.003), M0 (P=3.911e-05), and histological grade G3 (P=1.658e-04) subgroups. Cox regression analysis showed that reduced AKAP5 expression in NMSA is associated with age (HR =1.03, P=0.007), stage (HR =1.84 for stage I, II vs. stage III, IV, P=0.002) and M classification (HR =1.8 for M0 vs. M1, P=0.010). Gene sets related to cholesterol homeostasis, glycolysis, estrogen response late, adipogenesis, estrogen response early, notch signaling, and peroxisome were differentially enriched with the low AKAP5 expression phenotype.

CONCLUSIONS

Low expression of AKAP5 may be a potential molecular marker for predicting poor prognosis of NMSA. Besides, cholesterol homeostasis, glycolysis, estrogen response, adipogenesis, notch signaling, and peroxisome may be the key pathways regulated by AKAP5 in NMSA. It also suggested that AKAP5 might potentially have biological functions in the development of stomach adenocarcinoma.

Gene signatures and prognostic analyses of the Tob/BTG pituitary tumor-transforming gene (PTTG) family in clinical breast cancer patients

Breast cancer is the most common cancer type in females, and exploring the mechanisms of disease progression is playing a crucial role in the development of potential therapeutics. Pituitary tumor-transforming gene (PTTG) family members are well documented to be involved in cell-cycle regulation and mitosis, and contribute to cancer development by their involvement in cellular transformation in several tumor types. The critical roles of PTTG family members as crucial transcription factors in diverse types of cancers are recognized, but how they regulate breast cancer development still remains mostly unknown. Meanwhile, a holistic genetic analysis exploring whether PTTG family members regulate breast cancer progression via the cell cycle as well as the energy metabolism-related network is lacking. To comprehensively understand the messenger RNA expression profiles of PTTG proteins in breast cancer, we herein conducted a high-throughput screening approach by integrating information from various databases such as Oncomine, Kaplan-Meier Plotter, Metacore, ClueGo, and CluePedia. These useful databases and tools provide expression profiles and functional analyses. The present findings revealed that PTTG1 and PTTG3 are two important genes with high expressions in breast cancer relative to normal breast cells, implying their unique roles in breast cancer progression. Results of our coexpression analysis demonstrated that PTTG family genes were positively correlated with thiamine triphosphate (TTP), deoxycytidine triphosphate (dCTP) metabolic, glycolysis, gluconeogenesis, and cell-cycle related pathways. Meanwhile, through Cytoscape analyzed indicated that in addition to the metastasis markers AURKA, AURKB, and NDC80, many of the kinesin superfamily (KIF) members including KIFC1, KIF2C, KIF4A, KIF14, KIF20A, KIF23, were also correlated with PTTG family transcript expression. Finally, we revealed that high levels of PTTG1 and PTTG3 transcription predicted poor survival, which provided useful insights into prospective research of cancer associated with the PTTG family. Therefore, these members of the PTTG family would serve as distinct and essential prognostic biomarkers in breast cancer.

Single-cell RNA-seq variant analysis for exploration of genetic heterogeneity in cancer

Inter- and intra-tumour heterogeneity is caused by genetic and non-genetic factors, leading to severe clinical implications. High-throughput sequencing technologies provide unprecedented tools to analyse DNA and RNA in single cells and explore both genetic heterogeneity and phenotypic variation between cells in tissues and tumours. Simultaneous analysis of both DNA and RNA in the same cell is, however, still in its infancy. We have thus developed a method to extract and analyse information regarding genetic heterogeneity that affects cellular biology from single-cell RNA-seq data. The method enables both comparisons and clustering of cells based on genetic variation in single nucleotide variants, revealing cellular subpopulations corroborated by gene expression-based methods. Furthermore, the results show that lymph node metastases have lower levels of genetic heterogeneity compared to their original tumours with respect to variants affecting protein function. The analysis also revealed three previously unknown variants common across cancer cells in glioblastoma patients. These results demonstrate the power and versatility of scRNA-seq variant analysis and highlight it as a useful complement to already existing methods, enabling simultaneous investigations of both gene expression and genetic variation.

Depletion of dAKAP1-protein kinase A signaling islands from the outer mitochondrial membrane alters breast cancer cell metabolism and motility

Breast cancer screening and new precision therapies have led to improved patient outcomes. Yet, a positive prognosis is less certain when primary tumors metastasize. Metastasis requires a coordinated program of cellular changes that promote increased survival, migration, and energy consumption. These pathways converge on mitochondrial function, where distinct signaling networks of kinases, phosphatases, and metabolic enzymes regulate these processes. The protein kinase A-anchoring protein dAKAP1 compartmentalizes protein kinase A (PKA) and other signaling enzymes at the outer mitochondrial membrane and thereby controls mitochondrial function and dynamics. Modulation of these processes occurs in part through regulation of dynamin-related protein 1 (Drp1). Here, we report an inverse relationship between the expression of dAKAP1 and mesenchymal markers in breast cancer. Molecular, cellular, and in silico analyses of breast cancer cell lines confirmed that dAKAP1 depletion is associated with impaired mitochondrial function and dynamics, as well as with increased glycolytic potential and invasiveness. Furthermore, disruption of dAKAP1-PKA complexes affected cell motility and mitochondrial movement toward the leading edge in invasive breast cancer cells. We therefore propose that depletion of dAKAP1-PKA "signaling islands" from the outer mitochondrial membrane augments progression toward metastatic breast cancer.

Analysis of mutations in primary and metastatic synovial sarcoma

Synovial sarcoma is the most common pediatric non-rhabdomyosarcoma soft tissue sarcoma and accounts for about 8-10% of all soft tissue sarcoma in childhood and adolescence. The presence of a chromosomal translocation-associated SS18-SSX-fusion gene is causally linked to development of primary synovial sarcoma. Metastases occur in approximately 50-70% of synovial sarcoma cases with yet unknown mechanisms, which led to about 70-80% mortality rate in five years. To explore the possibilities to investigate metastatic mechanisms of synovial sarcoma, we carried out the first genome-wide search for potential genetic biomarkers and drivers associated with metastasis by comparative mutational profiling of 18 synovial sarcoma samples isolated from four patients carrying the primary tumors and another four patients carrying the metastatic tumors through whole exome sequencing. Selected from the candidates yielded from this effort, we examined the effect of the multiple missense mutations of ADAM17, which were identified solely in metastatic synovial sarcoma. The mutant alleles as well as the wild-type control were expressed in the mammalian cells harboring the SS18-SSX1 fusion gene. The ADAM17-P729H mutation was shown to enhance cell migration, a phenotype associated with metastasis. Therefore, like ADAM17-P729H, other mutations we identified solely in metastatic synovial sarcoma may also have the potential to serve as an entry point for unraveling the metastatic mechanisms of synovial sarcoma.

Genetic, Epigenetic, and Immunologic Profiling of MMR-Deficient Relapsed Glioblastoma

PURPOSE

In-depth characterization of recurrent glioblastoma (rGBM) might contribute to a better understanding of the mechanisms behind tumor progression and enable rGBM treatment with targeted drugs.Experimental Design: In this study, GBM samples were collected at diagnosis and recurrence from adult patients treated with Stupp protocol. Expression of mismatch repair (MMR) proteins was evaluated by IHC, followed by whole exome sequencing (WES) of tumor samples showing loss of MSH6 reactivity. Established genetic, epigenetic, and immunologic markers were assessed by standard methods and correlated with loss of MMR proteins and patient survival.

RESULTS

Expression of MMR proteins was partially or completely lost in 25.9% rGBM samples. Specifically, 12 samples showed partial or total MSH6 expression reduction. Conversely, 96.4% of GBM samples at diagnosis expressed MMR markers. WES disclosed lack of variants in MMR genes in primary samples, whereas two MSH6-negative rGBM samples shared a c.3438+1G>A* splicing MSH6 variant with a potential loss of function effect. MSH6-negative rGBM specimens had high tumor mutational burden (TMB), but no microsatellite instability. In contrast, GBM samples with partial loss of MMR proteins disclosed low TMB. MMR-deficient rGBM showed significant telomere shortening and MGMT methylation and are characterized by highly heterogeneous MHC class I expression.

CONCLUSIONS

Multilevel profiling of MMR-deficient rGBM uncovered hypermutated genotype uncoupled from enriched expression of immune-related markers. Assessment of MHC class I expression and TMB should be included in protocols aiming to identify rGBM patients potentially eligible for treatment with drugs targeting immune-checkpoint inhibitors.