- Open Access
Proteo-transcriptomics meta-analysis identifies SUMO2 as a promising target in glioblastoma multiforme therapeutics
Cancer Cell International volume 21, Article number: 575 (2021)
Glioblastoma multiforme (GBM) is a deadly brain tumour with minimal survival rates due to the ever-expanding heterogeneity, chemo and radioresistance. Kinases are known to crucially drive GBM pathology; however, a rationale therapeutic combination that can simultaneously inhibit multiple kinases has not yet emerged successfully.
Here, we analyzed the GBM patient data from several publicly available repositories and deduced hub GBM kinases, most of which were identified to be SUMOylated by SUMO2/3 isoforms. Not only the hub kinases but a significant proportion of GBM upregulated genes involved in proliferation, metastasis, invasion, epithelial-mesenchymal transition, stemness, DNA repair, stromal and macrophages maintenance were also identified to be the targets of SUMO2 isoform. Correlatively, high expression of SUMO2 isoform was found to be significantly associated with poor patient survival.
Although many natural products and drugs are evidenced to target general SUMOylation, however, our meta-analysis strongly calls for the need to design SUMO2/3 or even better SUMO2 specific inhibitors and also explore the SUMO2 transcription inhibitors for universally potential, physiologically non-toxic anti-GBM drug therapy.
The major highlights of this study are as follows:
Key upregulated hub kinases and coding genes in GBM are found to be targets of SUMO2 conjugation.
SUMO2 is significantly expressed in adult primary and recurrent GBMs as well as in pediatric GBM tumours.
Orthotropic xenografts from adult and pediatric GBMs confirm high expression of SUMO2 in GBM tumour samples.
SUMO2 is significantly associated with patient survival plot and pan-cancer cell fitness.
Rationale design of SUMO2 inhibitors or search for its transcriptional inhibitors is urgently required through industry-academia collaboration for an anti-GBM and potentially pan-cancer therapeutics.
Glioblastoma multiforme (GBM) is an astrocyte-derived brain tumour with minimal survival rates. It is the deadliest of all cancers thwarted with dismal hopes of any rationale therapeutics due to the scale and complexity of pathology across ages, genders, ethnicities and intracellular plasticity . While basic laboratory research through the “in vitro” cell culture studies and “in vivo” animal model screening of candidate drugs can lead to anti-GBM drug target discovery, such endeavours are most often not translatable in the clinical settings. The excavation of cancer cell-specific vulnerabilities and their progression into clinical utility is thus essentially required for GBM management. In this direction, kinases are identified as central players in the formation, maintenance and recurrence of aggressive GBM tumours . Both receptor kinases and non-receptor kinases are expressed in GBMs such as EGFR/mutant EGFRvIII, IGFR, PDGFR, VEGFR, FGFR, TGFβR, mutant PTEN, PI3K, AKT, MAPK, mTOR, Fyn, Src kinases (c-Src, Yes, Lyn and Lck), GSK3, Serum glucocorticoid kinases (SGK), ATM, casein kinase 2 (CK2), focal adhesion kinases (FAKs), sphingosine kinases (SPHK), Ephrins, polo-like kinases (PLK1), Aurora kinases (AURKA). Cyclin-dependent kinases (CDKs), LIMK1/2, ROCK1/2, PKD1/2, PAKs, to name a few [1,2,3].
Post-translational modifications in kinases crucially regulate their activity. Since several kinases are activated by phosphorylation, the attempts are on to target phosphorylation sites. However, multiple phosphorylation sites (identified via phosphoproteomics) and de novo mutations in phosphorylation sites are posing enormous difficulties in identifying druggable site-specific inhibitors . Besides, the rapid emergence of various somatic variants of these kinases brings forth another level of complexity in drug target site identification .
Moreover, multiple kinases need to be targeted simultaneously because there is a redundancy in their functions. For example, the PI3K/AKT/mTOR (PAM) pathway contributes significantly to oncogenesis and tumour progression. Multiple receptor tyrosine kinases such as EGFR/ EGFRvIII/IGFR/PDGFR generate redundant activation of phosphoinositide-3′-kinase (PI3K) signalling . Similarly, inhibition of mTOR signaling transactivates MAPK, which is linked to PI3K activation and can produce redundant functions [7, 8].
Other than multisite-phosphorylations, SUMOylation of kinases is crucially observed in various pathologies, especially in cancers. SUMO-conjugation of proteins has been evidenced to help cancer cells cope with microenvironment and metabolic stresses, thereby enabling survival [9, 10]. GBMs are evidenced to have prominent SUMOylation of both receptor and non-receptor kinases and other protumourigenic coding genes [11,12,13]. SUMOylation of CDK6, various cyclins, PI3K, AKT, ERK5 drive cell cycle progression and proliferation in GBMs . SUMOylation of several cytoskeletal proteins such as vimentin and collapsin response mediator protein 2 (CRMP2) are also reported to promote over-proliferation [14, 15]. SUMOylation of ERK5, IGFR and EGFR allows their trafficking to the nucleus, where they upregulate pro-proliferative cyclin D1 and beta-catenin transcription factors [16,17,18]. Many G2/M checkpoint kinases such as BUB1B/BubR1, MPS1/TTK, Aurora-B, polo-like kinase 1 (PLK) are modified by SUMOylation, only upon which they serve to promote mitotic progression and cell division . SUMOylation of Lats1 and AMPK antagonizes their tumour-suppressor functions [19, 20], whereas SUMO modification of NPM-ALK and Erythropoietin-producing hepatocellular carcinoma (Eph) family of kinases promote anchorage-independent growth .
However, since SUMOylation is also essential for maintaining cell homeostasis, a general inhibition of SUMOylation is not desirable. We, therefore, wanted to explore, in a clinical setting, the relevance of targeting a single SUMO isoform in suppressing multiple protumourigenic kinases and other genes vs global SUMOylation. This is required because general SUMOylation inhibitors like ginkgolic and anacardic acid fail to decrease global SUMOylation in GBM cell lines .
Meta-analysis and research weaving of clinical cancer data available in different public repositories and publications enables connecting the dots and bringing together an innovative description of cancer cell druggable vulnerabilities than is usually possible in a single study . Hence, these computational approaches can allow unique opportunities to prioritise “go versus no-go” decisions in developing bench-to-bedside approaches to conquer cancer cells. Therefore, we were motivated to utilize these bioinformatics approaches to explore if any SUMO isoform post-translation modifications can be targeted, in place of phosphorylation, to suppress the pro-tumorigenic functions of multiple GBM kinases.
Materials and methods
The Additional Tables/Additional files section provides detailed information on additional methods, individual data points, analyzed datasets, and website links associated with bioinformatics data analysis.
In this study, a non-redundant list of 723 kinases was created from various literature sources [22, 23]. The level 3 gene transcriptome profile data (RNAseq-HTSeq-Count) of 173 clinical samples of GDC TCGA Glioblastoma and associated metadata files were downloaded from the UCSC XENA browser (https://xenabrowser.net/datapages/?dataset=TCGA-GBM.htseq_counts.tsv&host=https%3A%2F%2Fgdc.xenahubs.net&removeHub=https%3A%2F%2Fxena.treehouse.gi.ucsc.edu%3A443) 07-19-2019 version . UCSC XENA has mRNA HTSeq counts normalized across the samples and is log2(Count + 1) transformed. Samples of solid normal tissue, primary tumour, and recurrent tumour were chosen for analysis (https://gdc-hub.s3.us-east-1.amazonaws.com/download/TCGA-GBM.GDC_phenotype.tsv.gz; https://docs.gdc.cancer.gov/Data_Dictionary/viewer/#?view=table-entity-listanchor=clinical). Please see Additional file 1: Table S1 for clinical data, and Additional file 2: Table S1 for mRNA HTSeq counts data of patient samples. Raw count matrix was obtained by back transforming the log2(Count + 1) data using R code for differential expression analysis (Additional file 3: Table S1). R code for back conversion is provided in the additional methods.
Identifying differentially expressed genes (DEGs)
Differential expression analysis was carried out with the DESeq2 package (version 1.24.0) in R Studio. Due to the high sensitivity and precision offered by the DESeq2 method, it is used for comparative analysis of transcriptomics data between test and control samples to obtain differentially expressed genes . Rows with low count genes (counts < 10) were pre-filtered to reduce the memory size and to increase the speed of transformation and testing functions within DESeq2. Genes with Padj < 0.05 and FC > = 2 were considered to be statistically significant DEGs. Differential expression analysis was performed for 3 sample types and three conditions, including primary tumour vs solid tissue normal, recurrent tumour vs solid tissue normal, and primary tumour vs recurrent tumour. From the kinases genes list, differentially expressed genes in all three conditions were manually curated.
Visualizing DEG’s with Heatmap and Volcano plot
Significantly expressed kinases in the three conditions (primary GBM vs adjacent normal samples, recurrent GBM vs adjacent normal samples and primary GBM vs recurrent GBM samples) were categorized into up-regulated (log2FC > = 1) and down-regulated (log2FC < = − 1) genes. Heatmap visualization of the high and low expressed kinases in all three conditions and the volcano plots were created using the Morpheus tool from Broad Institute (https://software.broadinstitute.org/morpheus/) and Graphpad Prism (version6.01), respectively.
Protein–protein interaction (PPI) network construction and hub genes selection
The Search Tool for the Retrieval of Interacting Genes (STRING version 11.0, https://string-db.org/) database is used to construct a PPI network with the predicted association for a group of genes . String predicts protein–protein association, which includes both physical and functional interactions. It provides the protein functional association ranked by a confidence score. Cytoscape software (version3.7.1) was used for PPI graphical network visualization. In the network, each node represents a gene/protein and edges represent the connection between them. Cytoscape is a freely available software package, majorly used for visualizing and analysing molecular & genetic interaction networks .Cytoscape support many algorithms for the network layout representation and many plugins for further network analysis. We used an Edge-weighted Spring Embedded layout for our network and Network Analyzer to efficiently compute the topological network parameters like degree and betweenness centrality .
Identifying hub genes from the PPI network
Cytohubba (version 0.1) plugin was used to get the top 30 hubs genes and sub-networks from the PPI network . Cytohubba uses different algorithms to identify the sub-network of the hub genes from a more extensive network, which includes Density of Maximum Neighbourhood Component (DMNC), Maximum Neighbourhood Component (MNC), Maximal Clique Centrality (MCC) and degree. These methods are local based methods that consider the neighborhood of a node. The result from DMNC methods, however, did not show consistency when compared to the other three methods. So the top 30 genes obtained from three methods MNC, MCC and degree, were overlapped, and the overlapping set of genes were considered as major hubs.
Gene set enrichment analysis
Gene set enrichment analysis was performed using GSEA software from Broad institute (version 4.0.3) [https://www.gsea-msigdb.org/gsea/index.jsp] . GSEA Pre ranked method was used to understand the enriched pathways and functional annotation of ranked differentially expressed genes. Upregulated and downregulated pre-ranked kinase genes from DESeq2 analysis was used for GSEA analysis. Curated geneset databases “c2.all.v7.1.symbols.gmt”and “Human_Symbol_with_Remapping_MSigDB.v7.1.chip” were set as chip platform for pathway enrichment analysis, also “c8.all.v7.2.symbols” [Cell type signature gene sets] was used to retrieve curated cluster markers for cell types identified in single-cell sequencing studies of human tissues using GSEA software (version 4.0.3). After 1000 permutations, Normalized enrichment score (NES) was calculated and gene set with p-value < 0.05, FDR < 0.25 were considered significant. Functional enrichment analysis was performed using DAVID [version 6.8, https://david.ncifcrf.gov/] ], and gene ontology over-representation analysis was performed with clusterProfiler [version 3.18] .
Survival analysis was performed using the PrognoScan database [http://www.prognoscan.org/; http://dna00.bio.kyutech.ac.jp/PrognoScan/] . The correlation between SUMO isoforms expression and survival in brain cancer patients were analysed, using PrognoScan and the Kaplan–Meier plot. The survival curves were plotted for high expression (red) and low expression (blue) groups dichotomized at the optimal cut-point. The significant corrected p-value threshold was adjusted at < 0.05.
Cell fitness analysis
Pan-cancer SUMO isoforms associated with cell fitness data were collected from the Cancer Dependency Map Dataset [https://depmap.org/portal/depmap/] . It comprises data of genome-wide CRISPR and shRNA screens to identify essential genes across hundreds of cancer cell lines.
Protein expression analysis from the human protein atlas (HPA)
The human protein atlas (https://www.proteinatlas.org/) provides a complete resource of proteomics and transcriptomics data generated from antibody-based microarray profiling and RNA sequencing. It can be majorly used to study protein co-localization and expression in human tissues and cells. HPAnalyze version 3.12, a freely available R package, was used to retrieve and visualize the data from the human protein atlas. hpaVisPatho() function from HPAnalyze was used to visualize the expression of the protein of interest in each cancer .
Other data download web links and data reference IDs used in the study
(i) UCSC Xena browser to compare TCGA tumour samples to GTEx normal samples to see if our gene or transcript is up- or down-regulated in one or more cancer types [https://xena.ucsc.edu/compare-tissue/]; (ii) Cancer Cell Line Encyclopedia (CCLE) [https://portals.broadinstitute.org/ccle]; (iii) RNA Seq data on subtype-specific differences in molecular and cellular composition at the margins of glioblastoma [Ref id: GEO-GSE59612, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE59612, PMID-25114226]; (iv) microarray data from GBM samples of patients showing radioresistance and chemoresistance [Ref id: GEO-GSE7696, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE7696, PMID-18565887, 21642372]; (iv) microarray data from the human glioblastoma cell culture resource (HGCC) [Ref id: GEO-GSE72217, https://www.hgcc.se/,PMID-26629530]; (vi) microarray data from patient-derived cell line and xenograft models of proneural, classical and mesenchymal glioblastoma [Ref id: GEO-GSE118793, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE118793; Ref id:, SRA-PRJNA508446, https://www.ncbi.nlm.nih.gov/bioproject/PRJNA508446/]; (vii) pediatric high and low grade glioma; [CBTTC,https://xenabrowser.net/datapages/?cohort=Pediatric%20Brain%20Tumour%20Atlas%3A%20CBTTC&removeHub=https%3A%2F% 2Fxena.treehouse.gi.ucsc.edu%3A443]; (viii) adult high and low grade glioma [UCSC Xena, High grade glioma:https://gdc-hub.s3.us-east-1.amazonaws.com/download/TCGA-GBM.htseq_counts.tsv.gz; Low grade glioma:https://gdc-hub.s3.us-east-1.amazonaws.com/download/TCGA-LGG.htseq_counts.tsv.gz]; (ix) microarray data for gene expression in immune cells [https://joycelab.shinyapps.io/braintime/]; (x) GBM RNA Seq data on cancer stem cells from IVYGAP [https://glioblastoma.alleninstitute.org/, https://glioblastoma.alleninstitute.org/api/v2/well_known_file_download/305873915]; (xi) protein mass spectrometry data from PDC commons [https://proteomic.datacommons.cancer.gov/pdc/https://proteomic.datacommons.cancer.gov/pdc/study/PDC000204]; (xii) RNA Seq data of developing/prenatal human brain [https://portal.brain-map.org/http://www.brainspan.org/static/download.html]; (xiii) microarray data from patient-derived pediatric brain tumour cell lines and tumour animal models [Ref id: GEO-GSE99961, https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE99961].
SUMO2 putatively modify upregulated kinases and coding genes to enable glioblastoma multiforme (GBM) pathology
Towards the inception of the multi-kinase targeting strategy for GBM, we first compiled a non-redundant list of 723 human kinases from various literature sources (Additional file 4: Table S1).
Next, to identify the expression of these kinases in GBM patients' tissue, the RNA-Seq samples from 168 GBM tumours (primary GBM n = 155, recurrent GBM n = 13) and 5 normal brain samples were downloaded from the UCSC XENA browser. The clinical data demonstrating the heterogeneity of the GBM samples is provided in the Additional file 1: Table S1.
Samples of solid normal tissues (non tumour adjacent reference sample), primary GBM tumours, and recurrent GBM tumours were chosen for differential transcriptomics analysis of GBM (Additional file 4: Table S2). Fold change FC ≥ 2, i.e., up-regulated (log2FC > = 1) and down-regulated (log2FC < = − 1) genes were considered to be statistically significant to the reference samples for all genes and kinases expression analysis (Additional file 4: Table S2-S4, Additional file 5: Figures S1, S2).
A total of 114 kinases were identified to be upregulated in primary tumours and 125 in recurrent tumours in comparison to the adjacent normal tissues (Fig. 1A and Additional file 4: Table S5). A comparison of upregulated kinases in primary GBMs and recurrent GBMs showed 102 commonly upregulated kinases (Fig. 1A, Additional file 4: Table S5). Functional gene set enrichment analysis (GSEA) of the upregulated kinases in primary and recurrent GBM samples were intriguingly associated with cell cycle and mitosis promoting processes (Fig. 1B and Additional file 4: Table S6). The top 30 differentially upregulated hub kinases were then screened in primary and recurrent samples (vs normal adjacent samples) by combining the three local‑based methods (MNC, MCC and degree) in the Cytoscape plugin cytoHubba (Additional file 4: Table S7). From the results obtained, 25 hub upregulated kinases were common in primary and recurrent GBM samples amongst the top 30 hubs obtained by individual methods (Fig. 1C and Additional file 4: Table S7).
RNA-Seq analysis showed high SUMO1, 2 and 3 isoforms in both primary and recurrent GBM tissue samples compared to unmatched normal brain control samples and tumour adjacent normal samples (Fig. 1D, Additional file 4: Table S8). We found that 88 percent [22 out of 25] of common hub kinases between primary and recurrent GBMs bear a capacity to be SUMOylated predominantly by SUMO2/3 isoforms (combined inferences from literature survey and mass-spectrometry studies) [Fig. 1C and Additional file 4: Table S9–11] .
Mass spectrometry data extraction confirmed 76% (19 out of 25) kinases as the targets of SUMO2 isoforms [Fig. 1C, Additional file 4: Table S12] . Not only hubs but overall, approximately more than 30% percent of upregulated kinases in primary and recurrent GBM samples showed potential to be modified by SUMO2/3 via mass-spectrometry based analysis (Additional file 5: Figure S3A and Additional file 4: Table S13).
Not just kinases, approximately 19% of all upregulated coding genes were identified to bear a potential to be modified by SUMO2/3 isoforms (Additional file 4: Table S14 and Additional file 5: Figure S3B). SUMO2 directed mass spectrometry data showed 36% of all upregulated kinases and 25% of all upregulated coding genes to be the targets of SUMO2 conjugation (Fig. 1E and F, Additional file 4: Table S15 and S16).
The GSEA analysis of all upregulated GBM (primary and recurrent) genes putatively modified by SUMO2/3, showed enrichment for major pro-tumourigenic processes such as invasion, metastasis, proliferation, stemness, EMT, drug resistance, hypoxia, plasticity, downregulation of senescence, adipogenesis, upregulation of protumourigenic genes associated with HOX, Myc, MLL, HNFIA, NPM1, KRT19 and RhoA and enrichment of inflammatory stromal and microglial cells (Fig. 1G and Additional file 4: Table S17).
Besides, supplementary platforms for identification of gene enrichment processes such as DAVID GO Clusters and Over-representation tests showed an increase in DNA repair, protumourigenic sonic hedgehog, src kinase, Hox9, cyclophilin, oxidative stress response and folic acid metabolism pathways (Additional file 4: Table S18, 19, Additional file 5: Figures S4 and S5). Hence, broadly putative SUMO2/3 modification of kinases and other upregulated genes could be strongly associated with primary and recurrent GBM maintenance and progression.
Indeed, even the proteomics data showed high upregulation of SUMO2 and SUMO2/3 as analyzed through Protein Data Commons and Human protein Atlas respectively (Additional file 4: Tables S20, 21, Additional file 5: Figure S6).
SUMO2 is significantly expressed in heterogenous cell types of GBM and is associated with poor prognosis
Analysis of RNA-Seq data by Gill et al.  [GSE59612] further showed that SUMO2 isoform was maximally upregulated in GBM core tissue and was also observed in tumour margins (Fig. 2A and Additional file 4: Table S22). Besides, pro-tumourigenic cancer stem cells (Fig. 2B, Additional file 4: Table S23), the tumour-promoting aberrant immune cells such as microglia and monocyte-derived macrophages , were also identified to express high levels of SUMO2 (Additional file 5: Figure S7 and Additional file 4: Table S24). Hence, SUMO2 was expressed throughout tumour heterogeneity.
Since the most aggressive GBMs arise from or contact the neural stem cell enriched lateral ventricular zones; analysis of the RNA-Seq data submitted by Mistry et al.  showed high expression of SUMO isoforms (except the SUMO4) in GBMs of both ventricular-sub ventricular neural stem cell regions (VSVZ +) and those that did not contact the ventricular stem cell linings (VSVZ-)(Fig. 2C and Additional file 4: Table S25).
SUMO2/3 high expression was also found to be predominantly associated with both chemo and radioresistant GBM samples, as per the analysis of microarray data submitted by Murat et al. [GSE7696]  (Fig. 2D and Additional file 4: Table S26), but amongst the two isoforms only SUMO2 was found to be significantly associated with the GBM patient reduced survival (Fig. 2E and Additional file 4: Table S27). This information suggests that identification of SUMO2 inhibitors and inhibition processes, which can target multiple GBM kinases/other upregulated pro-tumourigenic genes and render them non-functional or even degrade them, needs to be urgently initiated.
So at this juncture, we found it essential to consider whether a background toolbox is available for SUMO2 directed drug discovery against GBM. In order to facilitate the drug screening research, we found that both freshly derived patient GBM cell lines from HGCC (GEO-GSE72217, https://www.hgcc.se/, available for research) and other sources [Stringer et al., SRA-PRJNA508446]  as well as classically used GBM tumour cell lines (https://portals.broadinstitute.org/ccle), consistently upregulated transcriptome of SUMO2 (Fig. 2F and Additional file 4: Table S28), hence can be alternately used for anti-SUMO2 drug discovery.
The transcriptome of GBM patient tissue (TCGA) as well as patient derived cell lines (HGCC) under various molecularly distinct subtypes (classical, mesenchymal, neural and proneural) showed enhanced SUMO2 expression (Fig. 2G, H and Additional file 4: Table S29, 30).
SUMO2 isoform is significantly expressed in both adult and pediatric GBMs and is essential for cancer cell fitness
It is also worth noting that microarray-based transcriptomics analysis of GBM subtype tumours, freshly derived tumour cells from these cells and the respective orthotropic xenografts obtained from the freshly derived tumour cells, all showed similar high expression of SUMO2 [Stringer et al., GEO-GSE118793] (Fig. 3A–C and Additional file 4: Table S31).
Although pediatric GBM landscape is considered to be very different from the adult GBMs and the research is focused on finding different drugs for the two GBM types, we found that both pediatric low and high-grade gliomas expressed high levels of SUMO2 (Pediatric Brain Tumour Atlas: CBTTC from UCSC Xena) (Fig. 3D and Additional file 4: Table S32). Moreover, this was consistent with observed profiles in the adult primary and recurrent low and high gliomas (TCGA and CBTTC) (Additional file 5: Figure S8 and Additional file 4: Table S33). Besides, freshly derived cells from pediatric GBM tumours and orthotropic xenografts generated from them also showed consistently high SUMO2 expression [Fig. 3 E–G, cell lines are available for research; Additional file 4: Table S34] .
SUMOylation of cell cycle promoting EGFR kinase and SUMO2/3 involvement in astrocytic tumours have been documented; however, these were never considered or escalated seriously on translational platforms to become a part of GBM clinical management practice[12, 43].
There are some FDA-approved drugs (kerriamycin B, Spectinomycin B, topotecan) and plant-derived compounds (ginkgolic acid, tannic acid, anacardic acid) that target the SUMOylation process in general. Several other small molecule inhibitors are under development (GSK145A, 2- D08; SUMO deconjugating SENP1 inhibitors such as triptolide, momordine, compound J5, compound 4, compound 3, and compound 13 m) . The use of these compounds in the clinical management of GBM has not materialized as general inhibition of SUMOylation over an extended period of administration can have several overt effects on the physiology and functioning of other organs .
GBMs are evidenced to recapitulate developmental processes, and indeed, even though SUMO-2 and SUMO-3 display ~ 95% sequence homology, SUMO-2 expression was significantly higher than SUMO-3 during development as in GBMs [Additional file 5: Figure S9A] . Correlatively, it is also well documented that SUMO-2 knockout mice are embryonically lethal while SUMO-3 knockout mice are phenotypically normal .
As far as SUMO1 isoform is concerned, it has ~ 50% homology with SUMO-2/3 but has a lower expression than SUMO2 in developing mouse and human brain (Additional file 5: Figure S9A, Additional file 5: Figure S9B and Additional file 4: Table S35). Developmentally, SUMO2/3 can compensate for the loss of SUMO1; hence SUMO1 knockout mice are also viable. This suggests that even in GBMs, where the microenvironment recapitulates early brain development, SUMO1 may be overall less significant than SUMO2.
Indeed, although both SUMO1 and SUMO2 have been identified to be expressed in glioma stem cells (Fig. 2B) [44, 45], and since human glioblastoma cells reflect a return to a foetal/embryonic developmental cell state , SUMO2 being indispensable to the development but not in adult physiology is plausibly a more potential target in anti-GBM therapeutics than SUMO1 in targeting glioblastoma stem cells.
In glioma/glioblastomas, even at protein levels, SUMO2 is found to be significantly highly expressed than SUMO1 (Additional file 5: Figure S6). This further indicates that SUMO2 alone targeting via intelligent structural design or by better understanding its transcriptional process, or even exploring the possibility of genetic manipulation of the SUMO2 translation pathway via microRNAs must be explored. There is also a possibility that even with very high sequence homology of SUMO2 with SUMO3, their conformational folding and transcriptional stimuli may vary, assigning them under distinct spatio-temporal regulatory controls, which may further attribute these isoforms with different functional roles.
SUMO2 has also been confirmed to be majorly upregulated in various other cancer cell lines (Additional file 5: Figure S10). Interestingly, CRISPR/CAS9 high throughput screening databases (Fig. 3H, www.depMap.org) that allow gene-specific cancer cell fitness assessment have shown SUMO2 to be a pan-cancer fitness gene vs other isoforms. Hence, the use of SUMO2 inhibition strategies merits investigation in both high and low-grade astrocytomas/gliomas. In addition, SUMO2 may be a crucial drug target in other tumours too, as indicated by patients' survival plots (Additional file 5: Figure S11, S12, Additional file 4: Table S36).
We took proteo-transcriptomics approach to identify unique and promising drug target in GBM pathology and uncovered that (1) Key upregulated hub kinases and coding genes in GBM pathology controlling over-proliferation, hypoxia, stemness, ECM stiffening, inflammation, EMT, invasion, immune escape, chemo and radioresistance, are putative targets of SUMO2 conjugation, (2) SUMO2 is significantly expressed in adult primary and recurrent GBMs as well as in pediatric GBM tumours, (3) Orthotropic xenografts from adult and pediatric GBMs confirm high expression of SUMO2 in GBM tumour samples, (4) SUMO2 is significantly associated with patient survival plot and pan-cancer cell fitness, (5) Inhibition of SUMO2 isoform can potentially impair multiple oncogenic kinases as well as other key protumourigenic genes across GBM subtypes, in adult and pediatric patients alike, therefore, (6) Rationale design of SUMO2 inhibitors or search for its transcriptional inhibitors is urgently required for anti-GBM and potentially pan-cancer therapeutics.
Although we had set out on our exploration to target multiple pro-proliferation associated GBM kinases via inhibition of SUMOylation, we were intrigued to find that all key upregulated coding genes, involved in the generation and maintenance of the hallmarks of aggressive cancer, were SUMOylated. These observations jives with the fact that other than kinases; transcription factors (e.g. MITF, Myc, TFAP2A, FOXM1, FOXA2, RUNXs, HIF-1α, C/EBPbeta1, Ikaros, pRB, NDRG1, PCNA, IRF-1, NF-κB, ZEB1, p53, GATA3, VHL, SIRT1, BRCA1), signalling molecules (e.g. IQGAP1, Slug, RanGAP1,Rac1, MMPs), deubiquitinating enzyme (e.g. CYLD, Pontin, HDAC3), receptor proteins (e.g. TGFbR1, androgen receptor, estrogen receptors, progesterone receptor, glucocorticoid receptor, peroxisome proliferator activated receptor, retinoic acid receptor), DDR proteins (e.g. PML, DAXX, PRC1, MDC1, HERC2, RNF168, MRN complex, the homologous recombination proteins- Mre11, Rad50, Rad52 and Rad59, and the DNA damage checkpoint proteins Rad9 and Mrc1, Rad51, ATM, CHK1, DNA-PK and KU70/KU80), metabolic enzymes (e.g. Fatty-acid synthase, HMGS-1, PKM2), cell cycle regulators (e.g. CENP-A, CENP-E, BubR1, condensins, BLM helicase, and cohesin, Cyclin-E, Topoisomerase IIa, CHK1, Ki-67 and P53) are known substrates of SUMO isoforms [10, 47,48,49,50,51,52,53,54,55].
The novelty in our findings is that a significant proportion of the upregulated genes in GBM, involved in tumourigenesis, are found to be the targets of SUMO2/3 isoforms. However, since high expression of SUMO2 alone was found to be significantly associated with poor patient survival, this SUMO isoform could putatively regulate all GBM protumourigenic processes. The cancer community had almost lost hopes that such a ‘one size fits all’ kind of tumour target therapy may be possible with an exception to the CAR-T therapy (https://www.news-medical.net/news/20200214/Universal-One-Size-Fits-All-Cancer-Treatment.aspx). This is mainly due to the challenges posed by the heterogeneity of clinical and molecular data comprising of differences in age, gender, ethnicities, co-morbidities and mutations. In the context of GBMs, the primary tumours were held to be of a different molecular landscape than recurrent ones as the latter is evidenced to arise from more cancer stem cell-like populations [56, 57]. Therefore, a notion is held that primary and recurrent tumours will require different therapeutic regimes. Besides, the pediatric tumours are genetically and epigenetically identified to be different from that adult GBMs, again conforming to the opinion that exclusive treatment protocols needs to be developed for childhood GBMs . Our analysis suggests that SUMO2 therapeutic inhibition will be beneficial to GBM patients irrespective of age, tissue differentiation and the type (primary or recurrent).
SUMOylation of proteins is clearly evidenced in cellular stress such as hypoxia, wherein a positive correlation is derived between SUMOylation and cancer growth, angiogenesis, glucose metabolism and stemness. Besides, SUMOylation of ATR and NFκB is identified as a critical component in DNA damage repair that enables resistance to chemo and radiotherapy, respectively . DDX39B, a DExD/H-box RNA helicase involved in the downregulation of factors associated with the extracellular matrix, cellular migration, and angiogenesis, is degraded upon SUMOylation which in turn promotes resistance to alkylating chemotherapy in GBM . Hence, several researches have demonstrated that pharmacological inhibition of the SUMO pathway may be a practical approach in overcoming cancer cell resistance to treatments. However, general inhibition of SUMOylation can be highly deleterious to health as it is a cytoprotective process [60, 61]. Our results indicate that mere suppression of SUMO2 isoform may be sufficient in drug sensitization vs global SUMOylation targeting.
Glioblastoma tumour-initiating cells are identified to be major culprits in GBM aggressiveness and relapse . SUMOylation of Promyelocytic leukaemia protein (PML) facilitates its interaction with c-Myc, stabilising the protumourigenic c-Myc in glioma stem cells . SUMOylation modification is reported in different stemness marker proteins like Oct-4, Oct-1, Nestin etc. [63,64,65,66]. Also, SUMOylation of retinoid acid receptors in stem cells induces resistance to RA-mediated cancer stem cells differentiation . Since SUMO2 expression is found to be high in glioblastoma tumour-initiating stem cells, SUMO2 targeting is expected to greatly benefit in controlling tumour recurrence.
SUMO modified proteins such as Transducin β-like protein (TBL1) and Transducin β-like 1X-linked receptor 1 (TBLR1, coactivator for NF-κB-mediated transcription) are associated with inflammation-mediated promotion of tumourigenesis . SUMOylation also contributes to the negative regulation of NKG2D and DNAM-1, reducing NK cell-mediated surveillance against tumours . Since inflammation-driven immune escape and fibrosis (ECM stiffening) are key processes that drive glioblastoma pathogenesis, the high expression of SUMO2 in inflammation-promoting macrophages and stromal cells indicates that it may be crucial in their maintenance [69, 70]. Therefore, SUMO2 targeting may be helpful in reducing the numbers of such aberrant immune cells, which will invariably subside inflammatory secretions and cytokine storms.
Metastasis, migration and invasion associated cancer genes such as RanGAP1 (Ran-GTPase activating enzyme 1), Rho-like GTPase- Rac1, MMPs like MMP-9, MMP-14 and slug [10, 71] are found to be profoundly SUMOylated. Although metastasis is rare outside the central nervous system in GBMs, we find that SUMO2 is crucially associated with EMT genes, putatively supporting migration and invasion processes.
One of the major reasons for the extremely poor prognosis of GBMs is its asymptomatic nature in early stages, routine exosome profiling from blood serum can enable early diagnosis and thereby early staging of GBMs or precancerous lesions. Since our analysis suggests that SUMO2 specific conjugation of certain proteins putatively occurs in GBM, pull-down assays in combination with mass spectrometry may enable the identification of such SUMOylated tumour-specific biomarkers in exosomes. Given the reports that SUMO2 and SUMO2 conjugated proteins are identified in exosomes , SUMO2 and its target proteins can serve as a circulating biomarkers for determining cancer risk and for evaluating the chemotherapy response.
Further, most of the currently available anti-cancer drugs are weakly basic in nature, and the predominant acidic microenvironment of GBMs induces protonation of these drugs, inhibiting their entry into the tumour cells and the sites of action [73, 74]. Therefore, while using HTS library screening technologies to identify inhibitors of SUMO2 or its transcriptional inhibitors, due consideration should be laid on its cellular accessibility in the low pH GBM microenvironment. Additionally, combining SUMO2 inhibitors with FDA approved drugs and other anti-cancer phyto-compounds may efficiently target GBMs through drug synergism and may combat treatment resistance.
We found that SUMO2 isoform was significantly over-expressed in both primary and recurrent GBM tumour tissues. Cellularly, both bulk and glioblastoma cancer stem cell-like population and as well as pro-tumourigenic macrophages/stromal cells over-expressed SUMO2 isoform. Molecularly, all major kinases and other genes involved in GBM formation and maintenance were found to be targets of SUMO2 conjugation. Clinically, SUMO2 isoform over-expression was significantly correlated with the shorter survival rates in the GBM patients. The over-expression was found in the GBMs of both adult and pediatric origins and also in the low-grade gliomas. In summary, SUMO2 is identified as the master regulator of the processes involved in generating cancer hallmarks; hence it is a promising anti-GBM drug target (Fig. 3I). Thus, our study calls for an urgent industry-academia collaborations to develop and validate SUMO2 isoform-specific inhibition strategies, which can target multiple hub pro-tumorigenic genes in glioblastoma pathology, irrespective of patient age and other clinicopathological parameters.
Availability of data and materials
All data is available in the manuscript and its Additional material files.
AMP-activated protein kinase
Ataxia-telangiectasia mutated serine/threonine kinase
Ataxia telangiectasia and Rad3-related protein serine/threonine kinase
Budding Uninhibited By Benzimidazoles 1 (Yeast Homolog) Beta- Mitotic Checkpoint Serine/Threonine Kinase
BReast CAncer gene 1
Chimeric antigen receptor T cells
Children's Brain Tumour Tissue Consortium
Cancer Cell Line Encyclopedia
Centromere protein A
Checkpoint kinase 1
Clinical Proteomic Tumour Analysis Consortium
CYLD lysine 63 deubiquitinase
Database for Annotation, Visualization and Integrated Discovery
Death Domain Associated Protein
DNA Damage Response
DNA-dependent protein kinase
DNAX Accessory Molecule-1
Epidermal growth factor receptor
Extracellular signal-regulated kinase 5
Food and Drug Administration
False discovery rate
Fibroblast growth factor receptors
Forkhead Box M1
Fragments Per Kilobase of transcript per Million mapped reads
Genomic Data Commons Data Portal
Gene Expression Omnibus
Gene Set Enrichment Analysis
Glycogen synthase kinase 3
Genotype Tissue Expression portal
Human Glioblastoma Cell Culture
Histone Deacetylase 3
Hypoxia Inducible Factor
HECT and RLD domain containing E3 ubiquitin protein ligase 2
3-Hydroxy-3-Methylglutaryl-Coenzyme A Synthase HNF1A: Hepatocyte nuclear factor 1-alpha
Human Protein Atlas
IQ Motif Containing GTPase Activating Protein 1
Interferon Regulatory Factor 1
Insulin like growth factor 1
Large tumour suppressor kinase 1
LIM Domain Kinase 1
Mitogen-activated protein kinase
Mixed Lineage Leukemia
Melanocyte inducing transcription factor
Maximum neighborhood component
Maximal clique centrality
Mediator of DNA damage checkpoint 1
Monopolar spindle 1
Mammalian target of rapamycin
Mre11, Rad50 and Nbn gene complex
N-Myc Downstream Regulated 1
Natural killer cells
Natural killer group 2D
Nucleophosmin-anaplastic lymphoma kinase
Serine/threonine p21-activating kinases
Pediatric Brain Tumour
Protein Data Commons
Platelet-derived growth factor receptor alpha
Pediatric tumour derived orthotropic xenografts
Proliferating cell nuclear antigen
Polycystic kidney disease
Pyruvate kinase isozymes M1
Promyelocytic leukemia protein
Protein Regulator of cytokinesis 1
Phosphatase and tensin homolog
Ras homolog gene family, member A
Ring Finger Protein 168
Rho-associated protein kinase
RNA-seq by expectation–maximization
Runt-related transcription factor
Sestrin protease 1
Silent mating type information regulation 2 homolog 1
Small ubiquitin-like modifier
The Cancer Genome Atlas
Transforming growth factor beta receptor
Transcripts per million
Transcription Factor AP-2 Alpha
- UCSC Xena:
University of California Santa Cruz, data browser
Vascular endothelial growth factor receptor
Von Hippel–Lindau syndrome
Zinc Finger E-Box Binding Homeobox 1
Zhou A. Targeting protein kinases for the treatment of glioblastoma multiforme: linking basic studies to clinical applications. Curr Pharm Des. 2017;23:4290–302. https://doi.org/10.2174/1381612823666170710144325.
Pearson JRD, Regad T. Targeting cellular pathways in glioblastoma multiforme. Signal Transduct Target Ther. 2017;2:17040. https://doi.org/10.1038/sigtrans.2017.40.
Sathornsumetee S, Reardon DA. Targeting multiple kinases in glioblastoma multiforme. Expert Opin Investig Drugs. 2009;18:277–92. https://doi.org/10.1517/13543780802692603.
Radivojac P, Baenziger PH, Kann MG, Mort ME, Hahn MW, Mooney SD. Gain and loss of phosphorylation sites in human cancer. Bioinformatics. 2008;24:i241–7. https://doi.org/10.1093/bioinformatics/btn267.
Bleeker FE, Lamba S, Zanon C, Molenaar RJ, Hulsebos TJM, Troost D, van Tilborg AA, Vandertop WP, Leenstra S, van Noorden CJF, et al. Mutational profiling of kinases in glioblastoma. BMC Cancer. 2014;14:718. https://doi.org/10.1186/1471-2407-14-718.
Zhao JJ, Roberts TM. PI3 kinases in cancer: from oncogene artifact to leading cancer target. Sci STKE. 2006. https://doi.org/10.1126/stke.3652006pe52.
Albert L, Karsy M, Murali R, Jhanwar-Uniyal M. Inhibition of mTOR activates the MAPK pathway in glioblastoma multiforme. Cancer Genomics Proteomics. 2009;6:255–61.
Peña-Ortiz MÁ, Germán-Castelán L, González-Arenas A. Growth factors and kinases in glioblastoma growth. Adv Mod Oncol Res. 2016;2(5):248–60.
Hendriks IA, D’Souza RCJ, Yang B, Verlaan-de Vries M, Mann M, Vertegaal ACO. Uncovering global SUMOylation signaling networks in a site-specific manner. Nat Struct Mol Biol. 2014;21:927–36. https://doi.org/10.1038/nsmb.2890.
Kukkula A, Ojala VK, Mendez LM, Sistonen L, Elenius K, Sundvall M. Therapeutic potential of targeting the SUMO pathway in cancer. Cancers (Basel). 2021. https://doi.org/10.3390/cancers13174402.
Fox BM, Janssen A, Estevez-Ordonez D, Gessler F, Vicario N, Chagoya G, Elsayed G, Sotoudeh H, Stetler W, Friedman GK, et al. SUMOylation in glioblastoma: a novel therapeutic target. Int J Mol Sci. 2019;20:1853. https://doi.org/10.3390/ijms20081853.
Yang W, Wang L, Roehn G, Pearlstein RD, Ali-Osman F, Pan H, Goldbrunner R, Krantz M, Harms C, Paschen W. Small ubiquitin-like modifier 1–3 conjugation [corrected] is activated in human astrocytic brain tumours and is required for glioblastoma cell survival. Cancer Sci. 2013;104:70–7. https://doi.org/10.1111/cas.12047.
Bellail AC, Olson JJ, Hao C. SUMO1 modification stabilizes CDK6 protein and drives the cell cycle and glioblastoma progression. Nat Commun. 2014;5:4234. https://doi.org/10.1038/ncomms5234.
Wang L, Zhang J, Banerjee S, Barnes L, Barnes L, Sajja V, Liu Y, Guo B, Du Y, Agarwal MK, et al. Sumoylation of vimentin354 is associated with PIAS3 inhibition of glioma cell migration. Oncotarget. 2010;1:620–7. https://doi.org/10.18632/oncotarget.101101.
Wang L, Ji S. Inhibition of Ubc9-induced CRMP2 SUMOylation disrupts glioblastoma cell proliferation. J Mol Neurosci. 2019;69:391–8. https://doi.org/10.1007/s12031-019-01368-y.
Erazo T, Espinosa-Gil S, Diéguez-Martínez N, Gómez N, Lizcano JM. SUMOylation is required for ERK5 nuclear translocation and ERK5-mediated cancer cell proliferation. Int J Mol Sci. 2020. https://doi.org/10.3390/ijms21062203.
Sehat B, Tofigh A, Lin Y, Trocmé E, Liljedahl U, Lagergren J, Larsson O. SUMOylation mediates the nuclear translocation and signaling of the IGF-1 receptor. Sci Signal. 2010;3:ra10. https://doi.org/10.1126/scisignal.2000628.
Brand TM, Iida M, Li C, Wheeler DL. The nuclear epidermal growth factor receptor signaling network and its role in cancer. Discov Med. 2011;12:419–32.
Mei L, Yuan L, Shi W, Fan S, Tang C, Fan X, Yang W, Qian Y, Hussain M, Wu X. SUMOylation of large tumour suppressor 1 at Lys751 attenuates its kinase activity and tumour-suppressor functions. Cancer Lett. 2017;386:1–11. https://doi.org/10.1016/j.canlet.2016.11.009.
He X, Li C, Ke R, Luo L, Huang D. Down-regulation of adenosine monophosphate-activated protein kinase activity: a driver of cancer. Tumour Biol J Int Soc Oncodevelopmental Biol Med. 2017;39:1010428317697576. https://doi.org/10.1177/1010428317697576.
Nakagawa S, Koricheva J, Macleod M, Viechtbauer W. Introducing our series: research synthesis and meta-research in biology. BMC Biol. 2020;18:20. https://doi.org/10.1186/s12915-020-0755-0.
Eid S, Turk S, Volkamer A, Rippmann F, Fulle S. KinMap: a web-based tool for interactive navigation through human kinome data. BMC Bioinformatics. 2017;18:16. https://doi.org/10.1186/s12859-016-1433-7.
Wilson LJ, Linley A, Hammond DE, Hood FE, Coulson JM, MacEwan DJ, Ross SJ, Slupsky JR, Smith PD, Eyers PA, et al. New perspectives, opportunities, and challenges in exploring the human protein kinome. Cancer Res. 2018;78:15–29. https://doi.org/10.1158/0008-5472.CAN-17-2291.
Goldman MJ, Craft B, Hastie M, Repečka K, McDade F, Kamath A, Banerjee A, Luo Y, Rogers D, Brooks AN, et al. Visualizing and interpreting cancer genomics data via the Xena platform. Nat Biotechnol. 2020;38:675–8. https://doi.org/10.1038/s41587-020-0546-8.
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014;15:550. https://doi.org/10.1186/s13059-014-0550-8.
Szklarczyk D, Gable AL, Lyon D, Junge A, Wyder S, Huerta-Cepas J, Simonovic M, Doncheva NT, Morris JH, Bork P, et al. STRING v11: protein–protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2018;47:D607–13. https://doi.org/10.1093/nar/gky1131.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13:2498–504. https://doi.org/10.1101/gr.1239303.
Assenov Y, Ramírez F, Schelhorn S-E, Lengauer T, Albrecht M. Computing topological parameters of biological networks. Bioinformatics. 2008;24:282–4. https://doi.org/10.1093/bioinformatics/btm554.
Chin C-H, Chen S-H, Wu H-H, Ho C-W, Ko M-T, Lin C-Y. cytoHubba: identifying hub objects and sub-networks from complex interactome. BMC Syst Biol. 2014;8:S11. https://doi.org/10.1186/1752-0509-8-S4-S11.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005;102:15545–50. https://doi.org/10.1073/pnas.0506580102.
Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44–57. https://doi.org/10.1038/nprot.2008.211.
Yu G, Wang L-G, Han Y, He Q-Y. clusterProfiler: an R package for comparing biological themes among gene clusters. Omi A J Integr Biol. 2012;16:284–7. https://doi.org/10.1089/omi.2011.0118.
Mizuno H, Kitada K, Nakai K, Sarai A. PrognoScan: a new database for meta-analysis of the prognostic value of genes. BMC Med Genomics. 2009;2:18. https://doi.org/10.1186/1755-8794-2-18.
Dempster JM, Rossen J, Kazachkova M, Pan J, Kugener G, Root DE, Tsherniak A. Extracting biological insights from the project achilles genome-scale CRISPR screens in cancer cell lines. BioRxiv. 2019. https://doi.org/10.1101/720243.
Tran AN, Dussaq AM, Kennell T, Willey CD, Hjelmeland AB. HPAanalyze: an R package that facilitates the retrieval and analysis of the Human Protein Atlas data. BMC Bioinform. 2019;20:463. https://doi.org/10.1186/s12859-019-3059-z.
Hendriks IA, Lyon D, Su D, Skotte NH, Daniel JA, Jensen LJ, Nielsen ML. Site-specific characterization of endogenous SUMOylation across species and organs. Nat Commun. 2018;9:2456. https://doi.org/10.1038/s41467-018-04957-4.
Gill BJ, Pisapia DJ, Malone HR, Goldstein H, Lei L, Sonabend A, Yun J, Samanamud J, Sims JS, Banu M, et al. MRI-localized biopsies reveal subtype-specific differences in molecular and cellular composition at the margins of glioblastoma. Proc Natl Acad Sci. 2014;111:12550. https://doi.org/10.1073/pnas.1405839111.
Klemm F, Maas RR, Bowman RL, Kornete M, Soukup K, Nassiri S, Brouland J-P, Iacobuzio-Donahue CA, Brennan C, Tabar V, et al. Interrogation of the microenvironmental landscape in brain tumours reveals disease-specific alterations of immune cells. Cell. 2020;181:1643-1660.e17. https://doi.org/10.1016/j.cell.2020.05.007.
Mistry AM, Wooten DJ, Davis LT, Mobley BC, Quaranta V, Ihrie RA. Ventricular-Subventricular zone contact by glioblastoma is not associated with molecular signatures in bulk tumour data. Sci Rep. 2019;9:1842. https://doi.org/10.1038/s41598-018-37734-w.
Murat A, Migliavacca E, Gorlia T, Lambiv WL, Shay T, Hamou M-F, de Tribolet N, Regli L, Wick W, Kouwenhoven MCM, et al. Stem cell-related “self-renewal” signature and high epidermal growth factor receptor expression associated with resistance to concomitant chemoradiotherapy in glioblastoma. J Clin Oncol. 2008;26:3015–24. https://doi.org/10.1200/JCO.2007.15.7164.
Stringer BW, Day BW, D’Souza RCJ, Jamieson PR, Ensbey KS, Bruce ZC, Lim YC, Goasdoué K, Offenhäuser C, Akgül S, et al. A reference collection of patient-derived cell line and xenograft models of proneural, classical and mesenchymal glioblastoma. Sci Rep. 2019;9:4902. https://doi.org/10.1038/s41598-019-41277-z.
Brabetz S, Leary SES, Gröbner SN, Nakamoto MW, Şeker-Cin H, Girard EJ, Cole B, Strand AD, Bloom KL, Hovestadt V, et al. A biobank of patient-derived pediatric brain tumour models. Nat Med. 2018;24:1752–61. https://doi.org/10.1038/s41591-018-0207-3.
Horita H, Law A, Hong S, Middleton K. A simple toolset to identify endogenous post-translational modifications for a target protein: a snapshot of the EGFR signaling pathway. Biosci Rep. 2017;37:1.
Lichti CF, Wildburger NC, Shavkunov AS, Mostovenko E, Liu H, Sulman EP, Nilsson CL. The proteomic landscape of glioma stem-like cells. EuPA Open Proteom. 2015;8:85–93. https://doi.org/10.1016/j.euprot.2015.06.008.
Kozuka-Hata H, Nasu-Nishimura Y, Koyama-Nasu R, Ao-Kondo H, Tsumoto K, Akiyama T, Oyama M. Global proteome analysis of glioblastoma stem cells by high-resolution mass spectrometry. Curr Top Pept Protein Res. 2012;13:1–47.
La Manno G, Siletti K, Furlan A, Gyllborg D, Vinsland E, Mossi Albiach A, Mattsson Langseth C, Khven I, Lederer AR, Dratva LM, et al. Molecular architecture of the developing mouse brain. Nature. 2021;596:92–6. https://doi.org/10.1038/s41586-021-03775-x.
Apostolov A, Litim-Mecheri I, Oravecz A, Goepp M, Kirstetter P, Marchal P, Ittel A, Mauvieux L, Chan S, Kastner P. Sumoylation inhibits the growth suppressive properties of ikaros. PLoS ONE. 2016;11: e0157767. https://doi.org/10.1371/journal.pone.0157767.
Morris JR. SUMO in the mammalian response to DNA damage. Biochem Soc Trans. 2010;38:92–7. https://doi.org/10.1042/BST0380092.
Christopoulos PF, Msaouel P, Koutsilieris M. The role of the insulin-like growth factor-1 system in breast cancer. Mol Cancer. 2015;14:43. https://doi.org/10.1186/s12943-015-0291-7.
Abrieu A, Liakopoulos D. How Does SUMO Participate in Spindle Organization? Cells. 2019. https://doi.org/10.3390/cells8080801.
Floris A, Mazarei M, Yang X, Robinson AE, Zhou J, Barberis A, et al. SUMOylation Protects FASN against proteasomal degradation in breast cancer cells treated with grape leaf extract. Biomolecules. 2020. https://doi.org/10.3390/biom10040529.
Kobayashi T, Masoumi KC, Massoumi R. Deubiquitinating activity of CYLD is impaired by SUMOylation in neuroblastoma cells. Oncogene. 2015;34:2251–60. https://doi.org/10.1038/onc.2014.159.
Treuter E, Venteclef N. Transcriptional control of metabolic and inflammatory pathways by nuclear receptor SUMOylation. Biochim Biophys Acta. 2011;1812:909–18. https://doi.org/10.1016/j.bbadis.2010.12.008.
An S, Huang L, Miao P, Shi L, Shen M, Zhao X, Liu J, Huang G. Small ubiquitin-like modifier 1 modification of pyruvate kinase M2 promotes aerobic glycolysis and cell proliferation in A549 human lung cancer cells. Onco Targets Ther. 2018;11:2097–109. https://doi.org/10.2147/OTT.S156918.
Sapir A, Tsur A, Koorman T, Ching K, Mishra P, Bardenheier A, Podolsky L, Bening-Abu-Shach U, Boxem M, Chou T-F, et al. Controlled sumoylation of the mevalonate pathway enzyme HMGS-1 regulates metabolism during aging. Proc Natl Acad Sci. 2014;111:E3880–9. https://doi.org/10.1073/pnas.1414748111.
Birzu C, French P, Caccese M, Cerretti G, Idbaih A, Zagonel V, Lombardi G. Recurrent glioblastoma: from molecular landscape to new treatment perspectives. Cancers (Basel). 2020. https://doi.org/10.3390/cancers13010047.
Zhang L, Liu Z, Li J, Huang T, Wang Y, Chang L, Zheng W, Ma Y, Chen F, Gong X, et al. Genomic analysis of primary and recurrent gliomas reveals clinical outcome related molecular features. Sci Rep. 2019;9:16058. https://doi.org/10.1038/s41598-019-52515-9.
Sturm D, Bender S, Jones DTW, Lichter P, Grill J, Becher O, Hawkins C, Majewski J, Jones C, Costello JF, et al. Paediatric and adult glioblastoma: multiform (epi)genomic culprits emerge. Nat Rev Cancer. 2014;14:92–107. https://doi.org/10.1038/nrc3655.
Szymura SJ, Bernal GM, Wu L, Zhang Z, Crawley CD, Voce DJ, Campbell P-A, Ranoa DE, Weichselbaum RR, Yamini B. DDX39B interacts with the pattern recognition receptor pathway to inhibit NF-κB and sensitize to alkylating chemotherapy. BMC Biol. 2020;18:32. https://doi.org/10.1186/s12915-020-0764-z.
Bernstock JD, Yang W, Ye DG, Shen Y, Pluchino S, Lee Y-J, Hallenbeck JM, Paschen W. SUMOylation in brain ischemia: Patterns, targets, and translational implications. J Cereb blood flow Metab Off J Int Soc Cereb Blood Flow Metab. 2018;38:5–16. https://doi.org/10.1177/0271678X17742260.
Yang W, Sheng H, Wang H. Targeting the SUMO pathway for neuroprotection in brain ischaemia. Stroke Vasc Neurol. 2016;1:101–7. https://doi.org/10.1136/svn-2016-000031.
Gimple RC, Bhargava S, Dixit D, Rich JN. Glioblastoma stem cells: lessons from the tumour hierarchy in a lethal cancer. Genes Dev. 2019;33:591–609. https://doi.org/10.1101/gad.324301.119.
Bogachek MV, Park JM, De Andrade JP, Lorenzen AW, Kulak MV, White JR, Gu VW, Wu VT, Weigel RJ. Inhibiting the SUMO pathway represses the cancer stem cell population in breast and colorectal carcinomas. Stem Cell Rep. 2016;7:1140–51. https://doi.org/10.1016/j.stemcr.2016.11.001.
Du L, Li Y-J, Fakih M, Wiatrek RL, Duldulao M, Chen Z, Chu P, Garcia-Aguilar J, Chen Y. Role of SUMO activating enzyme in cancer stem cell maintenance and self-renewal. Nat Commun. 2016;7:12326. https://doi.org/10.1038/ncomms12326.
Lee JE, Kim JH. SUMO modification regulates the protein stability of NDRG1. Biochem Biophys Res Commun. 2015;459:161–5. https://doi.org/10.1016/j.bbrc.2015.02.090.
Yao Y, Lu Y, Chen W-C, Jiang Y, Cheng T, Ma Y, Lu L, Dai W. Cobalt and nickel stabilize stem cell transcription factor OCT4 through modulating its sumoylation and ubiquitination. PLoS ONE. 2014;9: e86620. https://doi.org/10.1371/journal.pone.0086620.
Rodriguez V, Bailey R, Larion M, Gilbert MR. Retinoid receptor turnover mediated by sumoylation, ubiquitination and the valosin-containing protein is disrupted in glioblastoma. Sci Rep. 2019;9:16250. https://doi.org/10.1038/s41598-019-52696-3.
Molfetta R, Zingoni A, Santoni A, Paolini R. Post-translational mechanisms regulating NK cell activating receptors and their ligands in cancer: potential targets for therapeutic intervention. Front Immunol. 2019;10:2557. https://doi.org/10.3389/fimmu.2019.02557.
Stauffer JK, Scarzello AJ, Jiang Q, Wiltrout RH. Chronic inflammation, immune escape, and oncogenesis in the liver: a unique neighborhood for novel intersections. Hepatology. 2012;56:1567–74. https://doi.org/10.1002/hep.25674.
Hibino S, Kawazoe T, Kasahara H, Itoh S, Ishimoto T, Sakata-Yanagimoto M, Taniguchi K. Inflammation-Induced Tumourigenesis and Metastasis. Int J Mol Sci. 2021. https://doi.org/10.3390/ijms22115421.
Bogachek MV, De Andrade JP, Weigel RJ. Regulation of epithelial-mesenchymal transition through SUMOylation of transcription factors. Cancer Res. 2015;75:11–5. https://doi.org/10.1158/0008-5472.CAN-14-2824.
Wang L, Wu J, Song S, Chen H, Hu Y, Xu B, Liu J. Plasma Exosome-Derived Sentrin SUMO-specific protease 1: a prognostic biomarker in patients with osteosarcoma. Front Oncol. 2021;11: 625109. https://doi.org/10.3389/fonc.2021.625109.
Yue Z, Li C, Voth GA, Swanson JMJ. Dynamic protonation dramatically affects the membrane permeability of drug-like molecules. J Am Chem Soc. 2019;141:13421–33. https://doi.org/10.1021/jacs.9b04387.
Adams DJ, Morgan LR. Tumour physiology and charge dynamics of anticancer drugs: implications for camptothecin-based drug development. Curr Med Chem. 2011;18:1367–72. https://doi.org/10.2174/092986711795029609.
This work was supported by the Department of Biotechnology, Ministry of Science and Technology, Govt. of India for the award of following Grants: BT/RLF/Re-entry/16/2011,BT/PR4959/MED/30/770/2012,BT/PR6331/GBD/27/403/2012, BT/PR26735/MED/31/378/2017 to RM. S.J is supported by Senior Research Fellowship from the Indian Council of Medical Research (Fellowship Award ID: 2014–26530; File no. BMS/FW/SCR/2014–26530/MAR-2015/04/KL/GOVT).
Ethics approval and consent to participate
Consent for publication
All authors give consent for the publication of the manuscript in Cancer Cell International.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
GDC TCGA GBM clinical data of 173 samples.
GDC TCGA GBM mRNA expression normalized count data of 173 samples.
mRNA raw count data of 173 samples clinical samples associated with GDC TCGA GBM study.
Tables associated with data shown in main and additional figures.
Materials and methods, addiitonal figures and legends.
About this article
Cite this article
Krishna, A.P., John, S., Shinde, P.L. et al. Proteo-transcriptomics meta-analysis identifies SUMO2 as a promising target in glioblastoma multiforme therapeutics. Cancer Cell Int 21, 575 (2021). https://doi.org/10.1186/s12935-021-02279-y
- Glioblastoma multiforme
- Cancer drug target discovery
- Kinases and cancer