USP44 overexpression drives a MYC-like gene expression program in neuroblastoma through epigenetic reprogramming.

Neuroblastoma is an embryonic cancer that contributes disproportionately to death in young children. Sequencing data have uncovered few recurrently mutated genes in this cancer, although epigenetic pathways have been implicated in disease pathogenesis. We used an expression-based computational screen that examined the impact of deubiquitinating enzymes on patient survival to identify potential new targets. We identified the histone H2B deubiquitinating enzyme USP44 as the enzyme with the greatest impact on survival in patients with neuroblastoma. High levels of USP44 significantly correlate with metastatic disease, unfavorable histology, advanced patient age, and MYCN-amplification. The subset of patients with tumors expressing high levels of USP44 had a significantly worse survival, including those with tumors lacking MYCN amplification. We showed experimentally that USP44 regulates neuroblastoma cell proliferation, migration, invasion, and neuronal development. Depletion of the histone H2B ubiquitin ligase subunit RNF20 resulted in similar findings, strongly implicating this histone mark as the target of USP44 activity in this disease. Integration of transcriptome and epigenome in analyses demonstrates a distinct set of genes that is regulated by USP44, including those in Hallmark MYC target genes in both murine embryonic fibroblasts and the SH-SY5Y neuroblastoma cell line. We conclude that USP44 is a novel epigenetic regulator that promotes aggressive features and may be a novel target in neuroblastoma. Implications: This study identifies a new genetic marker of aggressive neuroblastoma and identifies the mechanisms by which its overactivity contributes to pathophysiology in this disease.

MDM2 amplification in rod-shaped chromosomes provides clues to early stages of circularized gene amplification in liposarcoma.

Well-differentiated liposarcoma (WDLS) displays amplification of genes on chromosome 12 (Chr12) in supernumerary ring or giant marker chromosomes. These structures have been suggested to develop through chromothripsis, followed by circularization and breakage-fusion-bridge (BFB) cycles. To test this hypothesis, we compared WDLSs with Chr12 amplification in rod-shaped chromosomes with WDLSs with rings. Both types of amplicons share the same spectrum of structural variants (SVs), show higher SV frequencies in Chr12 than in co-amplified segments, have SVs that fuse the telomeric ends of co-amplified chromosomes, and lack interspersed deletions. Combined with the finding of cells with transient rod-shaped structures in tumors with ring chromosomes, this suggests a stepwise process starting with the gain of Chr12 material that, after remodeling which does not fit with classical chromothripsis, forms a dicentric structure with other chromosomes. Depending on if and when telomeres from other chromosomes are captured, circularized or linear gain of 12q sequences will predominate.

A consensus molecular subtypes classification strategy for clinical colorectal cancer tissues.

Consensus Molecular Subtype (CMS) classification of colorectal cancer (CRC) tissues is complicated by RNA degradation upon formalin-fixed paraffin-embedded (FFPE) preservation. Here, we present an FFPE-curated CMS classifier. The CMSFFPE classifier was developed using genes with a high transcript integrity in FFPE-derived RNA. We evaluated the classification accuracy in two FFPE-RNA datasets with matched fresh-frozen (FF) RNA data, and an FF-derived RNA set. An FFPE-RNA application cohort of metastatic CRC patients was established, partly treated with anti-EGFR therapy. Key characteristics per CMS were assessed. Cross-referenced with matched benchmark FF CMS calls, the CMSFFPE classifier strongly improved classification accuracy in two FFPE datasets compared with the original CMSClassifier (63.6% versus 40.9% and 83.3% versus 66.7%, respectively). We recovered CMS-specific recurrence-free survival patterns (CMS4 versus CMS2: hazard ratio 1.75, 95% CI 1.24-2.46). Key molecular and clinical associations of the CMSs were confirmed. In particular, we demonstrated the predictive value of CMS2 and CMS3 for anti-EGFR therapy response (CMS2&3: odds ratio 5.48, 95% CI 1.10-27.27). The CMSFFPE classifier is an optimized FFPE-curated research tool for CMS classification of clinical CRC samples.

The pluripotency factor NANOG contributes to mesenchymal plasticity and is predictive for outcome in esophageal adenocarcinoma.

BACKGROUND: Despite the advent of neoadjuvant chemoradiotherapy (CRT), overall survival rates of esophageal adenocarcinoma (EAC) remain low. A readily induced mesenchymal transition of EAC cells contributes to resistance to CRT. METHODS: In this study, we aimed to chart the heterogeneity in cell state transition after CRT and to identify its underpinnings. A panel of 12 esophageal cultures were treated with CRT and ranked by their relative epithelial-mesenchymal plasticity. RNA-sequencing was performed on 100 pre-treatment biopsies. After RNA-sequencing, Ridge regression analysis was applied to correlate gene expression to ranked plasticity, and models were developed to predict mesenchymal transitions in patients. Plasticity score predictions of the three highest significant predictive models were projected on the pre-treatment biopsies and related to clinical outcome data. Motif enrichment analysis of the genes associated with all three models was performed. RESULTS: This study reveals NANOG as the key associated transcription factor predicting mesenchymal plasticity in EAC. Expression of NANOG in pre-treatment biopsies is highly associated with poor response to neoadjuvant chemoradiation, the occurrence of recurrences, and median overall survival difference in EAC patients (>48 months). Perturbation of NANOG reduces plasticity and resensitizes cell lines, organoid cultures, and patient-derived in vivo grafts. CONCLUSIONS: In conclusion, NANOG is a key transcription factor in mesenchymal plasticity in EAC and a promising predictive marker for outcome.

Esophageal cancer is the sixth most common cause of cancer-related death worldwide. Although chemotherapy combined with radiotherapy (chemoradiotherapy) followed by surgery has improved survival, tumor recurrence and metastatic disease (that has spread to other parts of the body) are often observed after several months. In this study, we assessed the effect of chemoradiotherapy on esophageal cells in the lab to predict the effect in patients with esophageal cancer. To investigate this, genes were assessed from 12 different cell lines and 100 patient tissues. We revealed that levels of one of the genes, NANOG, associates with poor response in patients. NANOG could be a promising marker to predict outcome in patients with esophageal cancer. This knowledge might help clinicians to treat patients with esophageal cancer appropriately, or may lead to new or optimized treatments.

The neo-open reading frame peptides that comprise the tumor framome are a rich source of neoantigens for cancer immunotherapy.

Identification of immunogenic cancer neoantigens as targets for therapy is challenging. Here, we integrate cancer whole genome and long-read transcript sequencing to identify the collection of novel open reading frame peptides (NOPs) expressed in tumors, termed the framome. NOPs represent tumor-specific peptides that are different from wild-type proteins and may be strongly immunogenic. We describe an uncharacterized class of hidden NOPs, which derive from structural genomic variants involving an upstream protein coding gene driving expression and translation of non-coding regions of the genome downstream of a rearrangement breakpoint. NOPs represent a vast amount of possible neoantigens particularly in tumors with many (complex) structural genomic variants and a low number of missense mutations. We show that NOPs are immunogenic and epitopes derived from NOPs can bind to MHC class I molecules. Finally, we provide evidence for the presence of memory T-cells specific for hidden NOPs in lung cancer patient peripheral blood.

Integrative analysis of neuroblastoma by single-cell RNA sequencing identifies the NECTIN2-TIGIT axis as a target for immunotherapy.

Pediatric patients with high-risk neuroblastoma have poor survival rates and urgently need more effective treatment options with less side effects. Since novel and improved immunotherapies may fill this need, we dissect the immunoregulatory interactions in neuroblastoma by single-cell RNA-sequencing of 24 tumors (10 pre- and 14 post-chemotherapy, including 5 pairs) to identify strategies for optimizing immunotherapy efficacy. Neuroblastomas are infiltrated by natural killer (NK), T and B cells, and immunosuppressive myeloid populations. NK cells show reduced cytotoxicity and T cells have a dysfunctional profile. Interaction analysis reveals a vast immunoregulatory network and identifies NECTIN2-TIGIT as a crucial immune checkpoint. Combined blockade of TIGIT and PD-L1 significantly reduces neuroblastoma growth, with complete responses (CR) in vivo. Moreover, addition of TIGIT+PD-L1 blockade to standard relapse treatment in a chemotherapy-resistant Th-ALKF1174L/MYCN 129/SvJ syngeneic model induces CR. In conclusion, our integrative analysis provides promising targets and a rationale for immunotherapeutic combination strategies.

Germline specific genes increase DNA double-strand break repair and radioresistance in lung adenocarcinoma cells.

In principle, germline cells possess the capability to transmit a nearly unaltered set of genetic material to infinite future generations, whereas somatic cells are limited by strict growth constraints necessary to assure an organism's physical structure and eventual mortality. As the potential to replicate indefinitely is a key feature of cancer, we hypothesized that the activation of a "germline program" in somatic cells can contribute to oncogenesis. Our group recently described over one thousand germline specific genes that can be ectopically expressed in cancer, yet how germline specific processes contribute to the malignant properties of cancer is poorly understood. We here show that the expression of germ cell/cancer (GC) genes correlates with malignancy in lung adenocarcinoma (LUAD). We found that LUAD cells expressing more GC genes can repair DNA double strand breaks more rapidly, show higher rates of proliferation and are more resistant to ionizing radiation, compared to LUAD cells that express fewer GC genes. In particular, we identified the HORMA domain protein regulator TRIP13 to be predominantly responsible for this malignant phenotype, and that TRIP13 inhibition or expression levels affect the response to ionizing radiation and subsequent DNA repair. Our results demonstrate that GC genes are viable targets in oncology, as they induce increased radiation resistance and increased propagation in cancer cells. Because their expression is normally restricted to germline cells, we anticipate that GC gene directed therapeutic options will effectively target cancer, with limited side effects besides (temporary) infertility.

Consensus molecular subtype transition during progression of colorectal cancer.

The consensus molecular subtype (CMS) classification divides colorectal cancer (CRC) into four distinct subtypes based on RNA expression profiles. The biological differences between CMSs are already present in CRC precursor lesions, but not all CMSs pose the same risk of malignant transformation. To fully understand the path to malignant transformation and to determine whether CMS is a fixed entity during progression, genomic and transcriptomic data from two regions of the same CRC lesion were compared: the precursor region and the carcinoma region. In total, 24 patients who underwent endoscopic removal of T1-2 CRC were included. Regions were subtyped for CMS and DNA mutation analysis was performed. Additionally, a set of 85 benign adenomas was CMS-subtyped. This analysis revealed that almost all benign adenomas were classified as CMS3 (91.8%). In contrast, CMS2 was the most prevalent subtype in precursor regions (66.7%), followed by CMS3 (29.2%). CMS4 was absent in precursor lesions and originated at the carcinoma stage. Importantly, CMS switching occurred in a substantial number of cases and almost all (six out of seven) CMS3 precursor regions showed a shift to a different subtype in the carcinoma part of the lesion, which in four cases was classified as CMS4. In conclusion, our data indicate that CMS3 is related to a more indolent type of precursor lesion that less likely progresses to CRC and when this occurs, it is often associated with a subtype change that includes the more aggressive mesenchymal CMS4. In contrast, an acquired CMS2 signature appeared to be rather fixed during early CRC development. Combined, our data show that subtype changes occur during progression and that CMS3 switching is related to changes in the genomic background through acquisition of a novel driver mutation (TP53) or selective expansion of a clone, but also occurred independently of such genetic changes. © 2023 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Two opposing gene expression patterns within ATRX aberrant neuroblastoma.

Neuroblastoma is the most common extracranial solid tumor in children. A subgroup of high-risk patients is characterized by aberrations in the chromatin remodeller ATRX that is encoded by 35 exons. In contrast to other pediatric cancer where ATRX point mutations are most frequent, multi-exon deletions (MEDs) are the most frequent type of ATRX aberrations in neuroblastoma. 75% of these MEDs are predicted to produce in-frame fusion proteins, suggesting a potential gain-of-function effect compared to nonsense mutations. For neuroblastoma there are only a few patient-derived ATRX aberrant models. Therefore, we created isogenic ATRX aberrant models using CRISPR-Cas9 in several neuroblastoma cell lines and one tumoroid and performed total RNA-sequencing on these and the patient-derived models. Gene set enrichment analysis (GSEA) showed decreased expression of genes related to both ribosome biogenesis and several metabolic processes in our isogenic ATRX exon 2-10 MED model systems, the patient-derived MED models and in tumor data containing two patients with an ATRX exon 2-10 MED. In sharp contrast, these same processes showed an increased expression in our isogenic ATRX knock-out and exon 2-13 MED models. Our validations confirmed a role of ATRX in the regulation of ribosome homeostasis. The two distinct molecular expression patterns within ATRX aberrant neuroblastomas that we identified imply that there might be a need for distinct treatment regimens.

Iron-responsive element of Divalent metal transporter 1 (Dmt1) controls Notch-mediated cell fates.

Notch receptor activation is regulated by the intramembrane protease γ-secretase, which cleaves and liberates the Notch intracellular domain (Nicd) that regulates gene transcription. While γ-secretase cleavage is necessary, we demonstrate it is insufficient for Notch activation and requires vesicular trafficking. Here, we report Divalent metal transporter 1 (Dmt1, Slc11A2) as a novel and essential regulator of Notch signalling. Dmt1-deficient cells are defective in Notch signalling and have perturbed endolysosomal trafficking and function. Dmt1 encodes for two isoforms, with and without an iron response element (ire). We show that isoform-specific silencing of Dmt1-ire and Dmt1+ire has opposite consequences on Notch-dependent cell fates in cell lines and intestinal organoids. Loss of Dmt1-ire suppresses Notch activation and promotes differentiation, whereas loss of Dmt1+ire causes Notch activation and maintains stem-progenitor cell fates. Dmt1 isoform expression correlates with Notch and Wnt signalling in Apc-deficient intestinal organoids and human colorectal cancers. Consistently, Dmt1-ire silencing induces Notch-dependent differentiation in colorectal cancer cells. These data identify Dmt1 isoforms as binary switches controlling Notch cell fate decisions in normal and tumour cells.

Transcriptomic Differences Underlying the Activin-A Induced Large Osteoclast Formation in Both Healthy Control and Fibrodysplasia Ossificans Progressiva Osteoclasts.

Fibrodysplasia Ossificans Progressiva (FOP) is a very rare genetic disease characterized by progressive heterotopic ossification (HO) of soft tissues, leading to immobility and premature death. FOP is caused by a mutation in the Activin receptor Type 1 (ACVR1) gene, resulting in altered responsiveness to Activin-A. We recently revealed that Activin-A induces fewer, but larger and more active, osteoclasts regardless of the presence of the mutated ACVR1 receptor. The underlying mechanism of Activin-A-induced changes in osteoclastogenesis at the gene expression level remains unknown. Transcriptomic changes induced by Activin-A during osteoclast formation from healthy controls and patient-derived CD14-positive monocytes were studied using RNA sequencing. CD14-positive monocytes from six FOP patients and six age- and sex-matched healthy controls were differentiated into osteoclasts in the absence or presence of Activin-A. RNA samples were isolated after 14 days of culturing and analyzed by RNA sequencing. Non-supervised principal component analysis (PCA) showed that samples from the same culture conditions (e.g., without or with Activin-A) tended to cluster, indicating that the variability induced by Activin-A treatment was larger than the variability between the control and FOP samples. RNA sequencing analysis revealed 1480 differentially expressed genes induced by Activin-A in healthy control and FOP osteoclasts with p(adj) < 0.01 and a Log2 fold change of ≥±2. Pathway and gene ontology enrichment analysis revealed several significantly enriched pathways for genes upregulated by Activin-A that could be linked to the differentiation or function of osteoclasts, cell fusion or inflammation. Our data showed that Activin-A has a substantial effect on gene expression during osteoclast formation and that this effect occurred regardless of the presence of the mutated ACVR1 receptor causing FOP.

Vaccination provides superior in vivo recall capacity of SARS-CoV-2-specific memory CD8 T cells.

Memory CD8 T cells play an important role in the protection against breakthrough infections with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Whether the route of antigen exposure impacts these cells at a functional level is incompletely characterized. Here, we compare the memory CD8 T cell response against a common SARS-CoV-2 epitope after vaccination, infection, or both. CD8 T cells demonstrate comparable functional capacity when restimulated directly ex vivo, independent of the antigenic history. However, analysis of T cell receptor usage shows that vaccination results in a narrower scope than infection alone or in combination with vaccination. Importantly, in an in vivo recall model, memory CD8 T cells from infected individuals show equal proliferation but secrete less tumor necrosis factor (TNF) compared with those from vaccinated people. This difference is negated when infected individuals have also been vaccinated. Our findings shed more light on the differences in susceptibility to re-infection after different routes of SARS-CoV-2 antigen exposure.

Transcriptome analysis stratifies second-generation non-WNT/non-SHH medulloblastoma subgroups into clinically tractable subtypes.

Medulloblastoma (MB), one of the most common malignant pediatric brain tumor, is a heterogenous disease comprised of four distinct molecular groups (WNT, SHH, Group 3, Group 4). Each of these groups can be further subdivided into second-generation MB (SGS MB) molecular subgroups, each with distinct genetic and clinical characteristics. For instance, non-WNT/non-SHH MB (Group 3/4) can be subdivided molecularly into eight distinct and clinically relevant tumor subgroups. A further molecular stratification/summarization of these SGS MB would allow for the assignment of patients to risk-associated treatment protocols. Here, we performed DNA- and RNA-based analysis of 574 non-WNT/non-SHH MB and analyzed the clinical significance of various molecular patterns within the entire cohort and the eight SGS MB, with the aim to develop an optimal risk stratification of these tumors. Multigene analysis disclosed several survival-associated genes highly specific for each molecular subgroup within this non-WNT/non-SHH MB cohort with minimal inter-subgroup overlap. These subgroup-specific and prognostically relevant genes were associated with pathways that could underlie SGS MB clinical-molecular diversity and tumor-driving mechanisms. By combining survival-associated genes within each SGS MB, distinct metagene sets being appropriate for their optimal risk stratification were identified. Defined subgroup-specific metagene sets were independent variables in the multivariate models generated for each SGS MB and their prognostic value was confirmed in a completely non-overlapping validation cohort of non-WNT/non-SHH MB (n = 377). In summary, the current results indicate that the integration of transcriptome data in risk stratification models may improve outcome prediction for each non-WNT/non-SHH SGS MB. Identified subgroup-specific gene expression signatures could be relevant for clinical implementation and survival-associated metagene sets could be adopted for further SGS MB risk stratification. Future studies should aim at validating the prognostic role of these transcriptome-based SGS MB subtypes in prospective clinical trials.

Hsa-miR-323a-3p functions as a tumor suppressor and targets STAT3 in neuroblastoma cells.

Background: Studies conducted in the last decades have revealed a role for the non-coding microRNAs (miRNAs) in cancer development and progression. Several miRNAs within the chromosome region 14q32, a region commonly deleted in cancers, are associated with poor clinical outcome in the childhood cancer neuroblastoma. We have previously identified miR-323a-3p from this region to be downregulated in chemotherapy treated neuroblastoma cells compared to pre-treatment cells from the same patients. Furthermore, in neuroblastoma tumors, this miRNA is downregulated in advanced stage 4 disease compared to stage 1-2. In this study, we attempt to delineate the unknown functional roles of miR-323a-3p in neuroblastoma. Methods: Synthetic miRNA mimics were used to overexpress miR-323a-3p in neuroblastoma cell lines. To investigate the functional roles of miR-323a-3p, cell viability assay, flow cytometry, reverse transcription-quantitative polymerase chain reaction, luciferase reporter assay and western blot were conducted on the neuroblastoma cell lines Kelly, SH-SY5Y and SK-N-BE(2)-C. Results: Ectopic expression of miR-323a-3p resulted in marked reduction of cell viability in Kelly, SH-SY5Y and SK-N-BE(2)-C by causing G1-cell cycle arrest in Kelly and SH-SY5Y and apoptosis in all the cell lines tested. Furthermore, mRNA and protein levels of signal transducer and activator of transcription 3 (STAT3) were reduced upon miR-323a-3p overexpression. A direct binding of the miR-323a-3p to the 3'UTR of STAT3 was experimentally validated by luciferase reporter assay, where miR-323a-3p reduced luminescent signal from full length STAT3 3'UTR luciferase reporter, but not from a reporter with mutation in the predicted seed sequence. Conclusions: miR-323a-3p inhibits growth of neuroblastoma cell lines through G1-cell cycle arrest and apoptosis, and the well-known oncogene STAT3 is a direct target of this miRNA.

MEK inhibition causes BIM stabilization and increased sensitivity to BCL-2 family member inhibitors in RAS-MAPK-mutated neuroblastoma.

Introduction: Mutations affecting the RAS-MAPK pathway occur frequently in relapsed neuroblastoma tumors and are associated with response to MEK inhibition in vitro. However, these inhibitors alone do not lead to tumor regression in vivo, indicating the need for combination therapy. Methods and results: Via high-throughput combination screening, we identified that the MEK inhibitor trametinib can be combined with BCL-2 family member inhibitors, to efficiently inhibit growth of neuroblastoma cell lines with RAS-MAPK mutations. Suppressing the RAS-MAPK pathway with trametinib led to an increase in pro-apoptotic BIM, resulting in more BIM binding to anti-apoptotic BCL-2 family members. By favoring the formation of these complexes, trametinib treatment enhances sensitivity to compounds targeting anti-apoptotic BCL-2 family members. In vitro validation studies confirmed that this sensitizing effect is dependent on an active RAS-MAPK pathway. In vivo combination of trametinib with BCL-2 inhibitors led to tumor inhibition in NRAS-mutant and NF1-deleted xenografts. Conclusion: Together, these results show that combining MEK inhibition with BCL-2 family member inhibition could potentially improve therapeutic outcomes for RAS-MAPK-mutated neuroblastoma patients.

SOX11 regulates SWI/SNF complex components as member of the adrenergic neuroblastoma core regulatory circuitry.

The pediatric extra-cranial tumor neuroblastoma displays a low mutational burden while recurrent copy number alterations are present in most high-risk cases. Here, we identify SOX11 as a dependency transcription factor in adrenergic neuroblastoma based on recurrent chromosome 2p focal gains and amplifications, specific expression in the normal sympatho-adrenal lineage and adrenergic neuroblastoma, regulation by multiple adrenergic specific (super-)enhancers and strong dependency on high SOX11 expression in adrenergic neuroblastomas. SOX11 regulated direct targets include genes implicated in epigenetic control, cytoskeleton and neurodevelopment. Most notably, SOX11 controls chromatin regulatory complexes, including 10 SWI/SNF core components among which SMARCC1, SMARCA4/BRG1 and ARID1A. Additionally, the histone deacetylase HDAC2, PRC1 complex component CBX2, chromatin-modifying enzyme KDM1A/LSD1 and pioneer factor c-MYB are regulated by SOX11. Finally, SOX11 is identified as a core transcription factor of the core regulatory circuitry (CRC) in adrenergic high-risk neuroblastoma with a potential role as epigenetic master regulator upstream of the CRC.

Anti-GD2 immunoliposomes loaded with oxamate for neuroblastoma.

Oncometabolism can be targeted for the development of less myelotoxic oncotherapeutics. Lactate dehydrogenase A (LDHA) is central to the Warburg effect, a potential oncometabolic shift in neuroblastoma (NBL). Advanced surgical, cytotoxic and cell-differentiating therapies improved survival of children with NBL. Anti-GD2 monoclonal antibodies (mAb) effectively targeting NBL are also incorporated into complex therapies. However, poor clinical outcomes of high-risk NBL require improvements. Here, we verified the pre-reported prognostic value of LDHA expression in NBL using the R2 onco-genomics platform. Kaplan-Meier curves re-demonstrated that higher tumor LDHA expression correlates with worse survival. Multivariate statistics confirmed LDHA is independent from age, stage, and MYCN amplification. In conclusion, a molecular construct is proposed with anti-GD2 mAbs utilized for the targeted delivery of liposomes containing an LDHA inhibitor, Oxamate. Development and preclinical testing of this immunoliposome may validate targeted inhibition of the Warburg effect for NBL. IMPACT: Development of therapeutics against oncometabolism. Targeted specified drug-delivery with mAb. Sparing normal tissues from profound LDHA inhibition. Immunoliposome loaded with an anti-metabolite. If preclinically successful, has translational potential.

Comparison of transcriptome profiles between medulloblastoma primary and recurrent tumors uncovers novel variance effects in relapses.

Nowadays medulloblastoma (MB) tumors can be treated with risk-stratified approaches with up to 80% success rate. However, disease relapses occur in approximately 30% of patients and successful salvage treatment strategies at relapse remain scarce. Acquired copy number changes or TP53 mutations are known to occur frequently in relapses, while methylation profiles usually remain highly similar to those of the matching primary tumors, indicating that in general molecular subgrouping does not change during the course of the disease. In the current study, we have used RNA sequencing data to analyze the transcriptome profiles of 43 primary-relapse MB pairs in order to identify specific molecular features of relapses within various tumor groups. Gene variance analysis between primary and relapse samples demonstrated the impact of age in SHH-MB: the changes in gene expression relapse profiles were more pronounced in the younger patients (< 10 years old), which were also associated with increased DNA aberrations and somatic mutations at relapse probably driving this effect. For Group 3/4 MB transcriptome data analysis uncovered clear sets of genes either active or decreased at relapse that are significantly associated with survival, thus could be potential predictive markers. In addition, deconvolution analysis of bulk transcriptome data identified progression-associated differences in cell type enrichment. The proportion of undifferentiated progenitors increased in SHH-MB relapses with a concomitant decrease of differentiated neuron-like cells, while in Group 3/4 MB relapses cell cycle activity increases and differentiated neuron-like cells proportion decreases as well. Thus, our findings uncovered significant transcriptome changes in the molecular signatures of relapsed MB and could be potentially useful for further clinical purposes.