Sequential analysis of cfDNA reveals clonal evolution in neuroblastoma patients receiving ALK targeted therapy.

PURPOSE: The study of cell free DNA (cfDNA) enables sequential analysis of tumor cell-specific genetic alterations in neuroblastoma patients. EXPERIMENTAL DESIGN: Eighteen patients with relapsing neuroblastoma having received Lorlatinib, a 3rd generation ALK inhibitor, were identified (SACHA national registry and/or in the institution). cfDNA was analyzed at relapse for 9 patients, and sequentially for 5 patients (blood/bone marrow plasma) by performing WGS library construction followed by ALK-targeted ddPCR of the hotspot mutations (F1174L, R1275Q, I1170N) (variant allele fraction (VAF) detection limit 0.1%) and WES to evaluate disease burden and clonal evolution, following comparison with tumor/germline WES. RESULTS: Overall response rate to Lorlatinib was 33% (CI 13-59%), with response observed in 6/10 cases without versus 0/8 cases with MYCN amplification (MNA). ALK VAFs correlated with the overall clinical disease status, with a VAF<0.1% in clinical remission, versus higher VAFs (>30%) at progression. Importantly, sequential ALK ddPCR detected relapse earlier than clinical imaging. cfDNA WES revealed new SNVs, not seen in the primary tumor, in all instances of disease progression after Lorlatinib treatment, indicating clonal evolution, including alterations in genes linked to tumor aggressivity (TP53) or novel targets (EGFR). Gene pathway analysis revealed an enrichment for genes targeting cell differentiation in emerging clones, and cell adhesion in persistent clones. Evidence of clonal hematopoiesis could be observed in follow-up samples. CONCLUSION: We demonstrate the clinical utility of combining ALK cfDNA ddPCR for disease monitoring and cfDNA WES for the study of clonal evolution and resistance mechanisms in neuroblastoma patients receiving ALK targeted therapy.

CD74 supports accumulation and function of regulatory T cells in tumors.

Regulatory T cells (Tregs) are plastic cells playing a pivotal role in the maintenance of immune homeostasis. Tregs actively adapt to the microenvironment where they reside; as a consequence, their molecular and functional profiles differ among tissues and pathologies. In tumors, the features acquired by Tregs remains poorly characterized. Here, we observe that human tumor-infiltrating Tregs selectively overexpress CD74, the MHC class II invariant chain. CD74 has been previously described as a regulator of antigen-presenting cell biology, however its function in Tregs remains unknown. CD74 genetic deletion in human primary Tregs reveals that CD74KO Tregs exhibit major defects in the organization of their actin cytoskeleton and intracellular organelles. Additionally, intratumoral CD74KO Tregs show a decreased activation, a drop in Foxp3 expression, a low accumulation in the tumor, and consistently, they are associated with accelerated tumor rejection in preclinical models in female mice. These observations are unique to tumor conditions as, at steady state, CD74KO-Treg phenotype, survival, and suppressive capacity are unaffected in vitro and in vivo. CD74 therefore emerges as a specific regulator of tumor-infiltrating Tregs and as a target to interfere with Treg anti-tumor activity.

Deciphering the spatial landscape and plasticity of immunosuppressive fibroblasts in breast cancer.

Although heterogeneity of FAP+ Cancer-Associated Fibroblasts (CAF) has been described in breast cancer, their plasticity and spatial distribution remain poorly understood. Here, we analyze trajectory inference, deconvolute spatial transcriptomics at single-cell level and perform functional assays to generate a high-resolution integrated map of breast cancer (BC), with a focus on inflammatory and myofibroblastic (iCAF/myCAF) FAP+ CAF clusters. We identify 10 spatially-organized FAP+ CAF-related cellular niches, called EcoCellTypes, which are differentially localized within tumors. Consistent with their spatial organization, cancer cells drive the transition of detoxification-associated iCAF (Detox-iCAF) towards immunosuppressive extracellular matrix (ECM)-producing myCAF (ECM-myCAF) via a DPP4- and YAP-dependent mechanism. In turn, ECM-myCAF polarize TREM2+ macrophages, regulatory NK and T cells to induce immunosuppressive EcoCellTypes, while Detox-iCAF are associated with FOLR2+ macrophages in an immuno-protective EcoCellType. FAP+ CAF subpopulations accumulate differently according to the invasive BC status and predict invasive recurrence of ductal carcinoma in situ (DCIS), which could help in identifying low-risk DCIS patients eligible for therapeutic de-escalation.

The pelvic organs receive no parasympathetic innervation.

The pelvic organs (bladder, rectum, and sex organs) have been represented for a century as receiving autonomic innervation from two pathways - lumbar sympathetic and sacral parasympathetic - by way of a shared relay, the pelvic ganglion, conceived as an assemblage of sympathetic and parasympathetic neurons. Using single-cell RNA sequencing, we find that the mouse pelvic ganglion is made of four classes of neurons, distinct from both sympathetic and parasympathetic ones, albeit with a kinship to the former, but not the latter, through a complex genetic signature. We also show that spinal lumbar preganglionic neurons synapse in the pelvic ganglion onto equal numbers of noradrenergic and cholinergic cells, both of which therefore serve as sympathetic relays. Thus, the pelvic viscera receive no innervation from parasympathetic or typical sympathetic neurons, but instead from a divergent tail end of the sympathetic chains, in charge of its idiosyncratic functions.

Residual ANTXR1+ myofibroblasts after chemotherapy inhibit anti-tumor immunity via YAP1 signaling pathway.

Although cancer-associated fibroblast (CAF) heterogeneity is well-established, the impact of chemotherapy on CAF populations remains poorly understood. Here we address this question in high-grade serous ovarian cancer (HGSOC), in which we previously identified 4 CAF populations. While the global content in stroma increases in HGSOC after chemotherapy, the proportion of FAP+ CAF (also called CAF-S1) decreases. Still, maintenance of high residual CAF-S1 content after chemotherapy is associated with reduced CD8+ T lymphocyte density and poor patient prognosis, emphasizing the importance of CAF-S1 reduction upon treatment. Single cell analysis, spatial transcriptomics and immunohistochemistry reveal that the content in the ECM-producing ANTXR1+ CAF-S1 cluster (ECM-myCAF) is the most affected by chemotherapy. Moreover, functional assays demonstrate that ECM-myCAF isolated from HGSOC reduce CD8+ T-cell cytotoxicity through a Yes Associated Protein 1 (YAP1)-dependent mechanism. Thus, efficient inhibition after treatment of YAP1-signaling pathway in the ECM-myCAF cluster could enhance CD8+ T-cell cytotoxicity. Altogether, these data pave the way for therapy targeting YAP1 in ECM-myCAF in HGSOC.

Artificial intelligence-based pathology as a biomarker of sensitivity to atezolizumab-bevacizumab in patients with hepatocellular carcinoma: a multicentre retrospective study.

BACKGROUND: Clinical benefits of atezolizumab plus bevacizumab (atezolizumab-bevacizumab) are observed only in a subset of patients with hepatocellular carcinoma and the development of biomarkers is needed to improve therapeutic strategies. The atezolizumab-bevacizumab response signature (ABRS), assessed by molecular biology profiling techniques, has been shown to be associated with progression-free survival after treatment initiation. The primary objective of our study was to develop an artificial intelligence (AI) model able to estimate ABRS expression directly from histological slides, and to evaluate if model predictions were associated with progression-free survival. METHODS: In this multicentre retrospective study, we developed a model (ABRS-prediction; ABRS-P), which was derived from the previously published clustering-constrained attention multiple instance learning (or CLAM) pipeline. We trained the model fit for regression analysis using a multicentre dataset from The Cancer Genome Atlas (patients treated by surgical resection, n=336). The ABRS-P model was externally validated on two independent series of samples from patients with hepatocellular carcinoma (a surgical resection series, n=225; and a biopsy series, n=157). The predictive value of the model was further tested in a series of biopsy samples from a multicentre cohort of patients with hepatocellular carcinoma treated with atezolizumab-bevacizumab (n=122). All samples in the study were from adults (aged ≥18 years). The validation sets were sampled between Jan 1, 2008, to Jan 1, 2023. For the multicentre validation set, the primary objective was to assess the association of high versus low ABRS-P values, defined relative to cross-validation median split thresholds in the first biopsy series, with progression-free survival after treatment initiation. Finally, we performed spatial transcriptomics and matched prediction heatmaps with in situ expression profiles. FINDINGS: Of the 840 patients sampled, 641 (76%) were male and 199 (24%) were female. Across the development and validation datasets, hepatocellular carcinoma risk factors included alcohol intake, hepatitis B and C virus infections, and non-alcoholic steatohepatitis. Using cross-validation in the development series, the mean Pearson's correlation between ABRS-P values and ABRS score (mean expression of ABRS genes) was r=0·62 (SD 0·09; mean p<0·0001, SD<0·0001). The ABRS-P generalised well on the external validation series (surgical resection series, r=0·60 [95% CI 0·51-0·68], p<0·0001; biopsy series, r=0·53 [0·40-0·63], p<0·0001). In the 122 patients treated with atezolizumab-bevacizumab, those with ABRS-P-high tumours (n=74) showed significantly longer median progression-free survival than those with ABRS-P-low tumours (n=48) after treatment initiation (12 months [95% CI 7-not reached] vs 7 months [4-9]; p=0·014). Spatial transcriptomics showed significantly higher ABRS score, along with upregulation of various other immune effectors, in tumour areas with high ABRS-P values versus areas with low ABRS-P values. INTERPRETATION: Our study indicates that AI applied on hepatocellular carcinoma digital slides is able to serve as a biomarker for progression-free survival in patients treated with atezolizumab-bevacizumab. This approach could be used in the development of inexpensive and fast biomarkers for targeted therapies. The combination of AI heatmaps with spatial transcriptomics provides insight on the molecular features associated with predictions. This methodology could be applied to other cancers or diseases and improve understanding of the biological mechanisms that drive responses to treatments. FUNDING: Institut National du Cancer, Fondation ARC, China Scholarship Council, Ligue Contre le Cancer du Val de Marne, Fondation de l'Avenir, Ipsen, and Fondation Bristol Myers Squibb Pour la Recherche en Immuno-Oncologie.

Breast cancer remotely imposes a myeloid bias on haematopoietic stem cells by reprogramming the bone marrow niche.

Myeloid cell infiltration of solid tumours generally associates with poor patient prognosis and disease severity1-13. Therefore, understanding the regulation of myeloid cell differentiation during cancer is crucial to counteract their pro-tumourigenic role. Bone marrow (BM) haematopoiesis is a tightly regulated process for the production of all immune cells in accordance to tissue needs14. Myeloid cells differentiate during haematopoiesis from multipotent haematopoietic stem and progenitor cells (HSPCs)15-17. HSPCs can sense inflammatory signals from the periphery during infections18-21 or inflammatory disorders22-27. In these settings, HSPC expansion is associated with increased myeloid differentiation28,29. During carcinogenesis, the elevation of haematopoietic growth factors supports the expansion and differentiation of committed myeloid progenitors5,30. However, it is unclear whether cancer-related inflammation also triggers demand-adapted haematopoiesis at the level of multipotent HSPCs. In the BM, HSPCs reside within the haematopoietic niche which delivers HSC maintenance and differentiation cues31-35. Mesenchymal stem cells (MSCs) are a major cellular component of the BM niche and contribute to HSC homeostasis36-41. Modifications of MSCs in systemic disorders have been associated with HSC differentiation towards myeloid cells22,42. It is unknown if MSCs are regulated in the context of solid tumours and if their myeloid supportive activity is impacted by cancer-induced systemic changes. Here, using unbiased transcriptomic analysis and in situ imaging of HSCs and the BM niche during breast cancer, we show that both HSCs and MSCs are transcriptionally and spatially modified. We demonstrate that breast tumour can distantly remodel the cellular cross-talks in the BM niche leading to increased myelopoiesis.

KDM6B drives epigenetic reprogramming associated with lymphoid stromal cell early commitment and immune properties.

Mature lymphoid stromal cells (LSCs) are key organizers of immune responses within secondary lymphoid organs. Similarly, inflammation-driven tertiary lymphoid structures depend on immunofibroblasts producing lymphoid cytokines and chemokines. Recent studies have explored the origin and heterogeneity of LSC/immunofibroblasts, yet the molecular and epigenetic mechanisms involved in their commitment are still unknown. This study explored the transcriptomic and epigenetic reprogramming underlying LSC/immunofibroblast commitment. We identified the induction of lysine demethylase 6B (KDM6B) as the primary epigenetic driver of early immunofibroblast differentiation. In addition, we observed an enrichment for KDM6B gene signature in murine inflammatory fibroblasts and pathogenic stroma of patients with autoimmune diseases. Last, KDM6B was required for the acquisition of LSC/immunofibroblast functional properties, including the up-regulation of CCL2 and the resulting recruitment of monocytes. Overall, our results reveal epigenetic mechanisms that participate in the early commitment and immune properties of immunofibroblasts and support the use of epigenetic modifiers as fibroblast-targeting strategies in chronic inflammation.

Imaging and multi-omics datasets converge to define different neural progenitor origins for ATRT-SHH subgroups.

Atypical teratoid rhabdoid tumors (ATRT) are divided into MYC, TYR and SHH subgroups, suggesting diverse lineages of origin. Here, we investigate the imaging of human ATRT at diagnosis and the precise anatomic origin of brain tumors in the Rosa26-CreERT2::Smarcb1flox/flox model. This cross-species analysis points to an extra-cerebral origin for MYC tumors. Additionally, we clearly distinguish SHH ATRT emerging from the cerebellar anterior lobe (CAL) from those emerging from the basal ganglia (BG) and intra-ventricular (IV) regions. Molecular characteristics point to the midbrain-hindbrain boundary as the origin of CAL SHH ATRT, and to the ganglionic eminence as the origin of BG/IV SHH ATRT. Single-cell RNA sequencing on SHH ATRT supports these hypotheses. Trajectory analyses suggest that SMARCB1 loss induces a de-differentiation process mediated by repressors of the neuronal program such as REST, ID and the NOTCH pathway.

Molecular characterisation of tumours of the lacrimal apparatus.

AIMS: Malignant tumours of the lacrimal apparatus are rare and frequently show a poor prognosis, with no clear therapeutic standards. Characterisation of the genetic landscape of these rare tumours is sparse, and therefore therapeutics generally follow those of their common salivary gland counterparts. To further clarify the pathophysiology and discover potential therapeutic targets, we investigated the genetic landscape of eight tumours of the lacrimal apparatus. METHODS AND RESULTS: DNA and RNA sequencing were performed to identify genetic mutations and gene fusions. Immunohistochemistry, fluorescence in-situ hybridisation and reverse transcription-polymerase chain reaction followed by Sanger sequencing were performed to confirm the identified molecular alterations. Genetic alterations were detected in six tumours. Among five adenoid cystic carcinomas (ACC), four had confirmed alterations of MYB or MYBL1 genes, including a MYB::NFIB fusion, a MYBL1::NFIB fusion, a MYB amplification and a novel NFIB::THSD7B fusion. Mutations in genes encoding epigenetic modifiers, as well as NOTCH1, FGFR2 and ATM mutations, were also identified in ACCs. A carcinoma ex pleomorphic adenoma showed TP53 and CIC mutations and an amplification of ERBB2. A transitional cell carcinoma was associated with HPV16 infection. No genetic alteration was found for one adenocarcinoma, not otherwise specified. CONCLUSIONS: Our study highlights the variety of molecular alterations associated with lacrimal system tumours and emphasises the importance of molecular testing in these tumours, which can reveal potentially targetable mutations. Our results also reinforce the hypothesis of a common physiopathology of all ACCs, regardless of their primary location.

Multi-omics comparison of malignant and normal uveal melanocytes reveals molecular features of uveal melanoma.

Uveal melanoma (UM) is a rare cancer resulting from the transformation of melanocytes in the uveal tract. Integrative analysis has identified four molecular and clinical subsets of UM. To improve our molecular understanding of UM, we performed extensive multi-omics characterization comparing two aggressive UM patient-derived xenograft models with normal choroidal melanocytes, including DNA optical mapping, specific histone modifications, and DNA topology analysis using Hi-C. Our gene expression and cytogenetic analyses suggest that genomic instability is a hallmark of UM. We also identified a recurrent deletion in the BAP1 promoter resulting in loss of expression and associated with high risk of metastases in UM patients. Hi-C revealed chromatin topology changes associated with the upregulation of PRAME, an independent prognostic biomarker in UM, and a potential therapeutic target. Our findings illustrate how multi-omics approaches can improve our understanding of tumorigenesis and reveal two distinct mechanisms of gene expression dysregulation in UM.

Sweat Gland Tumors Arising on Acral Sites: A Molecular Survey.

Recurrent oncogenic drivers have been identified in a variety of sweat gland tumors. Recently, integration of human papillomavirus type 42 (HPV42) has been reported in digital papillary adenocarcinoma (DPA). The main objectives of the present study were (i) to provide an overview of the prevalence of previously identified oncogenic drivers in acral sweat gland tumors and (ii) to genetically characterize tumors in which no recurrent genetic alteration has been identified yet. Cases of acral sweat gland tumors were identified from the database of the French network CARADERM. After histologic review, the presence of previously identified genetic alterations was investigated in the entire cohort (n=79) using a combination of immunohistochemistry and targeted DNA and RNA sequencing. Tumor entities with no recurrent genetic alterations were submitted to whole-transcriptome sequencing. CRTC1::MAML2 fusion was identified in cases of hidradenoma and hidradenocarcinoma (n=9/12 and n=9/12). A p.V600E mutation of BRAF was observed in all cases of tubular adenoma (n=4). YAP1:MAML2 and YAP1::NUTM1 fusions were observed in poroid tumors (n=15/25). ETV6::NTRK3 and TRPS1::PLAG1 fusion transcripts were identified in secretory carcinoma (n=1/1) and cutaneous mixed tumors (n=3/4), respectively. The HPV42 genome was detected in most cases of DPA (n=10/11) and in 1 adnexal adenocarcinoma not otherwise specified. Finally, whole-transcriptome analysis revealed BRD3::NUTM1 or NSD3::NUTM1 fusions in 2 cases of NUT adnexal carcinoma and NCOA4::RET and CCDC6::RET fusion transcripts in 2 cystadenoma/hidrocystoma-like tumors. Our study confirms distinctive cytogenetic abnormalities in a wide number of acral adnexal neoplasms and supports the use of molecular analysis as a valuable aid in the diagnosis of these rare and often difficult to diagnose group of neoplasms.

Cell-Free DNA Extracted from CSF for the Molecular Diagnosis of Pediatric Embryonal Brain Tumors.

BACKGROUND: Liquid biopsies are revolutionary tools used to detect tumor-specific genetic alterations in body fluids, including the use of cell-free DNA (cfDNA) for molecular diagnosis in cancer patients. In brain tumors, cerebrospinal fluid (CSF) cfDNA might be more informative than plasma cfDNA. Here, we assess the use of CSF cfDNA in pediatric embryonal brain tumors (EBT) for molecular diagnosis. METHODS: The CSF cfDNA of pediatric patients with medulloblastoma (n = 18), ATRT (n = 3), ETMR (n = 1), CNS NB FOXR2 (n = 2) and pediatric EBT NOS (n = 1) (mean cfDNA concentration 48 ng/mL; range 4-442 ng/mL) and matched tumor genomic DNA were sequenced by WES and/or a targeted sequencing approach to determine single-nucleotide variations (SNVs) and copy number alterations (CNA). A specific capture covering transcription start sites (TSS) of genes of interest was also used for nucleosome footprinting in CSF cfDNA. RESULTS: 15/25 CSF cfDNA samples yielded informative results, with informative CNA and SNVs in 11 and 15 cases, respectively. For cases with paired tumor and CSF cfDNA WES (n = 15), a mean of 83 (range 1-160) shared SNVs were observed, including SNVs in classical medulloblastoma genes such as SMO and KMT2D. Interestingly, tumor-specific SNVs (mean 18; range 1-62) or CSF-specific SNVs (mean 5; range 0-25) were also observed, suggesting clonal heterogeneity. The TSS panel resulted in differential coverage profiles across all 112 studied genes in 7 cases, indicating distinct promoter accessibility. CONCLUSION: CSF cfDNA sequencing yielded informative results in 60% (15/25) of all cases, with informative results in 83% (15/18) of all cases analyzed by WES. These results pave the way for the implementation of these novel approaches for molecular diagnosis and minimal residual disease monitoring.

Oxidative phosphorylation is a metabolic vulnerability of endocrine therapy and palbociclib resistant metastatic breast cancers.

Resistance to endocrine treatments and CDK4/6 inhibitors is considered a near-inevitability in most patients with estrogen receptor positive breast cancers (ER + BC). By genomic and metabolomics analyses of patients' tumours, metastasis-derived patient-derived xenografts (PDX) and isogenic cell lines we demonstrate that a fraction of metastatic ER + BC is highly reliant on oxidative phosphorylation (OXPHOS). Treatment by the OXPHOS inhibitor IACS-010759 strongly inhibits tumour growth in multiple endocrine and palbociclib resistant PDX. Mutations in the PIK3CA/AKT1 genes are significantly associated with response to IACS-010759. At the metabolic level, in vivo response to IACS-010759 is associated with decreased levels of metabolites of the glutathione, glycogen and pentose phosphate pathways in treated tumours. In vitro, endocrine and palbociclib resistant cells show increased OXPHOS dependency and increased ROS levels upon IACS-010759 treatment. Finally, in ER + BC patients, high expression of OXPHOS associated genes predict poor prognosis. In conclusion, these results identify OXPHOS as a promising target for treatment resistant ER + BC patients.

Reversible transitions between noradrenergic and mesenchymal tumor identities define cell plasticity in neuroblastoma.

Noradrenergic and mesenchymal identities have been characterized in neuroblastoma cell lines according to their epigenetic landscapes and core regulatory circuitries. However, their relationship and relative contribution in patient tumors remain poorly defined. We now document spontaneous and reversible plasticity between the two identities, associated with epigenetic reprogramming, in several neuroblastoma models. Interestingly, xenografts with cells from each identity eventually harbor a noradrenergic phenotype suggesting that the microenvironment provides a powerful pressure towards this phenotype. Accordingly, such a noradrenergic cell identity is systematically observed in single-cell RNA-seq of 18 tumor biopsies and 15 PDX models. Yet, a subpopulation of these noradrenergic tumor cells presents with mesenchymal features that are shared with plasticity models, indicating that the plasticity described in these models has relevance in neuroblastoma patients. This work therefore emphasizes that intrinsic plasticity properties of neuroblastoma cells are dependent upon external cues of the environment to drive cell identity.

A druggable copper-signalling pathway that drives inflammation.

Inflammation is a complex physiological process triggered in response to harmful stimuli1. It involves cells of the immune system capable of clearing sources of injury and damaged tissues. Excessive inflammation can occur as a result of infection and is a hallmark of several diseases2-4. The molecular bases underlying inflammatory responses are not fully understood. Here we show that the cell surface glycoprotein CD44, which marks the acquisition of distinct cell phenotypes in the context of development, immunity and cancer progression, mediates the uptake of metals including copper. We identify a pool of chemically reactive copper(II) in mitochondria of inflammatory macrophages that catalyses NAD(H) redox cycling by activating hydrogen peroxide. Maintenance of NAD+ enables metabolic and epigenetic programming towards the inflammatory state. Targeting mitochondrial copper(II) with supformin (LCC-12), a rationally designed dimer of metformin, induces a reduction of the NAD(H) pool, leading to metabolic and epigenetic states that oppose macrophage activation. LCC-12 interferes with cell plasticity in other settings and reduces inflammation in mouse models of bacterial and viral infections. Our work highlights the central role of copper as a regulator of cell plasticity and unveils a therapeutic strategy based on metabolic reprogramming and the control of epigenetic cell states.

The mutational landscape of skull base and spinal chordomas and the identification of potential prognostic and theranostic biomarkers.

OBJECTIVE: Chordomas are rare bone neoplasms characterized by a high recurrence rate and no benefit from any approved medical treatment to date. However, the investigation of molecular alterations in chordomas could be essential to prognosticate, guide clinical decision-making, and identify theranostic biomarkers. The aim of this study was to provide a detailed genomic landscape of a homogeneous series of 64 chordoma samples, revealing driver events, theranostic markers, and outcome-related genomic features. METHODS: The authors conducted whole-exome sequencing (WES), targeted next-generation sequencing, and RNA sequencing of 64 skull base and spinal chordoma samples collected between December 2006 and September 2020. Clinical, histological, and radiological data were retrospectively analyzed and correlated to genetic findings. RESULTS: The authors identified homozygous deletions of CDKN2A/2B, PIK3CA mutations, and alterations affecting genes of SWI/SNF chromatin remodeling complexes (PBRM1 and ARID1A) as potential theranostic biomarkers. Using matched germline WES, they observed a higher frequency of a common genetic variant (rs2305089; p.(Gly177Asp)) in TBXT (97.8%, p < 0.001) compared to its distribution in the general population. PIK3CA mutation was identified as an independent biomarker of short progression-free survival (HR 10.68, p = 0.0008). Loss of CDKN2A/2B was more frequently observed in spinal tumors and recurrent tumors. CONCLUSIONS: In the current study, the authors identified driver events such as PBRM1 and PIK3CA mutations, TBXT alterations, or homozygous deletions of CDKN2A/2B, which could, for some, be considered potential theranostic markers and could allow for identifying novel therapeutic approaches. With the aim of a future biomolecular prognostication classification, alterations affecting PIK3CA and CDKN2A/2B could be considered as poor prognostic biomarkers.