Abstract Objective MCM5 is a DNA licensing factor involved in cell proliferation and has been previously established as an excellent biomarker in a number of malignancies. Nevertheless, the role of MCM5 in GBM has not been fully clarified. The present study aimed to investigate the potential roles of MCM5 in the treatment of GBM and to elucidate its underlying mechanism, which is beneficial for developing new therapeutic strategies and predicting prognosis. Methods Firstly, we obtained transcriptomic and proteomic data from the TCGA and CPTAC databases on glioma patients. Employing the DeSeq2 R package, we then identified genes with joint differential expression in GBM tissues subjected to chemotherapy. To develop a prognostic risk score model, we performed univariate and multivariate Cox regression analyses. In vitro knockdown and overexpression of MCM5 were used to further investigate the biological functions of GBM cells. Additionally, we also delved into the upstream regulation of MCM5, revealing associations with several transcription factors. Finally, we investigated differences in immune cell infiltration and drug sensitivity across diverse risk groups identified in the prognostic risk model. Results In this study, the chemotherapy-treated GBM samples exhibited consistent alterations in 46 upregulated and 94 downregulated genes at both the mRNA and protein levels. Notably, MCM5 emerged as a gene with prognostic significance as well as potential therapeutic relevance. In vitro experiments subsequently validated the role of increased MCM5 expression in promoting GBM cell proliferation and resistance to TMZ. Correlations with transcription factors such as CREB1, CTCF, NFYB, NRF1, PBX1, TEAD1, and USF1 were discovered during upstream regulatory analysis, enriching our understanding of MCM5 regulatory mechanisms. The study additionally delves into immune cell infiltration and drug sensitivity, providing valuable insights for personalized treatment approaches. Conclusion This study identifies MCM5 as a key player in GBM, demonstrating its prognostic significance and potential therapeutic relevance by elucidating its role in promoting cell proliferation and resistance to chemotherapy. Keywords: glioblastoma, minichromosome maintenance protein 5, chemotherapeutic drug resistance Introduction Gliomas are categorized into four grades (Grade I, II, III, and IV) based on their malignant nature in pathologic diagnosis, with glioblastoma multiforme (GBM) identified as Grade IV gliomas. Complete surgical removal of the invasive GBM parenchyma is often unattainable due to the tumor’s rapid growth and infiltration into normal brain tissue.[38]^1–3 Consequently, chemotherapy plays a crucial role in GBM treatment. Temozolomide (TMZ) is considered a critical chemotherapeutic agent for GBM, exerting its cytotoxic effects through DNA methylation at the N^7 and O^6 positions of guanine. The drug targets nucleic acids, proteins, and peptide nucleophilic regions within tumor cells, impacting various stages of cell division.[39]^4 Its primary objective is to eliminate rapidly proliferating cells, alter the growth and division patterns of cancer cells, and potentially induce apoptosis.[40]^5^,[41]^6 Despite its widespread use, GBM carries a bleak prognosis and demonstrates notable resistance to treatment.[42]^7 Previous studies have demonstrated that drug resistance mechanisms predominantly involve non-coding RNA,[43]^8^,[44]^9 DNA damage repair mechanisms,[45]^10^,[46]^11 aberrations in signaling pathways,[47]^11 the tumor microenvironment,[48]^12^,[49]^13 the glioma stem cell theory,[50]^12 tumor cell autophagy,[51]^14 and other associated factors. Ultimately, GBM frequently develops drug resistance through these molecular mechanisms concurrently, resulting in accelerated proliferation of tumor cells during the advanced stages of treatment.[52]^4 Hence, it is imperative to comprehend the molecular mechanisms underlying the accelerated proliferation of tumors following drug resistance. This study utilized proteomic analysis from the Clinical Proteomic Tumor Analysis Consortium (CPTAC) database in conjunction with transcriptomic analysis from the Cancer Genome Atlas (TCGA) database to examine alterations in GBM expression profiles. Our results indicate that increased levels of Mini chromosome Maintenance Protein 5 (MCM5) are correlated with chemotherapy resistance in patients with GBM. The MCM5 gene encodes a protein that functions as a constituent of the MCM2-7 complex, known as the MCM complex, which serves as the replicative helicase necessary for the initiation and elongation of DNA replication in eukaryotic cells.[53]^15 It is a critical component of the CDC45-MCM-GINS helicase, a molecular apparatus responsible for unwinding template DNA during replication and serving as the foundation for the replisome.[54]^16–18 Previous studies have found that MCM5 plays important roles in the progression of various cancers, including liver, lung, bladder, and colon cancers et al[55]^19–23 Meanwhile, Jiang et al previously reported that the ubiquitination and degradation of MCM5 protein might impede the repair of DNA damage caused by cisplatin.[56]^21 However, thus far, the expression, function and potential mechanisms of MCM5 in GBM chemotherapy resistance are still unclear. The current study aimed to investigate the involvement of MCM5 and its associated mechanisms in conferring chemotherapy resistance in GBM, potentially identifying a promising therapeutic target for the treatment of this disease. Methods and Materials Data Collection The transcriptome data and clinical information for GBM were obtained from the Cancer Genome Atlas (TCGA) Data Portal ([57]https://tcga-data.nci.nih.gov/tcga/). In the TCGA cohort, transcriptome data and clinical details were collected from 594 GBM patients and 5 normal individuals. For proteomic data and clinical information related to GBM, the source was the Clinical Proteomic Tumor Analysis Consortium (CPTAC) ([58]https://proteomics.cancer.gov/programs/cptac) (S047 S048). This dataset comprised a total of 302 tissues, including 10 normal brain tissues and 292 GBM tissues. Within the GBM patient group, 144 individuals did not undergo chemotherapy, while 48 received chemotherapy. Additionally, the chemotherapy status of the remaining 100 GBM patients was not recorded, exclusively utilized for the analysis of differences compared to normal brain tissue. All data in this study were collected from public databases (TCGA and CPTAC), and approved by the medical ethics committee of Shenzhen University School of Medicine. Chemotherapy-Related Genes with Prognostic Values Were Screened Differential gene expression analyses were conducted between GBM tissues (n = 594) and normal brain tissues (n = 5) utilizing the DeSeq2 R package (v1.22.0). Subsequent differential protein expression analyses included comparisons between GBM tissues (n = 100) and normal brain tissues (n = 10), as well as between GBM patients with chemotherapy (n = 48) and those without chemotherapy (n = 144), all performed with the DeSeq2 R package (v1.22.0). Criteria for identifying differential genes included a logFC > 2 and P < 0.05. Univariate and multivariate Cox regression analyses were used to screen independent prognostic factors for GBM. The predictive accuracy of 1,3,5 year was determined by receiver operating characteristic (ROC) curve analysis using the area under the ROC curve (AUC). Cell Culture and Cell Lentiviral Transduction The human glioma cell lines, U87 and U251, were procured from Procell Life Sciences & Technology Co. Ltd (Wuhan, Hubei, China). These cells were cultured in Dulbecco modified Eagle medium (DMEM) supplemented with 10% fetal bovine serum (FBS) obtained from Gibco (CA, USA). The cell cultures were maintained at 37 °C in 5% CO2-humidified air incubators provided by Thermo Fisher Scientific (Waltham, MA, USA). To knock down MCM5 expression, lentiviral short hairpin RNAs (shRNAs) were employed. These shRNAs were constructed by inserting the respective shRNA sequences into the pHBLV-U6-MCS-CMV-ZsGreen-PGK-PURO vector, sourced from Hanbio (Shanghai, China). The specific shRNA sequence used for MCM5 knockdown is as follows: GAAACTGAAGAACCGCTACAT. For stable overexpression of MCM5, the MCM5 sequences were cloned into the pHBLV-CMV-MCS-3FLAG-EF1-ZsGreen-T2A-PURO vector. Lentiviral transfection experiments were conducted following the manufacturer’s instructions, and stable transfected cells were selected using puromycin for a period of 2 weeks. Protein Extraction and Western Blotting Analysis Total protein extraction was performed using a radioimmunoprecipitation assay (RIPA) buffer obtained from Beyotime Biotechnology (Shanghai, China). After quantification of the extracted proteins, 30 μg of total protein was separated using 10–12% SDS-PAGE gels and subsequently electrophoretically transferred onto either a 0.2μm or 0.45μm PVDF membrane sourced from Millipore Sigma (Burlington, MA, USA). The membranes were then blocked with 5% fat-free milk for 1 hour at room temperature, followed by overnight incubation at 4°C with primary antibodies acquired from Abcam (Cambridge, MA, USA): MCM5 (ab75975, 1:500) and Actin (ab115777, 1:2000). Appropriate secondary antibodies (Cell Signaling Technology, Dallas, TX, USA) were incubated with the membranes for 2 hours at room temperature. Subsequently, the protein bands were visualized using chemiluminescence and quantified using ImageJ software from Media Cybernetics (Silver Springs, MD, USA). TMZ Chemosensitivity and Flow Cytometry Analysis After stable transfection, exponentially growing GBM cells were seeded into 96-well plates overnight. Then, the medium was changed with different concentrations of TMZ (25, 50, 100, 200, 400, 800 μg/mL), and the sensitivities of GBM cells to TMZ were evaluated by IC50 values using the Cell Counting Kit-8 (CCK-8) (Dojindo, Japan). GBM cells after stable transfection were treated with 200 μg/mL TMZ for 48 h, and then apoptotic cells were assessed using the FITC Annexin-V Apoptosis Detection Kit (BD Biosciences, Franklin Lakes, NJ, USA) following the manufacturer’s recommendations. In brief, GBM cells were isolated and incubated with 300 μL of binding buffer containing 5 μL Annexin V-FITC for 10 min in darkness. Then, 5 μL of PI solution was added and incubated under dark conditions for 5 min. Cells were immediately analyzed by BD FACS Calibur (BD Biosciences, MA, USA). Based on the staining fluorescence, viable cells (FITC^−PI^−), early apoptotic cells (FITC^+ PI^−), and late apoptotic or necrotic cells (FITC^+ PI^+) were identified. Cell Viability and Cell Cloning Assays Cell viability was also assessed at different times (24, 48, 72, 96, and 120 h) using a CCK-8 kit according to manufacturer instructions. The ability of GBM cells to proliferate after different treatments was determined by a colony formation assay. In brief, GBM cells were seeded in six-well plates at a density of 300 cells/well and were treated with 25 μg/mL TMZ. After 14 days, the colonies were fixed with glutaraldehyde (6.0% v/v), stained with crystal violet (0.5% w/v), and counted using a stereomicroscope. Statistical Analysis We performed statistical analysis using SPSS 21.0 software (SPSS Inc., Chicago, IL, USA). Data are expressed as the mean ± standard deviation (SD). To assess differences between two or more groups, we employed an unpaired Student’s t-test or v2-test and a one-way analysis of variance, respectively. Tukey’s multiple comparison test was applied when necessary. Each experiment was conducted with a minimum of three biological replicates, and statistical significance was defined as p-levels < 0.05. Results Differential Chemotherapy-Related Gene Expression in GBM Patient Samples To identify genes exhibiting consistent changes at both the mRNA and protein levels in GBM compared to normal brain tissue, a comprehensive analysis was conducted. Transcriptome and clinical data were obtained from the TCGA database, while proteome and clinical data were sourced from the CPTAC database. Differential expression mRNA analysis in GBM and normal brain tissue was performed using the DeSeq2 R package based on the TCGA database, revealing 6016 mRNAs with significant changes ([59]Figure 1a, [60]Table S1). Simultaneously, differential protein expression analysis was conducted using the DeSeq2 R package based on the CPTAC database, uncovering 3812 proteins with significant changes ([61]Figure 1b, [62]Table S2). To explore the impact of chemotherapy on proteins expression, a differential gene expression analysis was performed by comparing patients receiving chemotherapy against those without chemotherapy controls based on the CPTAC database, identifying 1415 significantly differentially expressed proteins ([63]Figure 1c, [64]Table S3). Figure 1. [65]Figure 1 [66]Open in a new tab Differential chemotherapy-related gene expression in GBM patient samples. (a) Differential expression mRNA analysis of GBM and normal brain tissue (Normal) based on the TCGA database. (b) Differential expression protein analysis of GBM and normal brain tissue (Normal) based on CPTAC database. (c) Differential expression protein analysis of chemotherapy and non-chemotherapy GBM Patients based on the CPTAC database. (d) The Venn chart of the differential gene among figures (a-c). Normal (N), GBM (G), chemotherapy (C) and non-chemotherapy (NC). (e and f) GO and KEGG pathway analysis of the differentially chemotherapy-related genes. (g) Presentation of known interactions between differentially chemotherapy-related genes using String. To pinpoint the intersection of genes consistently altered at both the mRNA and protein levels between patients with and without chemotherapy and normal brain tissues, a Venn diagram was constructed. The overlapping gene set from [67]Figures 1a-c revealed 46 upregulated and 94 downregulated genes ([68]Figure 1d). Subsequent pathway enrichment analysis of these 140 annotated genes identified significant associations with cellular response to drug, DNA replication initiation, serine/threonine protein kinase complex, transporter complex, phospholipase binding, histone kinase activity, DNA replication, cell cycle, p53 signaling pathway, ErbB signaling pathway, ECM-receptor interaction, synaptic vesicle cycle, and glioma ([69]Figure 1e-f). The interaction between these differential chemotherapy-related genes was further evaluated using the STRING database ([70]https://string-db.org/) ([71]Figure 1g). This integrative approach provides a comprehensive understanding of genes consistently altered at both the mRNA and protein levels in GBM, shedding light on potential pathways and interactions implicated in the disease and its response to chemotherapy. Identification of Chemotherapy-Related Genes with Prognostic Value in the GBM In our pursuit of potential therapeutic targets with clinical relevance, we conducted screening of chemotherapy-related genes through both univariate and multivariate analyses. The forest map displays 21 prognostic chemotherapy-related genes ([72]Figure 2a). Further independent prognostic prediction analysis highlighted MCM5 as having an independent prognostic value for GBM patients ([73]Figure 2b). For 1-, 3-, and 5-year survival, the area under the curve (AUC) of ROC curves was 0.861, 0.875, 0.619, and 0.826, respectively ([74]Figure 2c). Kaplan-Meier survival analysis of the CPTAC glioma dataset underscored that elevated MCM5 expression is associated with poor patient prognosis ([75]Figure 2d). Subsequent analysis revealed a positive correlation between MCM5 protein levels and glioma grade, indicating an increase in MCM5 expression with higher glioma grades ([76]Figure 2e). Additionally, MCM5 exhibited high expression in GBM cases subjected to radiation and chemotherapy ([77]Figure 2f and [78]g). In summary, these findings strongly suggest that heightened MCM5 expression can serve as a predictive marker for poor prognosis in GBM patients, emphasizing its potential as a therapeutic target with clinical significance. Figure 2. [79]Figure 2 [80]Open in a new tab Identification of chemotherapy-related genes with prognostic value in the GBM. (a) A forest plot showing the results of a multivariate COX regression model for identifying potential risk factors for GBM. (b) Forest plots for the relationships of abnormal MCM5 protein expression with clinicopathological characteristics and the prognosis of patients with GBM. (c) ROC curve analysis according to the 1, 3, and 5-year survival of the area under the AUC value. P<0.05 (d) Kaplan-Meier survival curve analysis based on the CPTAC database suggested the prognosis of GBM patients with higher or lower MCM5 levels. (e) Expression of MCM5 protein levels in different grades of gliomas. (f) Expression of MCM5 protein level between radiation and non-radiation GBM patients. (g) Expression of MCM5 protein level between chemotherapy and non-chemotherapy GBM patients. Effect of MCM5 on GBM Cell Apoptosis, Proliferation, and Temozolomide Sensitivity In the initial phase, we conducted lentiviral knockdown or overexpression of MCM5, and the transfection efficiency was verified by Western blot analysis ([81]Figure 3a and [82]e). To assess the impact of MCM5 on glioma cell sensitivity to temozolomide and apoptosis, we carried out CCK-8 assays and flow cytometry experiments. The results indicated that MCM5 knockdown heightened temozolomide sensitivity ([83]Figure 3b) and promoted cell apoptosis ([84]Figure 3c and [85]d). Conversely, MCM5 overexpression inhibited temozolomide sensitivity ([86]Figure 3f) and reduced cell apoptosis ([87]Figure 3g and [88]h). Further insights from clone formation and CCK-8 proliferation experiments revealed that the clone formation rate and cell proliferation rate were lower in MCM5 knockdown groups compared to control groups. In contrast, MCM5 overexpression groups induced the opposite effects ([89]Figure 3i-l). These results collectively demonstrate the critical role of MCM5 in modulating proliferation, temozolomide sensitivity, and apoptosis in glioma cells. Figure 3. [90]Figure 3 [91]Open in a new tab Effect of MCM5 on GBM cell apoptosis, proliferation, and chemosensitivity. (a) Detection of the expression of MCM5 in U87 and U251 cells after MCM5 knockdown transfection (sh-MCM5), compared to control group transfection (sh-NC) by WB. (b) The IC50 detected by the CCK8 assay in U87 and U251 cells after sh-MCM5, compared to sh-NC. ** p < 0.01. (c) The ratio of total apoptosis was examined by flow cytometry in U87 and U251 cells after sh-MCM5, compared to sh-NC. (d) The statistical figure of apoptotic activity in U87 and U251 cells after sh-MCM5, compared to sh-NC. *** p < 0.001. (e) Detection of the expression of MCM5 after MCM5 overexpression transfection (lv-MCM5), compared to control group transfection (lv-NC) by WB. (f) The IC50 detected by the CCK8 assay in U87 and U251 cells after lv-MCM5, compared to lv-NC. * p < 0.05, ** p < 0.01. (g) The ratio of total apoptosis was examined by flow cytometry in U87 and U251 cells after lv-MCM5, compared to lv-NC. (h) The statistical figure of apoptotic activity in U87 and U251 cells after lv-MCM5, compared to lv-NC. *** p < 0.001. (I and j) CCK8 assays were used to evaluate U87 and U251 cell proliferation after sh-MCM5 or lv-MCM5, compared to sh-NC or lv-NC, respectively. * p < 0.05, ** p < 0.01, *** p < 0.001. (k and l) Cell clone formation assay was used to evaluate U251 cell clone formation ability after sh-MCM5 or lv-MCM5, compared to sh-NC or lv-NC, respectively. * p < 0.05. The MCM5 Upstream Regulator Analysis Identifies the Upstream Transcriptional Regulators Transcription factors play a pivotal role in the regulation of gene expression. Utilizing the UCSC Genome Browser ([92]https://genome.ucsc.edu/), we predicted sequences in the promoter region of MCM5. Subsequently, the DAVID website ([93]http://david.abcc.ncifcrf.gov/) was employed to predict potential transcription factors for MCM5. To further investigate, we analyzed the expression correlation between MCM5 and the predicted transcription factors in GBM tissue samples based on the CPTAC database. The results revealed a positive correlation between MCM5 and several transcription factors, including CREB1, CTCF, NFYB, NRF1, PBX1, TEAD1, and USF1 ([94]Figure 4a). To gain insights into the transcription factor binding motif for MCM5, we generated motifs using WebLogo ([95]http://weblogo.berkeley.edu) ([96]Figure 4b-h). These findings collectively suggest that the regulation of MCM5 expression may be influenced by additional transcription factors. The logical progression from sequence prediction to correlation analysis and motif generation strengthens our understanding of the potential regulatory network governing MCM5 expression in GBM. Figure 4. [97]Figure 4 [98]Open in a new tab The MCM5 upstream regulator analytic identifies the upstream transcriptional regulators. (a) Correlation analysis of gene expression between MCM5 and transcription factors based on the CPTAC database. (b-h) Transcription factors binding to the MCM5 promoter Motif were predicted by using the WebLogo tool (Crooks et al 2004) available at [99]http://weblogo.berkeley.edu. GBM-Infiltrating Immune Cells Were Estimated in Different Risk Groups We extended our investigation to explore the role of immune cell infiltration in tumor growth. Utilizing the previously identified 21 prognostic chemotherapy-related genes ([100]Figure 2a), we developed a prognostic risk score model, categorizing patients into high-risk and low-risk groups based on the median of the risk score. Subsequently, we assessed immune cell infiltration using XCELL, TIMER, QUANTISEQ, MCPCOUNTER, and EPIC software. Among the immune cell types analyzed, including T cell CD4^+ Th2 (XCELL), T cell CD8^+ (TIMER), T cell CD4^+ (non-regulatory), Tregs (T cell regulatory), and uncharacterized cell populations (QUANTISEQ and EPIC), positive correlations with the risk score were evident. Conversely, the remaining cell types displayed negative correlations ([101]Figure 5a). These findings underscore the potential significance of immune cells in tumor growth and provide valuable insights for further research into the complex interplay between chemotherapy-related genes and the tumor microenvironment in GBM. Figure 5. [102]Figure 5 [103]Open in a new tab GBM-infiltrating immune cells and chemosensitivity were estimated using a risk assessment model. (a) Predictive maps of GBM-infiltrating immune cells using XCELL, TIMER, QUANTISEQ, MCPCOUNTER, and EPIC software based on the TCGA database. p < 0.05. (b-e) The model acted as a potential predictor for chemosensitivity, as high-risk scores were related to a higher IC50 for chemotherapeutics. * p < 0.05, ** p < 0.01, *** p < 0.001. (f-j) The model acted as a potential predictor for chemosensitivity, as high-risk scores were related to a lower IC50 for chemotherapeutics. ** p < 0.01. Chemosensitivity Was Estimated in Different Risk Groups We further investigated differences in drug sensitivity between the high- and low-risk groups by analyzing the IC50 of chemotherapeutic agents. The results revealed significant differences in drug sensitivities between the high- and low-risk groups. Patients in the low-risk group exhibited greater sensitivity to Bexarotene, Bryostatin.1, Gefitinib, and Imatinib ([104]Figure 5b-e). Conversely, patients in the high-risk group demonstrated higher sensitivity to Camptothecin, Cytarabine, Elesclomol, Metformin, and Thapsigargin ([105]Figure 5f-j). In summary, these findings underscore the potential pivotal role of immune cells in tumors and emphasize the importance of identifying drug sensitivities in tumors with different risk profiles. This knowledge is crucial for tailoring individualized therapies in the realm of cancer treatment. Discussion The drug resistance observed in GBM is a complex phenomenon influenced by various factors, including tumor stem cells, the microenvironment, cellular stress responses to chemotherapy, drug penetration into tumor tissues, and resistance mediated by DNA damage repair pathways.[106]^7 Conventional single transcriptome or proteomics analyses have limitations in fully capturing the changes in gene expression related to GBM resistance.[107]^24^,[108]^25 To address this, our study systematically investigated chemotherapy-related gene expression differences in GBM patients by synthesizing transcriptomic and proteomic data from the TCGA and CPTAC databases. Our findings revealed consistent changes at both mRNA and protein levels, particularly in biological processes related to key pathways such as cell response to drugs, DNA replication, and the cell cycle. Notably, MCM5 emerged as a gene with independent prognostic value, showing a significant correlation between its high expression and poor prognosis in GBM patients. The positive correlation of MCM5 expression with clinical survival and tumor grade underscored its crucial role in GBM. MCM5, a member of the MCM family of chromatin-binding proteins, is implicated in cell cycle regulation and is upregulated during the transition from G0 to G1/S phase.[109]^26^,[110]^27 Given the commonality of abnormal cell cycle regulation in tumors, targeting MCM5 function holds promise as a potential strategy to inhibit the excessive proliferation of tumor cells.[111]^20^,[112]^21^,[113]^23^,[114]^28–30 Importantly, we first reported that MCM5 promotes glioma cell U87 and U251 proliferation and resistance to TMZ, aligning with previous research indicating its role in cisplatin resistance in head and neck squamous cell carcinoma. However, further studies with the incorporation of more cell lines are necessary to confirm this scenario. In recent years, patient-derived xenograft (PDX) models have garnered significant attention in research due to their capacity to maintain the microenvironment and histological features of primary tumors, including colon cancer, stomach cancer, and esophageal cancer.[115]^31–34 In contrast to cancer cell lines, PDX models accurately represent the heterogeneity of tumors at both inter- and intra-tumoral levels.[116]^34 Consequently, utilizing the PDX model is deemed more robust for investigating the involvement of MCM5 in GBM chemotherapy resistance. Through upstream regulatory analysis, we identified several transcription factors positively correlated with MCM5 expression, including CREB1,[117]^35 CTCF,[118]^36 and NFYB,[119]^37 suggesting a complex regulatory network influencing MCM5 expression. NFYB plays an important role in drug resistance, especially in oxaliplatin-resistant colon cancer. It promotes oxaliplatin resistance by upregulating the expression of E2F1, thereby inhibiting apoptosis through Sirt1-dependent deacetylation. Furthermore, the NFYB-E2F1 axis upregulated the expression of CHK1, which is a key factor in inhibiting oxaliplatin-induced apoptosis and maintaining tumor formation in OR-CRC cells. Knockdown of E2F1 reduced the expression of Pol κ, which is required for CHK1 activation, thereby making OR-CRC cells more sensitive to oxaliplatin.[120]^37 At the same time, CTCF also plays a role in regulating drug resistance, while CREB1 is closely related to sensitivity to chemotherapy drugs.[121]^38–40 These findings reveal the complexity of MCM5 expression being regulated by multiple transcription factors and provide important clues for further studying the drug resistance mechanisms of tumors and developing targeted treatment strategies. Immune checkpoint inhibitors represent a widely utilized and highly effective form of immunotherapy for combating tumors.[122]^41–43 Alongside these inhibitors, alternative immunotherapies like CAR T-cell therapy and oncolytic virus therapy have gained FDA approval.[123]^44^,[124]^45 While TMZ-induced lymphopenia has shown promise in augmenting antitumor vaccination efforts by triggering passive immune responses, it’s crucial to acknowledge the associated risks.[125]^46^,[126]^47 Lymphopenic conditions, notably in GBM patients, have been associated with compromised immune surveillance and heightened susceptibility to opportunistic infections.[127]^48 Additionally, lymphopenia can hinder active anti-tumor immune responses by depleting both T and NK cells. This depletion stems from the reduced activity of O6-methylguanine-DNA-methyltransferase (MGMT), a DNA repair enzyme, induced by TMZ.[128]^49^,[129]^50 The present study utilized a set of 21 prognostic genes associated with chemotherapy to develop a risk assessment model, categorizing patients into high- and low-risk cohorts according to their respective risk scores. Subsequently, XCELL, TIMER, QUANTISEQ, MCPCOUNTER, and EPIC software tools were employed to assess immune cell infiltration in patients from high- and low-risk groups, offering valuable insights for tailoring immune-based therapeutic strategies in the context of chemotherapy for patients with varying risk profiles. While the role of MCM5 in regulating immune cells has not been reported, studies indicate its importance in the transcription activation of Stat1 target genes, raising the possibility of its involvement in tumor immune escape and potential implications for immunotherapy responses. These implied that MCM5 is a target of immunotherapy. Finally, our analysis of drug sensitivity within the high-risk group revealed an elevated sensitivity to specific chemotherapy drugs, offering important insights for tailored therapeutic strategies. Conclusions In summary, our study sheds light on the pivotal role of MCM5 in GBM, unveiling its significance as a prognostic marker and potential therapeutic target. Through exploration, we unraveled its involvement in fostering cell proliferation and conferring resistance to chemotherapy in GBM. The revelation of intricate upstream regulatory factors, coupled with a comprehensive examination of immune cell infiltration and drug sensitivity, enriches our understanding. Altogether, these findings intricately position MCM5 as a promising and nuanced therapeutic avenue for tailoring personalized treatment strategies in the complex landscape of GBM. Acknowledgment