Senescence Reprogramming by MTHFD2 Deficiency Facilitates Tumor Progression

Senescence Reprogramming by MTHFD2 Deficiency Facilitates Tumor Progression | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Article Senescence Reprogramming by MTHFD2 Deficiency Facilitates Tumor Progression Ping Wang, Zhou Fang, Wei Pei, Qi Wu, Tingting Niu, Chengyuan Dong, and 3 more This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-4353045/v1 This work is licensed under a CC BY 4.0 License Status: Posted Version 1 posted You are reading this latest preprint version Abstract Background Age is a critical risk factor for cancer, as its incidence and mortality increase with age. However, there's limited understanding of the molecular changes aging induces in tumors. Methods We explored demographic differences between young and old cancer patients and identified age sixty and above as pivotal in cancer prognosis. Subsequently, we developed an aging-related prognostic model based on genes to assess senescence's impact on aging-associated cancer. Results Our study revealed increased genomic instability and somatic mutations in tumors from older individuals. We also found alterations in carcinogenic signaling pathways, particularly immune responses, inflammatory pathways, and cell cycle arrest in susceptible populations. Single-cell RNA sequencing showed heightened frequencies of exhausted T cells, myeloid cells, and B cells in high-risk cohorts. Conclusion MTHFD2 emerged as a crucial molecular switch regulating senescence in cancer. Its deletion promoted tumor growth by inducing cell senescence and stimulating the senescence-associated secretory phenotype (SASP) in senescent tumor cells. This highlights the need for tailored methodologies in effective cancer management. Senescence Tumor Cell cycle MTHFD2 Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Figure 6 Background Aging, a complex and inevitable biological process, involves the gradual deterioration of physiological functions and increased vulnerability to aging-related diseases. In the realm of human biology, aging encompasses an intricate interplay of genetic and environmental factors[ 1 ], leading to a myriad of changes across biological, psychological, and social dimensions. The mechanisms underlying aging are diverse, with theories ranging from programmed genetic sequences to cellular damage accumulation[ 2 ], including DNA damage, mitochondrial dysfunction, and the impact of free radicals. Telomere shortening and epigenetic changes further contribute to the complex process of aging[ 3 ]. Aging is the most significant risk factor for cancer and is associated with an exponential increase in both incidence and mortality rates among various cancer types[ 4 ]. Nevertheless, the intricate correlation between aging and the molecular determinants of cancer remains to be fully characterized. The intricate relationship between the aging process and cancer signifies a pivotal domain within gerontology and oncology[ 1 , 3 – 5 ]. The underlying mechanisms underlying aging and cancer involve a complex interplay of genetic, cellular, and environmental factors. Mechanistically, aging is associated with a range of cellular and molecular changes that contribute to increased cancer susceptibility[ 4 ]. These changes, including genomic instability, telomere shortening, epigenetic alterations, and a dysfunctional microenvironment, may contribute to carcinogenesis[ 3 , 4 , 6 ]. Consequently, we hypothesized that, owing to the variances in selective pressures stemming from changes in tissue environments with age, tumors originating from individuals at different life stages may exhibit distinct molecular landscapes. Consequently, certain molecular alterations may be more or less prevalent in older or younger patients[ 4 ]. Aging influences cellular metabolism, and modifications in metabolic pathways play a significant role in the genesis of cancer[ 4 , 7 ]. For instance, a decrease in one-carbon metabolism becomes apparent in the aging process, and the perturbation of this metabolic pathway precipitates tumor proliferation and immune subversion. One-carbon metabolism, an integral biochemical network within cellular physiology, oversees the conveyance and utilization of one-carbon entities for a myriad of cellular processes[ 8 ]. This elaborate cascade encompasses a succession of interconnected biochemical reactions indispensable for the biogenesis of nucleotides, amino acids, and other molecules that are imperative for cellular functionality[ 9 , 10 ]. Methylenetetrahydrofolate dehydrogenase 2 (MTHFD2), a pivotal enzyme in cellular biochemistry, plays a pivotal role in one-carbon reactions[ 11 ]. Positioned at the nexus of one-carbon metabolism, this enzyme catalyzes the conversion of 5,10-methylenetetrahydrofolate to 5,10-methenyltetrahydrofolate, a pivotal juncture in the synthesis of purines and thymidylate, indispensable components of DNA[ 12 ]. MTHFD2 has been the subject of exhaustive scrutiny within the domain of cancer research because of its linkage to the metabolic adaptation of tumors. Elevated expression of MTHFD2 is recurrently noted in diverse cancer manifestations due to its overexpression concomitant with heightened cellular proliferation, invasive tendencies, and resistance to chemotherapy[ 13 ]. Although the age-related implications of MTHFD2 alterations in the genesis and prognosis of cancer remain to be fully clarified, MTHFD2 has emerged as a promising target for therapeutic interventions[ 14 ]. Here, we systematically examined aging-related disparities in genomic instability, somatic copy number alterations (SCNAs), somatic mutations, pathway modifications, and gene expression across various cancer types[ 6 , 15 , 16 ]. We investigated age-correlated indicators in cancers to precisely prognosticate outcomes[ 17 , 18 ]. Moreover, MTHFD2 plays a pivotal role as a prognostic factor in cancer, and MTHFD2 loss accelerates aging-like alterations to promote tumor growth[ 19 – 21 ]. Ultimately, we elucidated the pivotal significance of MTHFD2 as a prognostic determinant in cancer, where MTHFD2 deficiency expedites senescence-like alterations to foster tumor growth. Materials and methods Datasets The clinical pathological information from 1255 breast cancer patients treated at the Breast and Thyroid Surgery Department of Wuhan University People's Hospital between 2008 and 2016. RNA-sequencing expression matrix and clinical information of breast cancer samples and para-cancerous tissues were downloaded from the Cancer Genome Atlas (TCGA) database on UCSC Xena ( https://xena.ucsc.edu/ ). Two additional independent datasets (GSE20685 and GSE58812) and single-cell RNA-seq data and bulk RNA-seq data 24 breast tumors (GSE176078) were obtained from the GEO database ( https://www.ncbi.nlm.nih.gov/geo/ ). Aging-related genes were collected from SenMayo ( https://genomics.senescence.info/cells/ ), cell age ( https://genomics.senescence.info/cells/ ), SeneQuest ( http://Senequest.net ), and the aging atlas ( https://ngdc.cncb.ac.cn/aging/index ). Somatic mutation data were downloaded from the Genomic Data Commons (GDC) ( https://portal.gdc.cancer.gov/ ). The somatic mutation data, sorted in the form of Mutation Annotation Format (MAF), were analyzed and used to calculate the Tumor Mutation Burden (TMB) using the R package maftools. Construction and validation of a aging-related prognostic signature To identify genes associated with aging and construct a prognostic signature, we conducted two types of regression analyses: univariate Cox -LASSO regression. Through this analysis, we were able to identify 19 genes: MTHFD2, EIF4EBP1, SDC1, RAD54B, LIMCH1, CAB39L, ULBP2, CACNA1H, WT1, GATA4, SYT1, ELOVL2, DOK7, S100B, BCL2A1, IFNG, FOXE1, MAP2K6. These genes were used to develop an aging-related prognostic model based on these genes. To categorize breast cancer patients, we calculated the risk score for each patient in the training set using the following formula: Risk score = ∑ni = ∑(Coefi*xi) The cancer patients were classified into high-risk and low-risk groups based on the median of their risk scores. The R package survivalROC was utilized to estimate the predictive sensitivity of the risk score. The efficacy of the model was assessed in the validation set using the same coefficient and cutoff values that were employed in the training set. Biological functional analysis between high/low-risk group patients The DESeq2 R package was utilized to analyze differentially expressed genes (DEGs). DEGs were identified using a cutoff of an adjusted p-value of less than 0.05 and a fold change of |Log2| greater than 1. Gene set enrichment analysis (GSEA) was conducted using the clusterProfiler R package. Fisher's exact test was employed to determine significant indicators, with a false discovery rate (FDR)-corrected p-value threshold of less than 0.05. Single-sample gene set enrichment analysis was performed using the GSVA R package. Gene signatures of recurrent cancer cell states were obtained from a previous study. Cell culture We generated the B16F10/MC38-MTHFD2 cell line through knockdown experiments in B16F10/MC38 cells. We targeted two sites of the MTHFD2 gene and obtained two knockdown cell lines, which we named B16F10/MC38-MTHFD2 sh1 and B16F10/MC38-MTHFD2 sh2, respectively. For our knockdown experiments, we used B16F10/MC38-NC as the negative control. All cells were cultured in Dulbecco's modified Eagle's medium (DMEM) supplemented with 10% fetal bovine serum and 1% penicillin/streptomycin at 37°C in 5% CO 2 . The shRNA sequences are listed below: MTHFD2-sh1: GCTCATGAAGAACACCATTAT MTHFD2-sh2: CGGTCATCGATGTGGGAATAA Stable cell line generation To generate stable shRNA knockdown cells, lentiviruses were generated in 293T using the PLKO.1 lentivirus packaging system with gene-specific shRNAs. After 48 hours of lentivirus production, the media was collected and cell debris was removed using a 0.45 um microfiltration membrane. The lentivirus was immediately added to cells in a 6-well dish, along with 2 µg/ml polybrene. The cells were then diluted 1/20 and transferred to a 10 cm 2 dish 24 hours after transduction. They were subsequently treated with puromycin (4 µg/ml for DLD-1 and 1 µg/ml for SH-SY5Y) until single-cell colonies could be isolated and analyzed by RT-qPCR to confirm knockdown. Cell Proliferation Assays CCK-8 (beyotime, Cell Counting Kit-8) was added to complete culture medium of the same volume, cultured for the same duration as the experimental group, and the absorbance at 450nm was measured together. Cells in logarithmic growth phase with good condition were selected to prepare cell suspension and counted. Approximately 100µl of cell suspension was seeded per well based on appropriate cell seeding density, with 4–6 replicate wells per group. The culture plate was pre-incubated in a cell culture incubator (37°C, 5% CO 2 ) for 12–24 hours to allow cells to reach the logarithmic growth phase. 10µl of CCK-8 reagent was added to each well. The culture plate was returned to the cell culture incubator and further incubated for 0.5-4 hours. The absorbance at 450nm wavelength (OD value) was measured using a microplate reader. The experiment was repeated three times, and the average of the experimental results was taken as the final experimental result. Cell viability (%) = [(As-Ab) / (Ac-Ab)] × 100%; As: Absorbance of the experimental group (including cells, culture medium, CCK-8 solution, and drug solution); Ac: Absorbance of the control group (including cells, culture medium, CCK-8 solution, without drug); Ab: Absorbance of the blank group (including culture medium, CCK-8 solution, without cells or drug). RNA isolation and Real-time PCR Total RNA was extracted from the samples using TRIzol (Vazyme, Shanghai), following the instructions provided by the manufacturer. Subsequently, cDNA was synthesized using the Reverse Transcript Kit (Vazyme). Real-time PCR was then carried out in triplicate using the SYBR Green Master Mixture (Vazyme) on the Real-time PCR Detection System (Roche). Quantification was determined based on the cycle threshold (Ct) value and calculated using the 2 -ΔΔCt method. The primer sequences are listed below: MTHFD2-forward: ACTCCCAGAGCACATTGATG MTHFD2-reverse: CCAGCCACTACCACATTCTT VEGF-forward: TCAAACCTCACCAAAGCCAG VEGF-reverse: TCTGAACAAGGCTCACAGTG P21-forward: ACATCTCAGGGCCGAAAAC P21-reverse: TGGAGACTGGGAGAGGG Preparation of cells for flow cytometry Cells were seeded according to experimental requirements, and cells were harvested when they reached the desired density. The original culture medium was collected into centrifuge tubes, 1× trypsin digestion was added, and digestion was stopped by adding the original culture medium after the specified time, followed by centrifugation. The supernatant was removed, and the cells were resuspended in 1ml pre-chilled PBS buffer and transferred to 1.5ml Eppendorf tubes, then centrifuged at 4°C, 1000g for 5 minutes. The supernatant was removed, leaving approximately 50µl, gently tapping the bottom of the tube to loosely separate the cells. The dispersed cell suspension was added to 1ml pre-chilled 70% ethanol, gently mixed by pipetting, and fixed at 4°C for at least 4 hours in the refrigerator. The fixed cells were removed, centrifuged at 4°C, 1000g for 5 minutes, the supernatant was removed, 1ml pre-chilled PBS buffer was added to resuspend the cells, followed by centrifugation, removal of the supernatant, leaving 50µl PBS buffer, and gently tapping the tube bottom to separate the cells. Dye preparation: Dyes were prepared according to the instructions based on the number of samples, with the entire process conducted in the dark. 500µl propidium iodide staining solution was added to each sample tube, mixed slowly with a pipette gun, incubated in the dark at 37°C for 30 minutes, after completion, data was saved using a flow cytometer, and subsequent processing was carried out. Western blot analysis We electrophoresed equal amounts of lysates, ranging from 30 to 50 µg, onto polyvinylidene difluoride membranes. Subsequently, the membranes were blocked using PBST with 5% milk and probed with primary antibodies, specifically Actin (1:3000; Proteintech, P62736), and P21 (1:1000, Proteintech, P63000), overnight at 4°C. After washing thrice with PBST, the membranes were incubated for 1 h at room temperature with secondary antibodies, including goat anti-rabbit IgG-HRP (1:4000, Proteintech, SA00001-2) and goat anti-mouse IgG-HRP (1:2000, Proteintech, SA00001-1). Flow cytometry The stained cells were analyzed and sorted based on DNA-A and DNA-W of the Sytox Green fluorescence signal, as well as FSC and SSC light scattering. The analysis was conducted using an LSRII flow cytometer (Becton Dickinson, San Jose, CA, U.S.A.) with an excitation wavelength of 488 nm. The cells were sorted using a FACS Digital Vantage PE flow cytometer (Becton Dickinson) with the same excitation wavelength. The selected channels are as follows: Alexa Fluor™488: Excitation at 488 nm, emission collected in the 530/30 band; PI (Propidium Iodide): Excitation at 561 nm, emission collected in the 610/20 band. Animal experimentation Seven-week-old male wild-type C57BL/6 mice were kept in a controlled environment with a 12-hour light/dark cycle, ensuring a consistent temperature and pathogen-free conditions. They had free access to food and water. The mice were sacrificed either when the tumor size reached 200 mm 2 or when clear signs of discomfort were observed, in accordance with the guidelines established by the Institutional Animal Care and Use Committee of Tongji University Cancer Center, Shanghai Tenth People's Hospital, School of Medicine, Tongji University. Statistical analysis This study independently repeated all experiments three times or more, and all data were analyzed using GraphPad Prism 8. The data are presented as mean ± SEM. Statistical differences were tested using one-way ANOVA. A significance level of P < 0.05 indicates statistical significance. In the figures, *P < 0.05, **P < 0.01, ***P < 0.001, ****P < 0.0001, and "ns" indicates no statistical significance. Results Aging-Linked Gene Signatures Forecast Cancer Prognosis To study the impact of patient age on various clinicopathological characteristics and prognoses, a meticulous analysis was conducted using the clinical data of 1255 patients. The findings suggest that surpassing the age of 60 is a consequential determinant of prognosis in cancer patients ( Table 1 ). The investigation highlights a significant difference in outcomes among different age groups, specifically emphasizing the heightened vulnerability of individuals aged 60 and above. Within the cohort of 989 patients under the age of 60, an overwhelming majority of 879 individuals showed no cancer recurrence following systematic treatment, highlighting a positive prognosis in this age group. Conversely, among the 266 patients aged over 60, there was a heightened propensity for cancer recurrence, with 57 individuals experiencing relapse despite rigorous systematic treatment. This significant disparity in recurrence rates highlights an age of over 60 years as a key factor affecting cancer prognosis, emphasizing the need for customized interventions specific to age to enhance treatment effectiveness and reduce the increased risks linked to older age[ 22 ]. Following the initial analysis, a subsequent univariate Cox regression analysis was performed on the clinical data (p < 0.05) (Fig. 1 A), revealing compelling insights into the relationships between age and cancer development and prognosis. The findings unequivocally identified age over 60 years as a significant and independent risk factor[ 21 ]. To identify aging-associated gene alterations, we conducted a thorough examination of public databases and relevant literature, compiling a collection of four aging-related gene databases: SenMayo [ 23 ], CellAge, SeneQuest [ 24 ], and the Aging Atlas[ 1 ] (Fig. 1 B). This meticulously assembled collection comprises a total of 1861 genes associated with aging, of which 1794 genes were confirmed to be expressed within our designated training set. We identified 601 aging-related genes by differentially expressed gene (DEG) analysis that distinguished between normal and malignant tissues. The resulting volcano plot vividly depicted this dichotomy, illustrating 200 upregulated genes and 401 downregulated genes within the tumor samples compared to their normal tissue counterparts. (Fig. 1 C). Univariate Cox regression and least absolute shrinkage and selection operator (LASSO) regression algorithms were used to identify genes with the most robust prognostic significance within the aging-related gene set (Fig. 1 D). Grounded in the coefficients and expression levels of these identified signature genes, a risk score was computed, enabling the classification of collected samples into aging-related high-risk and low-risk cohorts. Prognostic analysis revealed a stark divergence, with the high-risk group showing a significantly poorer prognosis (p < 0.0001), as illustrated in Fig. 1 E. The high-risk group had a significantly greater incidence of mortality than did the low-risk group. The model's predictive performance was evaluated using the area under the curve (AUC) for 1-, 3-, and 5-year overall survival (OS), which were 0.77, 0.67, and 0.65, respectively. To validate the robustness of the findings, GSE20685 and GSE58812 served as independent verification cohorts. Coherently, samples in the validation cohorts, stratified into high- and low-risk groups based on the risk scoring method derived from the training set, displayed a more adverse prognosis in the high-risk category, accompanied by a higher incidence of mortality (Fig. 1 F-G). Further exploration of prognostic status among patients at different stages ensued. The findings indicated that the senescence-related signature score could accurately predict the prognosis of patients at stages I, II, and III. Patients in the high-risk category exhibited a more unfavorable prognosis (Fig. 1 H-K). These findings provide valuable insights into the nuanced relationship between aging and cancer outcomes, prompting further exploration of senescence-tailored therapeutic approaches for more effective cancer management in older populations[ 25 ]. Aging-associated mutational patterns in cancer Numerous investigations have consistently highlighted the pivotal role of mutational patterns in tumorigenesis. In this study, a comparative analysis was conducted to examine mutational landscapes across distinct strata of patient risk. Remarkably, significant disparities were observed in the mutation status of PIK3CA, TP53, and CDH1 between the low- and high-risk cohorts[ 16 ]. The prevalence of TP53 mutations notably increased within the high-risk demographic group, hinting at its potential involvement in the progression of cancer influenced by the aging process [ 26 ] (Fig. 2 A). A meticulous examination of variant classification revealed differences between these two cohorts (Fig. 2 B). PIK3CA, TP53, and CDH1, which are recognized as somatic driver mutations and small insertions/deletions (indels), have previously been implicated in the course of cancer evolution[ 27 , 28 ]. The observed variances in frameshift deletion mutations, ranking second in the high-risk category and surpassing nonsense mutations, and conversely in the low-risk category, may be attributed to aging-associated transformations. These findings emphasize the need for thorough investigations to uncover the effects of aging on these mutational patterns and their potential contributions to the development of cancer. Such revelations hold promise for a more nuanced understanding of the intricate interplay between aging and the molecular foundations of cancer development. Tumor mutation burden (TMB), a robust quantitative metric for assessing mutation levels[ 29 ], revealed a significantly elevated TMB within the high-risk group (Fig. 2 D). Notably, a positive correlation emerged between the risk score and TMB (Fig. 2 E), accentuating their mutual interdependence[ 31 ]. Prognostic analysis further elucidated the clinical relevance of the TMB, revealing that increased TMB correlated with an unfavorable prognosis[ 32 , 33 ] (Fig. 2 F). To determine the potential synergistic impact of the TMB on prognosis, a pioneering stratified prognostic analysis amalgamating the two parameters was conducted[ 30 , 34 ]. Intriguingly, this analysis revealed an augmented prognostic predictive capacity within the training cohort. Specifically, patients with elevated TMB demonstrated a strong association with an inferior prognosis, while those with low TMB exhibited a more favorable prognosis (P = 0.0085) (Fig. 2 G). These findings underscore the significance of incorporating TMB assessment for a comprehensive understanding of prognostic dynamics in cancer. This approach offers invaluable insights for personalized therapeutic strategies based on mutational profiles. Aging-associated molecular alterations in oncogenic signaling pathways To elucidate the mechanistic link between aging and cancer, we explored the underlying pathways through Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses based on differentially expressed genes (DEGs) between patients in the high- and low-risk groups. The GO and KEGG enrichment analyses revealed a distinct enrichment profile in the high-risk group, emphasizing immune responses such as leukocyte-mediated immunity, mononuclear cell differentiation, plasma membrane signaling receptor complex, T-cell receptor complex, antigen binding, and immune receptor activity (Fig. 3 A, B). Additionally, gene set enrichment analysis (GSEA) revealed the aggregation of gene sets related to the interferon gamma/alpha response and inflammation within the high-risk group, indicating a potential association between these pathways and aging-induced cancer progression (Fig. 3 C). These findings provide valuable insights into the molecular underpinnings of the impact of aging on cancer development, emphasizing the intricate involvement of immune-related processes and signaling pathways. To comprehensively assess molecular functional disparities between patients in the high- and low-risk groups, 16 cancer-related pathway activities were calculated[ 35 , 36 ]. This study revealed distinctive patterns of cell cycle dynamics, known for their impact on tissue regeneration, function, inflammation, and tumorigenesis. Specifically, patients in the high-risk group exhibited elevated cycle signature scores, indicative of increased activity within the cell cycle [ 37 ] (Fig. 3 F). This observation suggests a potential association between heightened cell cycle activity and an elevated risk profile, offering valuable insights into the potential drivers of cancer progression influenced by aging[ 37 ]. Aging-associated alterations in the tumor microenvironment In an endeavor to meticulously investigate the distinct cellular compositions of high- and low-risk patients, we collected publicly available single-cell RNA sequencing (scRNA-seq) data paired with bulk RNA-seq data from cancers. By integrating data from 24 samples with paired bulk and scRNA-seq information, we employed the mutual nearest neighbor (MNN) algorithm to mitigate batch effects. Rigorous quality control measures were applied to each individual sample, enabling the depiction of the cancer cellular landscape at single-cell resolution using uniform manifold approximation and projection (UMAP) visualization (Fig. 4 A). Using canonical lineage markers, we accurately annotated each cell population, which included epithelial cells, cycling cells, myeloid cells, T cells, B cells, plasma cells, cancer-associated fibroblasts (CAFs), endothelial cells, and pericytes (Fig. 4 B). For example, the expression of specific markers, such as CD79A, CD79B, and MS4A1 for B cells and EPCAM, KRT8, and KRT19 for epithelial cells, facilitated precise cell type identification[ 38 , 39 ]. The analysis of cellular subpopulations revealed distinctive compositions in the high- and low-risk groups, with significantly lower T-cell levels but higher frequencies of cycling cells, myeloid cells, and B cells in the high-risk group (Fig. 4 C). Further investigation of the functions of T cells within the high- to low-risk comparison revealed upregulation of the interferon γ and α response in these T cells. In contrast, the low-risk cohort exhibited a downregulation of mTORC1 signaling and TNFα signaling via NF-κB (Fig. 4 D-G). These findings suggest that aging-driven alterations in the tumor microenvironment contribute to cancer susceptibility [ 40 ]. MTHFD2 as an aging-associated factor in cancer In the aforementioned study, MTHFD2 emerged as both a factor accentuating aging vulnerability and a protective factor against cancer progression. Consequently, MTHFD2 was chosen for a more in-depth exploration of its intricate role in aging-related cancer. To comprehensively describe the expression of MTHFD2 across a spectrum of tissues and its correlation with clinicopathological features in cancer patients, we conducted an exhaustive exploration employing an online database. MTHFD2 was markedly upregulated in 31 tumors compared with 2 tumors (Fig. 5 A). The discernible expression patterns of MTHFD2 were found to be intricately linked with modifications in ten genes, including TP53, CDH1, and MAP3KI, whereas the mutation status significantly differed between the high- and low-expression groups of MTHFD2 (Fig. 5 B). Single-cell analysis provided additional insights, showing that MTHFD2 was expressed mainly in T cells and myeloid subsets (Fig. 5 C). These simultaneous alterations in gene expression provide invaluable insights into the potential functions of MTHFD2. We conducted GO and GSEA analyses of the differentially expressed genes in the MTHFD2 high- and low-expression cohorts. These results indicated that MTHFD2 participates in nuclear division and meiotic cell cycle processes, indicating that MTHFD2 regulates cell cycle dynamics (Fig. 5 E-G). GSEA further confirmed this finding, revealing upregulated pathways, including E2F targets and G2M checkpoints (Fig. 5 H), which are integral to mechanisms regulating the cell cycle. Given the pivotal role of the cell cycle in preserving cellular homeostasis and its intricate linkage to the aging process[ 41 ], our findings suggest that MTHFD2 may contribute to aging-related modifications involved in cell cycle regulation. Subsequently, by analyzing the expression level of MTHFD2 in the young group and the aging group, it was found that the expression level of MTHFD2 was greater in the young group (Fig. 5 I), which further confirmed that MTHFD2 may be involved in aging-related modifications. MTHFD2 loss drives senescence-like alterations to foster tumor growth In our meticulous exploration of the intricate interplay between MTHFD2 and tumor growth during the aging process, we initiated MTHFD2 knockdown experiments employing B16F10 and MC38 cells as our experimental model [ 42 ] (Fig. 6 A). Strikingly, we observed notable metamorphoses in cell morphology within the MTHFD2 knockdown group, which appeared enlarged and flattened in shape, with increased cytoplasmic granularity, vacuolization, and altered nuclear morphology. Simultaneously, the levels of β-galactosidase, which serves as a biomarker for senescence, were measured in cells subjected to MTHFD2 knockdown and CDDP treatment (positive control) [ 43 , 44 ] (Fig. 6 C). Moreover, the rate of cell proliferation significantly decreased under MTHFD2 knockdown conditions (Fig. 6 B). Additionally, the perturbation of the cell cycle in B16F10 and MC38 cells was meticulously examined through flow cytometry, which revealed a pronounced blockade in the G0/G1 phase within the MTHFD2 knockdown group (Fig. 6 D). This finding strongly suggested a potential link between MTHFD2 manipulation and impediments in the progression of the cell cycle. Concurrently, the mRNA and protein levels of P21, a pivotal checkpoint in the cell cycle, were increased in the MTHFD2 knockdown group (Fig. 6 E). Consistent results confirmed the presence of senescence-associated secretory factors such as VEGF in the MTHFD2 knockdown group (Fig. 6 F). We investigated this possibility by employing B16F10 cells with MTHFD2 knockdown in a subcutaneous tumor model. In contrast, MTHFD2 knockdown promoted faster tumor growth (Fig. 6 G, H). Immunohistochemical analysis revealed significantly greater P21, VEGF, and IL-8 staining in tumor tissues from the MTHFD2 knockdown group than in those from the control group. Additionally, there was a notable increase in CD8 + T-cell staining (Fig. 6 I). This multifaceted approach not only underscores the potential involvement of MTHFD2 loss in propelling a senescent state but also highlights its potential significance in modulating tumor behavior. Discussion Aging, an intricate biological process, is intricately intertwined with increased susceptibility to senescence-related diseases, notably cancer. Our detailed examination of senescence-related disparities across diverse cancer types revealed age to be a pivotal determinant of cancer prognosis, with individuals aged 60 years and older exhibiting increased vulnerability. By summarizing four different aging-related databases and using the Cox-LASSO algorithm, we obtained aging-associated gene alterations and aging-related risk scores, coupled with the integration of the TMB and exploration of oncogenic signaling pathways, which enriches our understanding of the molecular intricacies shaping cancer outcomes. The dissection of the cellular landscape within different risk groups reveals the complexities of the tumor microenvironment, providing a nuanced perspective. We identified MTHFD2 as both an aging risk factor and a protective factor in aging-related cancer. By analyzing the cell composition of the high- and low-risk groups, it was found that the T-cell composition of the high-risk group was significantly lower than that of the low-risk group, and T-cell exhaustion plays a pivotal role in rendering older individuals vulnerable to infections and cancer[ 45 ]. This decline is responsible for elevated susceptibility to infection and cancer[ 46 ]. Further investigations into the functions of T cells within the high- to low-risk group revealed the upregulation of the interferon γ and α response in these T cells, and the roles of IFN-γ and IFN-α are complex. While they contribute to immune surveillance against cancer, persistent activation or dysregulation of these interferons may also play a role in chronic inflammation, which is a hallmark of aging and a risk factor for cancer[ 47 , 48 ]. Methylenetetrahydrofolate dehydrogenase 2 (MTHFD2) is a crucial enzyme involved in cellular metabolism, particularly in the folate metabolic pathway. This enzyme plays a profound role in nucleotide synthesis, which is essential for DNA replication and cellular proliferation. MTHFD2 catalyzes the conversion of 5,10-methylenetetrahydrofolate to 5,10-methenyltetrahydrofolate, an essential step in the synthesis of purines, amino acids, and other important molecules within the cell. MTHFD2 has garnered attention for its involvement in various cellular processes linked to cancer progression. Multiple investigations have underscored the propensity for increased MTHFD2 expression across diverse cancer types, including breast, lung, colorectal, and pancreatic malignancies[ 19 ]. Notably, MTHFD2 is preferentially upregulated in undifferentiated or poorly differentiated tumors[ 49 ], suggesting its cancer-specific expression pattern. In addition to its metabolic functions, MTHFD2 has been implicated in promoting cancer immune evasion. By steering the folate cycle toward the maintenance of adequate UDP-GlcNAc levels, MTHFD2 propels the O-GlcNAcylation of cMYC, thereby bolstering cMYC stability and PD-L1 transcription.[ 50 ] Nevertheless, the functional repertoire of MTHFD2 is intricate and multifaceted. Although MTHFD2 is primarily expressed within the mitochondria, it is also present in the nucleus. TH9619 effectively inhibits both the dehydrogenase and cyclohydrolase activities of MTHFD1/2. MTHFD1, found in the cytoplasm, shares similar functions with MTHFD2. TH9619 inhibits MTHFD1 activity, preventing the incorporation of formate produced by the mitochondria into dTMP. It selectively targets nuclear MTHFD2 without affecting mitochondrial MTHFD2, leading to an overflow of formate from the mitochondria. Consequently, cancer cells undergo cell death despite high MTHFD2 expression[ 51 ]. In addition, recent studies have indicated that MTHFD2 inhibition induces apoptosis solely in UQCR11-null cells, while it does not affect UQCR11-intact cells[ 52 ]. This suggests that the apoptotic effect of MTHFD2 inhibition is specific to cancer cells. Inversely, under conditions of oxygen deprivation or impaired electron transport chain activity, MTHFD2 maintains a substantial supply of NADH through the promotion of serine catabolism, leading to cell death. In hypoxic cells with compromised respiration, inhibition of MTHFD2 partially restores NADH levels and facilitates cell proliferation[ 49 ]. These findings suggest a potential anticancer role for MTHFD2. Our findings reveal that deletion of MTHFD2 contributes to cellular senescence. Similarly, prior research has revealed that the colocalization of MTHFD2 with DNA replication sites in the nucleus promotes cell cycle progression. Deletion of MTHFD2 leads to S-phase cell cycle arrest and fosters a senescence-like state[ 53 ]. Additionally, MTHFD2 is expressed in developing embryos but is notably absent in most healthy adult tissues, including proliferating ones[ 54 ]. This observation implies a decrease in MTHFD2 expression with age. MTHFD2 deletion increases the infiltration of CD8 + T cells in B16F10 tumors. The potential mechanisms are manifold. Initially, senescent cancer cells manifest hyperploidy, rendering them genomically unstable and enabling the presentation of tumor antigens to activate immunosurveillance, which entails the recruitment of immune effectors such as B, NK, NKT, and T cells[ 55 , 56 ]. Furthermore, the senescence-associated secretory phenotype (SASP) factor VEGF fosters angiogenesis, thereby facilitating T cell infiltration. Nonetheless, the buildup of senescent cells in tumor-bearing mice precipitates T-cell exhaustion[ 57 ]. For instance, IL-8, elevated in MTHFD2-depleted cells, augments PD-1 expression in CD8 + T cells, leading to T-cell exhaustion. Therefore, the deletion of MTHFD2 facilitates the infiltration of CD8 + T cells into tumors but compromises their functional activity. Hence, we propose a novel perspective suggesting that MTHFD2 exerts its anticancer effect by retarding the aging process. However, some limitations are noted in this study. First, numerous studies have demonstrated that aging is a significant contributor to the development and progression of cancer. Based on data analysis related to breast cancer and animal experiments, this study revealed that downregulation of MTHFD2 can promote cellular senescence, thereby accelerating tumor growth. As such, our study may lack sufficient novelty. Second, the mechanisms involved in this study were not further investigated. For instance, we did not explore how downregulation of MTHFD2 promotes cellular senescence and its impact on tumor immunity during this process. The prognostic role of MTHFD2 in cancer patients requires validation in additional cohorts to confirm its prognostic significance. Finally, the exact effect of MTHFD2 downregulation on the tumor microenvironment remains unclear due to the lack of single-cell sequencing or multi-parameter flow cytometry analyses. Conclusions In summary, the process of aging leads to heightened genetic mutations within tumors, triggers the activation of pathways essential for tumor proliferation and spread, and promotes an immunosuppressive microenvironment. The identification of Mthfd2 as a suppressor in tumors associated with aging implies that maintaining its expression and functionality may offer a potential therapeutic approach for addressing such tumors. Declarations Ethics approval and consent to participate The animal experiment of the study was approved by the Institutional Animal Care and Use Committee of Tongji University Cancer Center, Shanghai Tenth People's Hospital, School of Medicine, Tongji University. (Issue No. 22KN151). Consent for publication Not applicable. Availability of data and materials The original contributions presented in the study are included in the article/additional file. Competing interests The author (s) declare no potential conficts of interest with respect to the research, authorship, and/or publication of this article. Funding This project was supported by the National Natural Science Foundation of China (82203629 and 82372780), Shanghai Pujiang Program (22PJD054 and 22PJ1412000) and The Fundamental Research Funds for the Central University (22120240320). Author contribution: PW designed the study. PW, ZF, and WP performed the experiments. PW and ZF analyzed the data. PW drafted the manuscript. QW revised the manuscript. All authors read and approved the final manuscript. Acknowledgements We are grateful to all the editors and reviewers for their hard work, and all the authors for their efforts. References López-Otín C, Blasco MA, Partridge L, Serrano M, Kroemer G. The Hallmarks of Aging. Cell. 2013;153:1194–217. Niccoli T, Partridge L. Ageing as a Risk Factor for Disease. Curr Biol. 2012;22:R741–52. Vijg J, Suh Y. Genetics of Longevity and Aging. Annu Rev Med. 2005;56:193–212. Campisi J. Aging, Cellular Senescence, and Cancer. Annu Rev Physiol. 2013;75:685–705. Hurria A, Jones L, Muss HB. Cancer Treatment as an Accelerated Aging Process: Assessment, Biomarkers, and Interventions. 2016; Hanahan D, Weinberg RA. Hallmarks of Cancer: The Next Generation. Cell. 2011;144:646–74. Franceschi C, Garagnani P, Parini P, Giuliani C, Santoro A. Inflammaging: a new immune–metabolic viewpoint for age-related diseases. Nat Rev Endocrinol. 2018;14:576–90. Lionaki E, Ploumi C, Tavernarakis N. One-Carbon Metabolism: Pulling the Strings behind Aging and Neurodegeneration. Cells. 2022;11:214. Tedeschi PM, Markert EK, Gounder M, Lin H, Dvorzhinski D, Dolfi SC, et al. Contribution of serine, folate and glycine metabolism to the ATP, NADPH and purine requirements of cancer cells. Cell Death Dis. 2013;4:e877–e877. Yang M, Vousden KH. Serine and one-carbon metabolism in cancer. Nat Rev Cancer. 2016;16:650–62. Ducker GS, Chen L, Morscher RJ, Ghergurovich JM, Esposito M, Teng X, et al. Reversal of Cytosolic One-Carbon Flux Compensates for Loss of the Mitochondrial Folate Pathway. Cell Metab. 2016;23:1140–53. Nilsson R, Jain M, Madhusudhan N, Sheppard NG, Strittmatter L, Kampf C, et al. Metabolic enzyme expression highlights a key role for MTHFD2 and the mitochondrial folate pathway in cancer. Nat Commun. 2014;5:3128. Fan J, Ye J, Kamphorst JJ, Shlomi T, Thompson CB, Rabinowitz JD. Quantitative flux analysis reveals folate-dependent NADPH production. Nature. 2014;510:298–302. Green NH, Galvan DL, Badal SS, Chang BH, LeBleu VS, Long J, et al. MTHFD2 links RNA methylation to metabolic reprogramming in renal cell carcinoma. Oncogene. 2019;38:6211–25. Galili T. dendextend: an R package for visualizing, adjusting and comparing trees of hierarchical clustering. Bioinformatics. 2015;31:3718–20. Lawrence MS, Stojanov P, Polak P, Kryukov GV, Cibulskis K, Sivachenko A, et al. Mutational heterogeneity in cancer and the search for new cancer-associated genes. Nature. 2013;499:214–8. Zhou M, Zhao H, Wang Z, Cheng L, Yang L, Shi H, et al. Identification and validation of potential prognostic lncRNA biomarkers for predicting survival in patients with multiple myeloma. J Exp Clin Cancer Res. 2015;34:102. Tang Z, Li C, Kang B, Gao G, Li C, Zhang Z. GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 2017;45:W98–102. Nilsson R, Jain M, Madhusudhan N, Sheppard NG, Strittmatter L, Kampf C, et al. Metabolic enzyme expression highlights a key role for MTHFD2 and the mitochondrial folate pathway in cancer. Nat Commun. 2014;5:3128. Yang Y, Lane AN, Ricketts CJ, Sourbier C, Wei M-H, Shuch B, et al. Metabolic Reprogramming for Producing Energy and Reducing Power in Fumarate Hydratase Null Cells from Hereditary Leiomyomatosis Renal Cell Carcinoma. Tan M, editor. PLoS ONE. 2013;8:e72179. Zhang B, Zheng A, Hydbring P, Ambroise G, Ouchida AT, Goiny M, et al. PHGDH Defines a Metabolic Subtype in Lung Adenocarcinomas with Poor Prognosis. Cell Rep. 2017;19:2289–303. Balducci L, Ershler WB. Cancer and ageing: a nexus at several levels. Nat Rev Cancer. 2005;5:655–62. Saul D, Kosinsky RL, Atkinson EJ, Doolittle ML, Zhang X, LeBrasseur NK, et al. A new gene set identifies senescent cells and predicts senescence-associated pathways across tissues. Nat Commun. 2022;13:4827. Gorgoulis V, Adams PD, Alimonti A, Bennett DC, Bischof O, Bishop C, et al. Cellular Senescence: Defining a Path Forward. Cell. 2019;179:813–27. Hurria A, Togawa K, Mohile SG, Owusu C, Klepin HD, Gross CP, et al. Predicting Chemotherapy Toxicity in Older Adults With Cancer: A Prospective Multicenter Study. J Clin Oncol. 2011;29:3457–65. Muller PAJ, Vousden KH. p53 mutations in cancer. Nat Cell Biol. 2013;15:2–8. Futreal PA, Coin L, Marshall M, Down T, Hubbard T, Wooster R, et al. A census of human cancer genes. Nat Rev Cancer. 2004;4:177–83. Kandoth C, McLellan MD, Vandin F, Ye K, Niu B, Lu C, et al. Mutational landscape and significance across 12 major cancer types. Nature. 2013;502:333–9. Sha D, Jin Z, Budzcies J, Kluck K, Stenzinger A, Sinicrope FA. Tumor Mutational Burden (TMB) as a Predictive Biomarker in Solid Tumors. Cancer Discov. 2020;10:1808–25. Samstein RM, Lee C-H, Shoushtari AN, Hellmann MD, Shen R, Janjigian YY, et al. Tumor mutational load predicts survival after immunotherapy across multiple cancer types. Nat Genet. 2019;51:202–6. Chalmers ZR, Connelly CF, Fabrizio D, Gay L, Ali SM, Ennis R, et al. Analysis of 100,000 human cancer genomes reveals the landscape of tumor mutational burden. Genome Med. 2017;9:34. Chan TA, Yarchoan M, Jaffee E, Swanton C, Quezada SA, Stenzinger A, et al. Development of tumor mutation burden as an immunotherapy biomarker: utility for the oncology clinic. Ann Oncol. 2019;30:44–56. Hellmann MD, Nathanson T, Rizvi H, Creelan BC, Sanchez-Vega F, Ahuja A, et al. Genomic Features of Response to Combination Immunotherapy in Patients with Advanced Non-Small-Cell Lung Cancer. Cancer Cell. 2018;33:843-852.e4. Hellmann MD, Ciuleanu T-E, Pluzanski A, Lee JS, Otterson GA, Audigier-Valette C, et al. Nivolumab plus Ipilimumab in Lung Cancer with a High Tumor Mutational Burden. N Engl J Med. 2018;378:2093–104. Hanahan D, Weinberg RA. Hallmarks of Cancer: The Next Generation. Cell. 2011;144:646–74. Vogelstein B, Papadopoulos N, Velculescu VE, Zhou S, Diaz LA, Kinzler KW. Cancer Genome Landscapes. Science. 2013;339:1546–58. Malumbres M, Barbacid M. Cell cycle, CDKs and cancer: a changing paradigm. Nat Rev Cancer. 2009;9:153–66. Palanichamy A, Apeltsin L, Kuo TC, Sirota M, Wang S, Pitts SJ, et al. Immunoglobulin class-switched B cells form an active immune axis between CNS and periphery in multiple sclerosis. Sci Transl Med [Internet]. 2014 [cited 2024 Jan 11];6. Available from: https://www.science.org/doi/10.1126/scitranslmed.3008930 Falini B, Fizzotti M, Pucciarini A, Bigerna B, Marafioti T, Gambacorta M, et al. A monoclonal antibody (MUM1p) detects expression of the MUM1/IRF4 protein in a subset of germinal center B cells, plasma cells, and activated T cells. Blood. 2000;95:2084–92. De Visser KE, Joyce JA. The evolving tumor microenvironment: From cancer initiation to metastatic outgrowth. Cancer Cell. 2023;41:374–403. Hartl FU. Cellular Homeostasis and Aging. Annu Rev Biochem. 2016;85:1–4. Boccaletto P, Machnicka MA, Purta E, Piątkowski P, Bagiński B, Wirecki TK, et al. MODOMICS: a database of RNA modification pathways. 2017 update. Nucleic Acids Res. 2018;46:D303–7. Lee BY, Han JA, Im JS, Morrone A, Johung K, Goodwin EC, et al. Senescence‐associated β‐galactosidase is lysosomal β‐galactosidase. Aging Cell. 2006;5:187–95. Debacq-Chainiaux F, Erusalimsky JD, Campisi J, Toussaint O. Protocols to detect senescence-associated beta-galactosidase (SA-βgal) activity, a biomarker of senescent cells in culture and in vivo. Nat Protoc. 2009;4:1798–806. Goronzy JJ, Weyand CM. Mechanisms underlying T cell ageing. Nat Rev Immunol. 2019;19:573–83. Mittelbrunn M, Kroemer G. Hallmarks of T cell aging. Nat Immunol. 2021;22:687–98. Takaoka A, Mitani Y, Sato M, Yokochi T, Noguchi S, Tanaka N, et al. Cross Talk Between Interferon-y and - d P Signaling Components in Caveolar Membrane Domains. 2000;288. Coelho MA, Cooper S, Strauss ME, Karakoc E, Bhosle S, Gonçalves E, et al. Base editing screens map mutations affecting interferon-γ signaling in cancer. Cancer Cell. 2023;41:288-303.e6. Yang L, Canaveras JCG, Chen Z, Wang L, Liang L, Jang C, et al. Serine Catabolism Feeds NADH when Respiration Is Impaired. Cell Metab. 2020;31:809-821.e6. Shang M, Yang H, Yang R, Chen T, Fu Y, Li Y, et al. The folate cycle enzyme MTHFD2 induces cancer immune evasion through PD-L1 up-regulation. Nat Commun. 2021;12:1940. Green AC, Marttila P, Kiweler N, Chalkiadaki C, Wiita E, Cookson V, et al. Formate overflow drives toxic folate trapping in MTHFD1 inhibited cancer cells. Nat Metab. 2023;5:642–59. Achreja A, Yu T, Mittal A, Choppara S, Animasahun O, Nenwani M, et al. Metabolic collateral lethal target identification reveals MTHFD2 paralogue dependency in ovarian cancer. Nat Metab. 2022;4:1119–37. Bonagas N, Gustafsson NMS, Henriksson M, Marttila P, Gustafsson R, Wiita E, et al. Pharmacological targeting of MTHFD2 suppresses acute myeloid leukemia by inducing thymidine depletion and replication stress. Nat Cancer. 2022;3:156–72. Nilsson R, Jain M, Madhusudhan N, Sheppard NG, Strittmatter L, Kampf C, et al. Metabolic enzyme expression highlights a key role for MTHFD2 and the mitochondrial folate pathway in cancer. Nat Commun. 2014;5:3128. Senovilla L, Vitale I, Martins I, Tailler M, Pailleret C, Michaud M, et al. An Immunosurveillance Mechanism Controls Cancer Cell Ploidy. Science. 2012;337:1678–84. Aranda F, Chaba K, Bloy N, Garcia P, Bordenave C, Martins I, et al. Immune effectors responsible for the elimination of hyperploid cancer cells. OncoImmunology. 2018;7:e1463947. Maggiorani D, Le O, Lisi V, Landais S, Moquin-Beaudry G, Lavallée VP, et al. Senescence drives immunotherapy resistance by inducing an immunosuppressive tumor microenvironment. Nat Commun. 2024;15:2435. Tables Table1：Clinical information of breast cancer patients Variable No recurrence Recurrence Age at diagnosis, years ≤60 879 110 >60 209 57 Histopathological grade 1 55 1 2 307 56 3 165 11 T stage T1+T2 693 72 T3+T4 48 7 N stage N0 477 40 N1+N2+N3 57 39 Table2：Abbreviations SASP Senescence-associated secretory phenotype MTHFD2 Methylenetetrahydrofolate dehydrogenase 2 SCNAs Somatic copy number alterations MAF Mutation Annotation Format GDC Genomic Data Commons TMB Tumor Mutation Burden DEGs Differentially expressed genes FDR False discovery rate GO Gene Ontology KEGG Kyoto Encyclopedia of Genes and Genomes GSEA Gene set enrichment analysis scRNA-seq Single-cell RNA sequencing MNN Mutual nearest neighbor UMAP Uniform manifold approximation and projection CAFs Cancer-associated fibroblasts Additional Declarations No competing interests reported. Cite Share Download PDF Status: Posted Version 1 posted You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-4353045","acceptedTermsAndConditions":true,"allowDirectSubmit":true,"archivedVersions":[],"articleType":"Research Article","associatedPublications":[],"authors":[{"id":300282396,"identity":"c7b5cc36-5a39-4bb7-bb10-5ab543756f45","order_by":0,"name":"Ping Wang","email":"","orcid":"","institution":"Anhui University of Science and Technology","correspondingAuthor":false,"prefix":"","firstName":"Ping","middleName":"","lastName":"Wang","suffix":""},{"id":300282397,"identity":"0af950ad-c74d-4520-b79a-83f134712c9f","order_by":1,"name":"Zhou Fang","email":"","orcid":"","institution":"Department of Breast and Thyroid Surgery, Renmin Hospital of Wuhan University, Wuhan, Hubei, China","correspondingAuthor":false,"prefix":"","firstName":"Zhou","middleName":"","lastName":"Fang","suffix":""},{"id":300282398,"identity":"67813db9-a7a8-470e-ac42-2ee55e51c92b","order_by":2,"name":"Wei Pei","email":"","orcid":"","institution":"Tongji University Cancer Center, Tongji University","correspondingAuthor":false,"prefix":"","firstName":"Wei","middleName":"","lastName":"Pei","suffix":""},{"id":300282399,"identity":"8577cea5-6949-456b-802f-de2479cc83dc","order_by":3,"name":"Qi Wu","email":"","orcid":"","institution":"Tongji University Cancer Center, Shanghai Tenth People’s Hospital, School of Medicine, Tongji University, Shanghai, China","correspondingAuthor":false,"prefix":"","firstName":"Qi","middleName":"","lastName":"Wu","suffix":""},{"id":300282400,"identity":"f6f9e6d3-f415-4c2a-bfe3-3f15bb0299be","order_by":4,"name":"Tingting Niu","email":"","orcid":"","institution":"Tongji University Cancer Center, Shanghai Tenth People’s Hospital, School of Medicine, Tongji University, Shanghai, China","correspondingAuthor":false,"prefix":"","firstName":"Tingting","middleName":"","lastName":"Niu","suffix":""},{"id":300282401,"identity":"205645bf-7e1b-4249-8831-e69f122912bb","order_by":5,"name":"Chengyuan Dong","email":"","orcid":"","institution":"Anhui University of Science and Technology","correspondingAuthor":false,"prefix":"","firstName":"Chengyuan","middleName":"","lastName":"Dong","suffix":""},{"id":300282402,"identity":"430db467-dc91-45a4-b925-144b3d86bc81","order_by":6,"name":"Mingkang Wu","email":"","orcid":"","institution":"Tongji University Cancer Center, Tongji University","correspondingAuthor":false,"prefix":"","firstName":"Mingkang","middleName":"","lastName":"Wu","suffix":""},{"id":300282403,"identity":"06f11416-44cf-4d01-aa6f-0b566021a6dd","order_by":7,"name":"Bei Li","email":"","orcid":"","institution":"Department of Pathology, Renmin Hospital of Wuhan University, Wuhan, Hubei, P. R. China","correspondingAuthor":false,"prefix":"","firstName":"Bei","middleName":"","lastName":"Li","suffix":""},{"id":300282404,"identity":"7c8dac48-3f80-4913-b7b8-622af6b01c31","order_by":8,"name":"Zhijie Gao","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAA1UlEQVRIiWNgGAWjYBAC/gYG9g8VBgxyDBIgLhsRWiQOMDAYSBgwGBOvxQCEgaoTG4jXwn74QIJEgU36hts9Bgwfyg4z8M9uIKCFJy3hgIRBWu6GO2cMGGecO8wgcecAfi2GB3IMGyQMDuduu5FjwMzbdhjozAQCtpx/A/S5wf90M5CWv0RpuZFjxihhcCABrIWRGC0SN56lMW4wSDbcfyOt4GDPuXQeiRsEtPD3Jx//8OCPnbzkjOSND36UWcvxzyCgBQUcAGIeEtSPglEwCkbBKMAFAKTiQ5uFYNl4AAAAAElFTkSuQmCC","orcid":"","institution":"Department of Breast and Thyroid Surgery, Renmin Hospital of Wuhan University, Wuhan, Hubei, China","correspondingAuthor":true,"prefix":"","firstName":"Zhijie","middleName":"","lastName":"Gao","suffix":""}],"badges":[],"createdAt":"2024-05-01 07:52:05","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-4353045/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-4353045/v1","draftVersion":[],"editorialEvents":[],"editorialNote":"","failedWorkflow":false,"files":[{"id":56477050,"identity":"1af3ddcd-55e4-4d36-b79f-e563c70f7181","added_by":"auto","created_at":"2024-05-14 17:46:08","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":1217926,"visible":true,"origin":"","legend":"\u003cp\u003eAging-Linked Gene Signatures Forecast Cancer Prognosis\u003c/p\u003e\n\u003cp\u003eA. Univariate Cox regression analysis based on cross-validation and minimum partial likelihood deviation of clinical information. B. Venn diagram of four aging-related gene databases. C. Volcano plot construction using fold-change values and p-adjustment. Red dots represent upregulated genes; blue dots represent downregulated genes; gray dots represent non-significant genes. D. Cox-LASSO regression algorithm identifies genes with the most robust prognostic significance in the aging-related gene set. E-F. Kaplan-Meier analysis of overall survival (OS) curves for high/low-risk subgroups of patients in the training cohort. G-J. Kaplan-Meier analysis of overall survival (OS) curves for high/low-risk groups of patients in the GEO training cohort. K-N. Survival analysis of high/low-risk groups at different stages.\u003c/p\u003e","description":"","filename":"Fig1.png","url":"https://assets-eu.researchsquare.com/files/rs-4353045/v1/938cf986c49556fc938428ef.png"},{"id":56477055,"identity":"fa49e7e9-5770-4827-95f8-ae6d643d3346","added_by":"auto","created_at":"2024-05-14 17:46:08","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":1150809,"visible":true,"origin":"","legend":"\u003cp\u003eAging-associated mutational patterns in cancer\u003c/p\u003e\n\u003cp\u003eA. Comparison of mutation profiles between high and low-risk groups. B-C. Summary plots of the cohort displaying the distribution of mutated strains based on mutation type and SNV classification, with stacked bar graphs showing the top 10 mutated genes.\u003c/p\u003e\n\u003cp\u003eD. Boxplot illustrating the correlation between risk scores and TMB in the breast cancer cohort. E. Correlation between risk scores and TMB levels. F. Kaplan-Meier survival analysis of TMB, risk scores, and OS in the breast cancer cohort. G. Kaplan-Meier survival analysis of TMB and OS in the TCGA breast cancer cohort.\u003c/p\u003e","description":"","filename":"Fig2.png","url":"https://assets-eu.researchsquare.com/files/rs-4353045/v1/a2eb681b16ca3dacad40474e.png"},{"id":56477053,"identity":"5f13d0dd-a840-40bc-824d-055bd2616138","added_by":"auto","created_at":"2024-05-14 17:46:08","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":977075,"visible":true,"origin":"","legend":"\u003cp\u003eAging-associated molecular alterations in oncogenic signalling pathways\u003c/p\u003e\n\u003cp\u003eA. GO enrichment analysis of upregulated aging-related DEGs. B. KEGG pathways of upregulated aging-related DEGs. C-E. Gene set enrichment analysis (GSEA) of aging-related prognostic models. F. Boxplot of label scores for 16 cancer cell states in high/low-risk groups based on GSVA scores. Paired two-sided Wilcoxon test.\u003c/p\u003e","description":"","filename":"Fig3.png","url":"https://assets-eu.researchsquare.com/files/rs-4353045/v1/cf859755fd67fbaee7f531fc.png"},{"id":56477051,"identity":"a19f2ab0-02f8-4307-a06f-c5bdd5b294e3","added_by":"auto","created_at":"2024-05-14 17:46:08","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":712873,"visible":true,"origin":"","legend":"\u003cp\u003eAging-associated alterations in tumor microenvironment\u003c/p\u003e\n\u003cp\u003eA. UMAP plot depicting major cellular subgroups in cancer. B. Bubble heatmap displaying expression levels of specific feature genes in cancer. The size of the circles represents the proportion of expressing cells, colored based on standardized expression levels. C. Relative proportions of different cell types in high/low-risk tumors.\u003c/p\u003e\n\u003cp\u003eD-G. GSEA analysis of T cells in high/low-risk groups.\u003c/p\u003e","description":"","filename":"Fig4.png","url":"https://assets-eu.researchsquare.com/files/rs-4353045/v1/070a4eba29eee6e71f73e3c8.png"},{"id":56477052,"identity":"c42fe76a-dab7-45bd-8379-3c744bab4f13","added_by":"auto","created_at":"2024-05-14 17:46:08","extension":"png","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":1388349,"visible":true,"origin":"","legend":"\u003cp\u003eMTHFD2 as an aging-associated factor in cancer\u003c/p\u003e\n\u003cp\u003eA. Differential expression of MTHFD2 in different cancer types based on pan-cancer analysis. B. Oncoplot displaying the somatic landscape of the tumor cohort. Genes are sorted by their mutation frequency, and samples are sorted by MTHFD2 expression, as indicated by the color bar (bottom). The side bar graph shows the -log10-transformed q values estimated using MutSigCV. The waterfall plot shows the mutation information for each gene in each sample. Color annotations for various cancer types are displayed at the bottom. The barplot above the legend shows the number of mutations per sample.\u003c/p\u003e\n\u003cp\u003eC. Boxplot showing the expression levels of MTHFD2 in tumor tissues compared to normal cells. D. Volcano plot constructed using fold-change values and p-adjusted values. Red dots represent upregulated genes; blue dots represent downregulated genes; gray dots represent non-significant genes. E-G. GO and GSEA analysis of differentially expressed genes in MTHFD2 high-expressing and low-expressing cohorts. H. Expression of MTHFD2 in different age groups. I. Expression of MTHFD2 in different staging groups.\u003c/p\u003e","description":"","filename":"Fig5.png","url":"https://assets-eu.researchsquare.com/files/rs-4353045/v1/798895a227a01220a5c9c146.png"},{"id":56477056,"identity":"1ebfcf89-83c8-46e2-b02e-7d9e11b152ad","added_by":"auto","created_at":"2024-05-14 17:46:08","extension":"png","order_by":6,"title":"Figure 6","display":"","copyAsset":false,"role":"figure","size":4480623,"visible":true,"origin":"","legend":"\u003cp\u003eMTHFD2 loss drives aging-like alterations to foster tumor growth\u003c/p\u003e\n\u003cp\u003eA. RT-qPCR analysis of MTHFD2 knockdown. B. Growth rate of MTHFD2 knockdown cells. C. β-galactosidase staining level in MTHFD2 knockdown group. D. Flow cytometry analysis of cell cycle in MTHFD2 knockdown cells. E. Protein level of P21 in MTHFD2 knockdown group. F. mRNA levels of P21 and VEGF in MTHFD2 knockdown group. G. Tumor growth curve of C57BL/6J mice injected with MTHFD2 knockdown B16F10 cells. H. Photographs of dissected tumors (n ≥ 6). I. Immunohistochemical staining in tumor tissues with MTHFD2 knockdown.\u003c/p\u003e","description":"","filename":"Fig6.png","url":"https://assets-eu.researchsquare.com/files/rs-4353045/v1/1098601652e4c4bd297b4f35.png"},{"id":56624586,"identity":"44c3ccbb-e7e4-4147-8151-ab61a38dc6e6","added_by":"auto","created_at":"2024-05-16 20:01:39","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":9295556,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-4353045/v1/fb166733-bbb0-4fb6-bb4a-c9a0e8ba806f.pdf"}],"financialInterests":"No competing interests reported.","formattedTitle":"Senescence Reprogramming by MTHFD2 Deficiency Facilitates Tumor Progression","fulltext":[{"header":"Background","content":"\u003cp\u003eAging, a complex and inevitable biological process, involves the gradual deterioration of physiological functions and increased vulnerability to aging-related diseases. In the realm of human biology, aging encompasses an intricate interplay of genetic and environmental factors[\u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e], leading to a myriad of changes across biological, psychological, and social dimensions. The mechanisms underlying aging are diverse, with theories ranging from programmed genetic sequences to cellular damage accumulation[\u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2\u003c/span\u003e], including DNA damage, mitochondrial dysfunction, and the impact of free radicals. Telomere shortening and epigenetic changes further contribute to the complex process of aging[\u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e3\u003c/span\u003e]. Aging is the most significant risk factor for cancer and is associated with an exponential increase in both incidence and mortality rates among various cancer types[\u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e4\u003c/span\u003e]. Nevertheless, the intricate correlation between aging and the molecular determinants of cancer remains to be fully characterized.\u003c/p\u003e \u003cp\u003eThe intricate relationship between the aging process and cancer signifies a pivotal domain within gerontology and oncology[\u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e, \u003cspan additionalcitationids=\"CR4\" citationid=\"CR3\" class=\"CitationRef\"\u003e3\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e5\u003c/span\u003e]. The underlying mechanisms underlying aging and cancer involve a complex interplay of genetic, cellular, and environmental factors. Mechanistically, aging is associated with a range of cellular and molecular changes that contribute to increased cancer susceptibility[\u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e4\u003c/span\u003e]. These changes, including genomic instability, telomere shortening, epigenetic alterations, and a dysfunctional microenvironment, may contribute to carcinogenesis[\u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e3\u003c/span\u003e, \u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e4\u003c/span\u003e, \u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e6\u003c/span\u003e]. Consequently, we hypothesized that, owing to the variances in selective pressures stemming from changes in tissue environments with age, tumors originating from individuals at different life stages may exhibit distinct molecular landscapes. Consequently, certain molecular alterations may be more or less prevalent in older or younger patients[\u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e4\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eAging influences cellular metabolism, and modifications in metabolic pathways play a significant role in the genesis of cancer[\u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e4\u003c/span\u003e, \u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e7\u003c/span\u003e]. For instance, a decrease in one-carbon metabolism becomes apparent in the aging process, and the perturbation of this metabolic pathway precipitates tumor proliferation and immune subversion. One-carbon metabolism, an integral biochemical network within cellular physiology, oversees the conveyance and utilization of one-carbon entities for a myriad of cellular processes[\u003cspan citationid=\"CR8\" class=\"CitationRef\"\u003e8\u003c/span\u003e]. This elaborate cascade encompasses a succession of interconnected biochemical reactions indispensable for the biogenesis of nucleotides, amino acids, and other molecules that are imperative for cellular functionality[\u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e9\u003c/span\u003e, \u003cspan citationid=\"CR10\" class=\"CitationRef\"\u003e10\u003c/span\u003e]. Methylenetetrahydrofolate dehydrogenase 2 (MTHFD2), a pivotal enzyme in cellular biochemistry, plays a pivotal role in one-carbon reactions[\u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e11\u003c/span\u003e]. Positioned at the nexus of one-carbon metabolism, this enzyme catalyzes the conversion of 5,10-methylenetetrahydrofolate to 5,10-methenyltetrahydrofolate, a pivotal juncture in the synthesis of purines and thymidylate, indispensable components of DNA[\u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e12\u003c/span\u003e]. MTHFD2 has been the subject of exhaustive scrutiny within the domain of cancer research because of its linkage to the metabolic adaptation of tumors. Elevated expression of MTHFD2 is recurrently noted in diverse cancer manifestations due to its overexpression concomitant with heightened cellular proliferation, invasive tendencies, and resistance to chemotherapy[\u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e13\u003c/span\u003e]. Although the age-related implications of MTHFD2 alterations in the genesis and prognosis of cancer remain to be fully clarified, MTHFD2 has emerged as a promising target for therapeutic interventions[\u003cspan citationid=\"CR14\" class=\"CitationRef\"\u003e14\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eHere, we systematically examined aging-related disparities in genomic instability, somatic copy number alterations (SCNAs), somatic mutations, pathway modifications, and gene expression across various cancer types[\u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e6\u003c/span\u003e, \u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e15\u003c/span\u003e, \u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e16\u003c/span\u003e]. We investigated age-correlated indicators in cancers to precisely prognosticate outcomes[\u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e17\u003c/span\u003e, \u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e18\u003c/span\u003e]. Moreover, MTHFD2 plays a pivotal role as a prognostic factor in cancer, and MTHFD2 loss accelerates aging-like alterations to promote tumor growth[\u003cspan additionalcitationids=\"CR20\" citationid=\"CR19\" class=\"CitationRef\"\u003e19\u003c/span\u003e\u0026ndash;\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e]. Ultimately, we elucidated the pivotal significance of MTHFD2 as a prognostic determinant in cancer, where MTHFD2 deficiency expedites senescence-like alterations to foster tumor growth.\u003c/p\u003e"},{"header":"Materials and methods","content":"\u003cdiv id=\"Sec3\" class=\"Section2\"\u003e \u003ch2\u003eDatasets\u003c/h2\u003e \u003cp\u003eThe clinical pathological information from 1255 breast cancer patients treated at the Breast and Thyroid Surgery Department of Wuhan University People's Hospital between 2008 and 2016. RNA-sequencing expression matrix and clinical information of breast cancer samples and para-cancerous tissues were downloaded from the Cancer Genome Atlas (TCGA) database on UCSC Xena (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://xena.ucsc.edu/\u003c/span\u003e\u003cspan address=\"https://xena.ucsc.edu/\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e). Two additional independent datasets (GSE20685 and GSE58812) and single-cell RNA-seq data and bulk RNA-seq data 24 breast tumors (GSE176078) were obtained from the GEO database (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://www.ncbi.nlm.nih.gov/geo/\u003c/span\u003e\u003cspan address=\"https://www.ncbi.nlm.nih.gov/geo/\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e). Aging-related genes were collected from SenMayo (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://genomics.senescence.info/cells/\u003c/span\u003e\u003cspan address=\"https://genomics.senescence.info/cells/\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e), cell age (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://genomics.senescence.info/cells/\u003c/span\u003e\u003cspan address=\"https://genomics.senescence.info/cells/\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e), SeneQuest (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttp://Senequest.net\u003c/span\u003e\u003cspan address=\"http://Senequest.net\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e), and the aging atlas (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://ngdc.cncb.ac.cn/aging/index\u003c/span\u003e\u003cspan address=\"https://ngdc.cncb.ac.cn/aging/index\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e). Somatic mutation data were downloaded from the Genomic Data Commons (GDC) (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://portal.gdc.cancer.gov/\u003c/span\u003e\u003cspan address=\"https://portal.gdc.cancer.gov/\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e). The somatic mutation data, sorted in the form of Mutation Annotation Format (MAF), were analyzed and used to calculate the Tumor Mutation Burden (TMB) using the R package maftools.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec4\" class=\"Section2\"\u003e \u003ch2\u003eConstruction and validation of a aging-related prognostic signature\u003c/h2\u003e \u003cp\u003eTo identify genes associated with aging and construct a prognostic signature, we conducted two types of regression analyses: univariate Cox -LASSO regression. Through this analysis, we were able to identify 19 genes: MTHFD2, EIF4EBP1, SDC1, RAD54B, LIMCH1, CAB39L, ULBP2, CACNA1H, WT1, GATA4, SYT1, ELOVL2, DOK7, S100B, BCL2A1, IFNG, FOXE1, MAP2K6. These genes were used to develop an aging-related prognostic model based on these genes. To categorize breast cancer patients, we calculated the risk score for each patient in the training set using the following formula:\u003c/p\u003e \u003cp\u003eRisk score = \u0026sum;ni = \u0026sum;(Coefi*xi)\u003c/p\u003e \u003cp\u003eThe cancer patients were classified into high-risk and low-risk groups based on the median of their risk scores. The R package survivalROC was utilized to estimate the predictive sensitivity of the risk score. The efficacy of the model was assessed in the validation set using the same coefficient and cutoff values that were employed in the training set.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec5\" class=\"Section2\"\u003e \u003ch2\u003eBiological functional analysis between high/low-risk group patients\u003c/h2\u003e \u003cp\u003eThe DESeq2 R package was utilized to analyze differentially expressed genes (DEGs). DEGs were identified using a cutoff of an adjusted p-value of less than 0.05 and a fold change of |Log2| greater than 1. Gene set enrichment analysis (GSEA) was conducted using the clusterProfiler R package. Fisher's exact test was employed to determine significant indicators, with a false discovery rate (FDR)-corrected p-value threshold of less than 0.05. Single-sample gene set enrichment analysis was performed using the GSVA R package. Gene signatures of recurrent cancer cell states were obtained from a previous study.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec6\" class=\"Section2\"\u003e \u003ch2\u003eCell culture\u003c/h2\u003e \u003cp\u003eWe generated the B16F10/MC38-MTHFD2 cell line through knockdown experiments in B16F10/MC38 cells. We targeted two sites of the MTHFD2 gene and obtained two knockdown cell lines, which we named B16F10/MC38-MTHFD2 sh1 and B16F10/MC38-MTHFD2 sh2, respectively. For our knockdown experiments, we used B16F10/MC38-NC as the negative control. All cells were cultured in Dulbecco's modified Eagle's medium (DMEM) supplemented with 10% fetal bovine serum and 1% penicillin/streptomycin at 37\u0026deg;C in 5% CO\u003csup\u003e2\u003c/sup\u003e. The shRNA sequences are listed below:\u003c/p\u003e \u003cp\u003eMTHFD2-sh1: GCTCATGAAGAACACCATTAT\u003c/p\u003e \u003cp\u003eMTHFD2-sh2: CGGTCATCGATGTGGGAATAA\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec7\" class=\"Section2\"\u003e \u003ch2\u003eStable cell line generation\u003c/h2\u003e \u003cp\u003eTo generate stable shRNA knockdown cells, lentiviruses were generated in 293T using the PLKO.1 lentivirus packaging system with gene-specific shRNAs. After 48 hours of lentivirus production, the media was collected and cell debris was removed using a 0.45 um microfiltration membrane. The lentivirus was immediately added to cells in a 6-well dish, along with 2 \u0026micro;g/ml polybrene. The cells were then diluted 1/20 and transferred to a 10 cm\u003csup\u003e2\u003c/sup\u003e dish 24 hours after transduction. They were subsequently treated with puromycin (4 \u0026micro;g/ml for DLD-1 and 1 \u0026micro;g/ml for SH-SY5Y) until single-cell colonies could be isolated and analyzed by RT-qPCR to confirm knockdown.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec8\" class=\"Section2\"\u003e \u003ch2\u003eCell Proliferation Assays\u003c/h2\u003e \u003cp\u003eCCK-8 (beyotime, Cell Counting Kit-8) was added to complete culture medium of the same volume, cultured for the same duration as the experimental group, and the absorbance at 450nm was measured together. Cells in logarithmic growth phase with good condition were selected to prepare cell suspension and counted. Approximately 100\u0026micro;l of cell suspension was seeded per well based on appropriate cell seeding density, with 4\u0026ndash;6 replicate wells per group. The culture plate was pre-incubated in a cell culture incubator (37\u0026deg;C, 5% CO\u003csup\u003e2\u003c/sup\u003e) for 12\u0026ndash;24 hours to allow cells to reach the logarithmic growth phase. 10\u0026micro;l of CCK-8 reagent was added to each well. The culture plate was returned to the cell culture incubator and further incubated for 0.5-4 hours. The absorbance at 450nm wavelength (OD value) was measured using a microplate reader. The experiment was repeated three times, and the average of the experimental results was taken as the final experimental result. Cell viability (%) = [(As-Ab) / (Ac-Ab)] \u0026times; 100%; As: Absorbance of the experimental group (including cells, culture medium, CCK-8 solution, and drug solution); Ac: Absorbance of the control group (including cells, culture medium, CCK-8 solution, without drug); Ab: Absorbance of the blank group (including culture medium, CCK-8 solution, without cells or drug).\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec9\" class=\"Section2\"\u003e \u003ch2\u003eRNA isolation and Real-time PCR\u003c/h2\u003e \u003cp\u003eTotal RNA was extracted from the samples using TRIzol (Vazyme, Shanghai), following the instructions provided by the manufacturer. Subsequently, cDNA was synthesized using the Reverse Transcript Kit (Vazyme). Real-time PCR was then carried out in triplicate using the SYBR Green Master Mixture (Vazyme) on the Real-time PCR Detection System (Roche). Quantification was determined based on the cycle threshold (Ct) value and calculated using the 2\u003csup\u003e-ΔΔCt\u003c/sup\u003e method. The primer sequences are listed below:\u003c/p\u003e \u003cp\u003eMTHFD2-forward: ACTCCCAGAGCACATTGATG\u003c/p\u003e \u003cp\u003eMTHFD2-reverse: CCAGCCACTACCACATTCTT\u003c/p\u003e \u003cp\u003eVEGF-forward: TCAAACCTCACCAAAGCCAG\u003c/p\u003e \u003cp\u003eVEGF-reverse: TCTGAACAAGGCTCACAGTG\u003c/p\u003e \u003cp\u003eP21-forward: ACATCTCAGGGCCGAAAAC\u003c/p\u003e \u003cp\u003eP21-reverse: TGGAGACTGGGAGAGGG\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec10\" class=\"Section2\"\u003e \u003ch2\u003ePreparation of cells for flow cytometry\u003c/h2\u003e \u003cp\u003eCells were seeded according to experimental requirements, and cells were harvested when they reached the desired density. The original culture medium was collected into centrifuge tubes, 1\u0026times; trypsin digestion was added, and digestion was stopped by adding the original culture medium after the specified time, followed by centrifugation. The supernatant was removed, and the cells were resuspended in 1ml pre-chilled PBS buffer and transferred to 1.5ml Eppendorf tubes, then centrifuged at 4\u0026deg;C, 1000g for 5 minutes. The supernatant was removed, leaving approximately 50\u0026micro;l, gently tapping the bottom of the tube to loosely separate the cells. The dispersed cell suspension was added to 1ml pre-chilled 70% ethanol, gently mixed by pipetting, and fixed at 4\u0026deg;C for at least 4 hours in the refrigerator. The fixed cells were removed, centrifuged at 4\u0026deg;C, 1000g for 5 minutes, the supernatant was removed, 1ml pre-chilled PBS buffer was added to resuspend the cells, followed by centrifugation, removal of the supernatant, leaving 50\u0026micro;l PBS buffer, and gently tapping the tube bottom to separate the cells. Dye preparation: Dyes were prepared according to the instructions based on the number of samples, with the entire process conducted in the dark. 500\u0026micro;l propidium iodide staining solution was added to each sample tube, mixed slowly with a pipette gun, incubated in the dark at 37\u0026deg;C for 30 minutes, after completion, data was saved using a flow cytometer, and subsequent processing was carried out.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec11\" class=\"Section2\"\u003e \u003ch2\u003eWestern blot analysis\u003c/h2\u003e \u003cp\u003eWe electrophoresed equal amounts of lysates, ranging from 30 to 50 \u0026micro;g, onto polyvinylidene difluoride membranes. Subsequently, the membranes were blocked using PBST with 5% milk and probed with primary antibodies, specifically Actin (1:3000; Proteintech, P62736), and P21 (1:1000, Proteintech, P63000), overnight at 4\u0026deg;C. After washing thrice with PBST, the membranes were incubated for 1 h at room temperature with secondary antibodies, including goat anti-rabbit IgG-HRP (1:4000, Proteintech, SA00001-2) and goat anti-mouse IgG-HRP (1:2000, Proteintech, SA00001-1).\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec12\" class=\"Section2\"\u003e \u003ch2\u003eFlow cytometry\u003c/h2\u003e \u003cp\u003eThe stained cells were analyzed and sorted based on DNA-A and DNA-W of the Sytox Green fluorescence signal, as well as FSC and SSC light scattering. The analysis was conducted using an LSRII flow cytometer (Becton Dickinson, San Jose, CA, U.S.A.) with an excitation wavelength of 488 nm. The cells were sorted using a FACS Digital Vantage PE flow cytometer (Becton Dickinson) with the same excitation wavelength. The selected channels are as follows: Alexa Fluor\u0026trade;488: Excitation at 488 nm, emission collected in the 530/30 band; PI (Propidium Iodide): Excitation at 561 nm, emission collected in the 610/20 band.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec13\" class=\"Section2\"\u003e \u003ch2\u003eAnimal experimentation\u003c/h2\u003e \u003cp\u003eSeven-week-old male wild-type C57BL/6 mice were kept in a controlled environment with a 12-hour light/dark cycle, ensuring a consistent temperature and pathogen-free conditions. They had free access to food and water. The mice were sacrificed either when the tumor size reached 200 mm\u003csup\u003e2\u003c/sup\u003e or when clear signs of discomfort were observed, in accordance with the guidelines established by the Institutional Animal Care and Use Committee of Tongji University Cancer Center, Shanghai Tenth People's Hospital, School of Medicine, Tongji University.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec14\" class=\"Section2\"\u003e \u003ch2\u003eStatistical analysis\u003c/h2\u003e \u003cp\u003eThis study independently repeated all experiments three times or more, and all data were analyzed using GraphPad Prism 8. The data are presented as mean\u0026thinsp;\u0026plusmn;\u0026thinsp;SEM. Statistical differences were tested using one-way ANOVA. A significance level of P\u0026thinsp;\u0026lt;\u0026thinsp;0.05 indicates statistical significance. In the figures, *P\u0026thinsp;\u0026lt;\u0026thinsp;0.05, **P\u0026thinsp;\u0026lt;\u0026thinsp;0.01, ***P\u0026thinsp;\u0026lt;\u0026thinsp;0.001, ****P\u0026thinsp;\u0026lt;\u0026thinsp;0.0001, and \"ns\" indicates no statistical significance.\u003c/p\u003e \u003c/div\u003e"},{"header":"Results","content":"\u003cdiv id=\"Sec16\" class=\"Section2\"\u003e \u003ch2\u003eAging-Linked Gene Signatures Forecast Cancer Prognosis\u003c/h2\u003e \u003cp\u003eTo study the impact of patient age on various clinicopathological characteristics and prognoses, a meticulous analysis was conducted using the clinical data of 1255 patients. The findings suggest that surpassing the age of 60 is a consequential determinant of prognosis in cancer patients (\u003cb\u003eTable\u0026nbsp;1\u003c/b\u003e). The investigation highlights a significant difference in outcomes among different age groups, specifically emphasizing the heightened vulnerability of individuals aged 60 and above. Within the cohort of 989 patients under the age of 60, an overwhelming majority of 879 individuals showed no cancer recurrence following systematic treatment, highlighting a positive prognosis in this age group. Conversely, among the 266 patients aged over 60, there was a heightened propensity for cancer recurrence, with 57 individuals experiencing relapse despite rigorous systematic treatment. This significant disparity in recurrence rates highlights an age of over 60 years as a key factor affecting cancer prognosis, emphasizing the need for customized interventions specific to age to enhance treatment effectiveness and reduce the increased risks linked to older age[\u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e22\u003c/span\u003e]. Following the initial analysis, a subsequent univariate Cox regression analysis was performed on the clinical data (p\u0026thinsp;\u0026lt;\u0026thinsp;0.05) (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eA), revealing compelling insights into the relationships between age and cancer development and prognosis. The findings unequivocally identified age over 60 years as a significant and independent risk factor[\u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e21\u003c/span\u003e].\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eTo identify aging-associated gene alterations, we conducted a thorough examination of public databases and relevant literature, compiling a collection of four aging-related gene databases: SenMayo [\u003cspan citationid=\"CR23\" class=\"CitationRef\"\u003e23\u003c/span\u003e], CellAge, SeneQuest [\u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e24\u003c/span\u003e], and the Aging Atlas[\u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e1\u003c/span\u003e] (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eB). This meticulously assembled collection comprises a total of 1861 genes associated with aging, of which 1794 genes were confirmed to be expressed within our designated training set. We identified 601 aging-related genes by differentially expressed gene (DEG) analysis that distinguished between normal and malignant tissues. The resulting volcano plot vividly depicted this dichotomy, illustrating 200 upregulated genes and 401 downregulated genes within the tumor samples compared to their normal tissue counterparts. (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eC). Univariate Cox regression and least absolute shrinkage and selection operator (LASSO) regression algorithms were used to identify genes with the most robust prognostic significance within the aging-related gene set (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eD). Grounded in the coefficients and expression levels of these identified signature genes, a risk score was computed, enabling the classification of collected samples into aging-related high-risk and low-risk cohorts. Prognostic analysis revealed a stark divergence, with the high-risk group showing a significantly poorer prognosis (p\u0026thinsp;\u0026lt;\u0026thinsp;0.0001), as illustrated in Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eE. The high-risk group had a significantly greater incidence of mortality than did the low-risk group. The model's predictive performance was evaluated using the area under the curve (AUC) for 1-, 3-, and 5-year overall survival (OS), which were 0.77, 0.67, and 0.65, respectively. To validate the robustness of the findings, GSE20685 and GSE58812 served as independent verification cohorts. Coherently, samples in the validation cohorts, stratified into high- and low-risk groups based on the risk scoring method derived from the training set, displayed a more adverse prognosis in the high-risk category, accompanied by a higher incidence of mortality (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eF-G). Further exploration of prognostic status among patients at different stages ensued. The findings indicated that the senescence-related signature score could accurately predict the prognosis of patients at stages I, II, and III. Patients in the high-risk category exhibited a more unfavorable prognosis (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003eH-K). These findings provide valuable insights into the nuanced relationship between aging and cancer outcomes, prompting further exploration of senescence-tailored therapeutic approaches for more effective cancer management in older populations[\u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e25\u003c/span\u003e].\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec17\" class=\"Section2\"\u003e \u003ch2\u003eAging-associated mutational patterns in cancer\u003c/h2\u003e \u003cp\u003eNumerous investigations have consistently highlighted the pivotal role of mutational patterns in tumorigenesis. In this study, a comparative analysis was conducted to examine mutational landscapes across distinct strata of patient risk. Remarkably, significant disparities were observed in the mutation status of PIK3CA, TP53, and CDH1 between the low- and high-risk cohorts[\u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e16\u003c/span\u003e]. The prevalence of TP53 mutations notably increased within the high-risk demographic group, hinting at its potential involvement in the progression of cancer influenced by the aging process [\u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e26\u003c/span\u003e] (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eA). A meticulous examination of variant classification revealed differences between these two cohorts (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eB). PIK3CA, TP53, and CDH1, which are recognized as somatic driver mutations and small insertions/deletions (indels), have previously been implicated in the course of cancer evolution[\u003cspan citationid=\"CR27\" class=\"CitationRef\"\u003e27\u003c/span\u003e, \u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e28\u003c/span\u003e]. The observed variances in frameshift deletion mutations, ranking second in the high-risk category and surpassing nonsense mutations, and conversely in the low-risk category, may be attributed to aging-associated transformations. These findings emphasize the need for thorough investigations to uncover the effects of aging on these mutational patterns and their potential contributions to the development of cancer. Such revelations hold promise for a more nuanced understanding of the intricate interplay between aging and the molecular foundations of cancer development.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eTumor mutation burden (TMB), a robust quantitative metric for assessing mutation levels[\u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e29\u003c/span\u003e], revealed a significantly elevated TMB within the high-risk group (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eD). Notably, a positive correlation emerged between the risk score and TMB (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eE), accentuating their mutual interdependence[\u003cspan citationid=\"CR31\" class=\"CitationRef\"\u003e31\u003c/span\u003e]. Prognostic analysis further elucidated the clinical relevance of the TMB, revealing that increased TMB correlated with an unfavorable prognosis[\u003cspan citationid=\"CR32\" class=\"CitationRef\"\u003e32\u003c/span\u003e, \u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e33\u003c/span\u003e] (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eF). To determine the potential synergistic impact of the TMB on prognosis, a pioneering stratified prognostic analysis amalgamating the two parameters was conducted[\u003cspan citationid=\"CR30\" class=\"CitationRef\"\u003e30\u003c/span\u003e, \u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e34\u003c/span\u003e]. Intriguingly, this analysis revealed an augmented prognostic predictive capacity within the training cohort. Specifically, patients with elevated TMB demonstrated a strong association with an inferior prognosis, while those with low TMB exhibited a more favorable prognosis (P\u0026thinsp;=\u0026thinsp;0.0085) (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eG). These findings underscore the significance of incorporating TMB assessment for a comprehensive understanding of prognostic dynamics in cancer. This approach offers invaluable insights for personalized therapeutic strategies based on mutational profiles.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec18\" class=\"Section2\"\u003e \u003ch2\u003eAging-associated molecular alterations in oncogenic signaling pathways\u003c/h2\u003e \u003cp\u003eTo elucidate the mechanistic link between aging and cancer, we explored the underlying pathways through Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses based on differentially expressed genes (DEGs) between patients in the high- and low-risk groups. The GO and KEGG enrichment analyses revealed a distinct enrichment profile in the high-risk group, emphasizing immune responses such as leukocyte-mediated immunity, mononuclear cell differentiation, plasma membrane signaling receptor complex, T-cell receptor complex, antigen binding, and immune receptor activity (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003eA, B). Additionally, gene set enrichment analysis (GSEA) revealed the aggregation of gene sets related to the interferon gamma/alpha response and inflammation within the high-risk group, indicating a potential association between these pathways and aging-induced cancer progression (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003eC). These findings provide valuable insights into the molecular underpinnings of the impact of aging on cancer development, emphasizing the intricate involvement of immune-related processes and signaling pathways. To comprehensively assess molecular functional disparities between patients in the high- and low-risk groups, 16 cancer-related pathway activities were calculated[\u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e35\u003c/span\u003e, \u003cspan citationid=\"CR36\" class=\"CitationRef\"\u003e36\u003c/span\u003e]. This study revealed distinctive patterns of cell cycle dynamics, known for their impact on tissue regeneration, function, inflammation, and tumorigenesis. Specifically, patients in the high-risk group exhibited elevated cycle signature scores, indicative of increased activity within the cell cycle [\u003cspan citationid=\"CR37\" class=\"CitationRef\"\u003e37\u003c/span\u003e] (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003eF). This observation suggests a potential association between heightened cell cycle activity and an elevated risk profile, offering valuable insights into the potential drivers of cancer progression influenced by aging[\u003cspan citationid=\"CR37\" class=\"CitationRef\"\u003e37\u003c/span\u003e].\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec19\" class=\"Section2\"\u003e \u003ch2\u003eAging-associated alterations in the tumor microenvironment\u003c/h2\u003e \u003cp\u003eIn an endeavor to meticulously investigate the distinct cellular compositions of high- and low-risk patients, we collected publicly available single-cell RNA sequencing (scRNA-seq) data paired with bulk RNA-seq data from cancers. By integrating data from 24 samples with paired bulk and scRNA-seq information, we employed the mutual nearest neighbor (MNN) algorithm to mitigate batch effects. Rigorous quality control measures were applied to each individual sample, enabling the depiction of the cancer cellular landscape at single-cell resolution using uniform manifold approximation and projection (UMAP) visualization (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003eA). Using canonical lineage markers, we accurately annotated each cell population, which included epithelial cells, cycling cells, myeloid cells, T cells, B cells, plasma cells, cancer-associated fibroblasts (CAFs), endothelial cells, and pericytes (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003eB). For example, the expression of specific markers, such as CD79A, CD79B, and MS4A1 for B cells and EPCAM, KRT8, and KRT19 for epithelial cells, facilitated precise cell type identification[\u003cspan citationid=\"CR38\" class=\"CitationRef\"\u003e38\u003c/span\u003e, \u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e39\u003c/span\u003e]. The analysis of cellular subpopulations revealed distinctive compositions in the high- and low-risk groups, with significantly lower T-cell levels but higher frequencies of cycling cells, myeloid cells, and B cells in the high-risk group (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003eC). Further investigation of the functions of T cells within the high- to low-risk comparison revealed upregulation of the interferon γ and α response in these T cells. In contrast, the low-risk cohort exhibited a downregulation of mTORC1 signaling and TNFα signaling via NF-κB (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003eD-G). These findings suggest that aging-driven alterations in the tumor microenvironment contribute to cancer susceptibility [\u003cspan citationid=\"CR40\" class=\"CitationRef\"\u003e40\u003c/span\u003e].\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec20\" class=\"Section2\"\u003e \u003ch2\u003eMTHFD2 as an aging-associated factor in cancer\u003c/h2\u003e \u003cp\u003eIn the aforementioned study, MTHFD2 emerged as both a factor accentuating aging vulnerability and a protective factor against cancer progression. Consequently, MTHFD2 was chosen for a more in-depth exploration of its intricate role in aging-related cancer. To comprehensively describe the expression of MTHFD2 across a spectrum of tissues and its correlation with clinicopathological features in cancer patients, we conducted an exhaustive exploration employing an online database. MTHFD2 was markedly upregulated in 31 tumors compared with 2 tumors (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003eA). The discernible expression patterns of MTHFD2 were found to be intricately linked with modifications in ten genes, including TP53, CDH1, and MAP3KI, whereas the mutation status significantly differed between the high- and low-expression groups of MTHFD2 (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003eB). Single-cell analysis provided additional insights, showing that MTHFD2 was expressed mainly in T cells and myeloid subsets (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003eC). These simultaneous alterations in gene expression provide invaluable insights into the potential functions of MTHFD2. We conducted GO and GSEA analyses of the differentially expressed genes in the MTHFD2 high- and low-expression cohorts. These results indicated that MTHFD2 participates in nuclear division and meiotic cell cycle processes, indicating that MTHFD2 regulates cell cycle dynamics (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003eE-G). GSEA further confirmed this finding, revealing upregulated pathways, including E2F targets and G2M checkpoints (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003eH), which are integral to mechanisms regulating the cell cycle. Given the pivotal role of the cell cycle in preserving cellular homeostasis and its intricate linkage to the aging process[\u003cspan citationid=\"CR41\" class=\"CitationRef\"\u003e41\u003c/span\u003e], our findings suggest that MTHFD2 may contribute to aging-related modifications involved in cell cycle regulation. Subsequently, by analyzing the expression level of MTHFD2 in the young group and the aging group, it was found that the expression level of MTHFD2 was greater in the young group (Fig.\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003eI), which further confirmed that MTHFD2 may be involved in aging-related modifications.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec21\" class=\"Section2\"\u003e \u003ch2\u003eMTHFD2 loss drives senescence-like alterations to foster tumor growth\u003c/h2\u003e \u003cp\u003eIn our meticulous exploration of the intricate interplay between MTHFD2 and tumor growth during the aging process, we initiated MTHFD2 knockdown experiments employing B16F10 and MC38 cells as our experimental model [\u003cspan citationid=\"CR42\" class=\"CitationRef\"\u003e42\u003c/span\u003e] (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003eA). Strikingly, we observed notable metamorphoses in cell morphology within the MTHFD2 knockdown group, which appeared enlarged and flattened in shape, with increased cytoplasmic granularity, vacuolization, and altered nuclear morphology. Simultaneously, the levels of β-galactosidase, which serves as a biomarker for senescence, were measured in cells subjected to MTHFD2 knockdown and CDDP treatment (positive control) [\u003cspan citationid=\"CR43\" class=\"CitationRef\"\u003e43\u003c/span\u003e, \u003cspan citationid=\"CR44\" class=\"CitationRef\"\u003e44\u003c/span\u003e] (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003eC). Moreover, the rate of cell proliferation significantly decreased under MTHFD2 knockdown conditions (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003eB). Additionally, the perturbation of the cell cycle in B16F10 and MC38 cells was meticulously examined through flow cytometry, which revealed a pronounced blockade in the G0/G1 phase within the MTHFD2 knockdown group (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003eD). This finding strongly suggested a potential link between MTHFD2 manipulation and impediments in the progression of the cell cycle. Concurrently, the mRNA and protein levels of P21, a pivotal checkpoint in the cell cycle, were increased in the MTHFD2 knockdown group (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003eE). Consistent results confirmed the presence of senescence-associated secretory factors such as VEGF in the MTHFD2 knockdown group (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003eF). We investigated this possibility by employing B16F10 cells with MTHFD2 knockdown in a subcutaneous tumor model. In contrast, MTHFD2 knockdown promoted faster tumor growth (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003eG, H). Immunohistochemical analysis revealed significantly greater P21, VEGF, and IL-8 staining in tumor tissues from the MTHFD2 knockdown group than in those from the control group. Additionally, there was a notable increase in CD8\u003csup\u003e+\u003c/sup\u003e T-cell staining (Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003eI). This multifaceted approach not only underscores the potential involvement of MTHFD2 loss in propelling a senescent state but also highlights its potential significance in modulating tumor behavior.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e"},{"header":"Discussion","content":"\u003cp\u003eAging, an intricate biological process, is intricately intertwined with increased susceptibility to senescence-related diseases, notably cancer. Our detailed examination of senescence-related disparities across diverse cancer types revealed age to be a pivotal determinant of cancer prognosis, with individuals aged 60 years and older exhibiting increased vulnerability. By summarizing four different aging-related databases and using the Cox-LASSO algorithm, we obtained aging-associated gene alterations and aging-related risk scores, coupled with the integration of the TMB and exploration of oncogenic signaling pathways, which enriches our understanding of the molecular intricacies shaping cancer outcomes. The dissection of the cellular landscape within different risk groups reveals the complexities of the tumor microenvironment, providing a nuanced perspective. We identified MTHFD2 as both an aging risk factor and a protective factor in aging-related cancer.\u003c/p\u003e \u003cp\u003eBy analyzing the cell composition of the high- and low-risk groups, it was found that the T-cell composition of the high-risk group was significantly lower than that of the low-risk group, and T-cell exhaustion plays a pivotal role in rendering older individuals vulnerable to infections and cancer[\u003cspan citationid=\"CR45\" class=\"CitationRef\"\u003e45\u003c/span\u003e]. This decline is responsible for elevated susceptibility to infection and cancer[\u003cspan citationid=\"CR46\" class=\"CitationRef\"\u003e46\u003c/span\u003e]. Further investigations into the functions of T cells within the high- to low-risk group revealed the upregulation of the interferon γ and α response in these T cells, and the roles of IFN-γ and IFN-α are complex. While they contribute to immune surveillance against cancer, persistent activation or dysregulation of these interferons may also play a role in chronic inflammation, which is a hallmark of aging and a risk factor for cancer[\u003cspan citationid=\"CR47\" class=\"CitationRef\"\u003e47\u003c/span\u003e, \u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e48\u003c/span\u003e].\u003c/p\u003e \u003cp\u003eMethylenetetrahydrofolate dehydrogenase 2 (MTHFD2) is a crucial enzyme involved in cellular metabolism, particularly in the folate metabolic pathway. This enzyme plays a profound role in nucleotide synthesis, which is essential for DNA replication and cellular proliferation. MTHFD2 catalyzes the conversion of 5,10-methylenetetrahydrofolate to 5,10-methenyltetrahydrofolate, an essential step in the synthesis of purines, amino acids, and other important molecules within the cell. MTHFD2 has garnered attention for its involvement in various cellular processes linked to cancer progression. Multiple investigations have underscored the propensity for increased MTHFD2 expression across diverse cancer types, including breast, lung, colorectal, and pancreatic malignancies[\u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e19\u003c/span\u003e]. Notably, MTHFD2 is preferentially upregulated in undifferentiated or poorly differentiated tumors[\u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e49\u003c/span\u003e], suggesting its cancer-specific expression pattern. In addition to its metabolic functions, MTHFD2 has been implicated in promoting cancer immune evasion. By steering the folate cycle toward the maintenance of adequate UDP-GlcNAc levels, MTHFD2 propels the O-GlcNAcylation of cMYC, thereby bolstering cMYC stability and PD-L1 transcription.[\u003cspan citationid=\"CR50\" class=\"CitationRef\"\u003e50\u003c/span\u003e] Nevertheless, the functional repertoire of MTHFD2 is intricate and multifaceted. Although MTHFD2 is primarily expressed within the mitochondria, it is also present in the nucleus. TH9619 effectively inhibits both the dehydrogenase and cyclohydrolase activities of MTHFD1/2. MTHFD1, found in the cytoplasm, shares similar functions with MTHFD2. TH9619 inhibits MTHFD1 activity, preventing the incorporation of formate produced by the mitochondria into dTMP. It selectively targets nuclear MTHFD2 without affecting mitochondrial MTHFD2, leading to an overflow of formate from the mitochondria. Consequently, cancer cells undergo cell death despite high MTHFD2 expression[\u003cspan citationid=\"CR51\" class=\"CitationRef\"\u003e51\u003c/span\u003e]. In addition, recent studies have indicated that MTHFD2 inhibition induces apoptosis solely in UQCR11-null cells, while it does not affect UQCR11-intact cells[\u003cspan citationid=\"CR52\" class=\"CitationRef\"\u003e52\u003c/span\u003e]. This suggests that the apoptotic effect of MTHFD2 inhibition is specific to cancer cells. Inversely, under conditions of oxygen deprivation or impaired electron transport chain activity, MTHFD2 maintains a substantial supply of NADH through the promotion of serine catabolism, leading to cell death. In hypoxic cells with compromised respiration, inhibition of MTHFD2 partially restores NADH levels and facilitates cell proliferation[\u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e49\u003c/span\u003e]. These findings suggest a potential anticancer role for MTHFD2. Our findings reveal that deletion of MTHFD2 contributes to cellular senescence. Similarly, prior research has revealed that the colocalization of MTHFD2 with DNA replication sites in the nucleus promotes cell cycle progression. Deletion of MTHFD2 leads to S-phase cell cycle arrest and fosters a senescence-like state[\u003cspan citationid=\"CR53\" class=\"CitationRef\"\u003e53\u003c/span\u003e]. Additionally, MTHFD2 is expressed in developing embryos but is notably absent in most healthy adult tissues, including proliferating ones[\u003cspan citationid=\"CR54\" class=\"CitationRef\"\u003e54\u003c/span\u003e]. This observation implies a decrease in MTHFD2 expression with age. MTHFD2 deletion increases the infiltration of CD8\u003csup\u003e+\u003c/sup\u003e T cells in B16F10 tumors. The potential mechanisms are manifold. Initially, senescent cancer cells manifest hyperploidy, rendering them genomically unstable and enabling the presentation of tumor antigens to activate immunosurveillance, which entails the recruitment of immune effectors such as B, NK, NKT, and T cells[\u003cspan citationid=\"CR55\" class=\"CitationRef\"\u003e55\u003c/span\u003e, \u003cspan citationid=\"CR56\" class=\"CitationRef\"\u003e56\u003c/span\u003e]. Furthermore, the senescence-associated secretory phenotype (SASP) factor VEGF fosters angiogenesis, thereby facilitating T cell infiltration. Nonetheless, the buildup of senescent cells in tumor-bearing mice precipitates T-cell exhaustion[\u003cspan citationid=\"CR57\" class=\"CitationRef\"\u003e57\u003c/span\u003e]. For instance, IL-8, elevated in MTHFD2-depleted cells, augments PD-1 expression in CD8\u0026thinsp;+\u0026thinsp;T cells, leading to T-cell exhaustion. Therefore, the deletion of MTHFD2 facilitates the infiltration of CD8\u0026thinsp;+\u0026thinsp;T cells into tumors but compromises their functional activity. Hence, we propose a novel perspective suggesting that MTHFD2 exerts its anticancer effect by retarding the aging process.\u003c/p\u003e \u003cp\u003eHowever, some limitations are noted in this study. First, numerous studies have demonstrated that aging is a significant contributor to the development and progression of cancer. Based on data analysis related to breast cancer and animal experiments, this study revealed that downregulation of MTHFD2 can promote cellular senescence, thereby accelerating tumor growth. As such, our study may lack sufficient novelty. Second, the mechanisms involved in this study were not further investigated. For instance, we did not explore how downregulation of MTHFD2 promotes cellular senescence and its impact on tumor immunity during this process. The prognostic role of MTHFD2 in cancer patients requires validation in additional cohorts to confirm its prognostic significance. Finally, the exact effect of MTHFD2 downregulation on the tumor microenvironment remains unclear due to the lack of single-cell sequencing or multi-parameter flow cytometry analyses.\u003c/p\u003e"},{"header":"Conclusions","content":"\u003cp\u003eIn summary, the process of aging leads to heightened genetic mutations within tumors, triggers the activation of pathways essential for tumor proliferation and spread, and promotes an immunosuppressive microenvironment. The identification of Mthfd2 as a suppressor in tumors associated with aging implies that maintaining its expression and functionality may offer a potential therapeutic approach for addressing such tumors.\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003eEthics approval and consent to participate\u003c/p\u003e\n\u003cp\u003eThe animal experiment of the study was approved by\u0026nbsp;the Institutional Animal Care and Use Committee of Tongji University Cancer Center, Shanghai Tenth People\u0026apos;s Hospital, School of Medicine, Tongji University.\u0026nbsp;(Issue No. 22KN151).\u003c/p\u003e\n\u003cp\u003eConsent for publication\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eNot applicable.\u003c/p\u003e\n\u003cp\u003eAvailability of data and materials\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eThe original contributions presented in the study are included in the article/additional file.\u003c/p\u003e\n\u003cp\u003eCompeting interests\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eThe author (s) declare no potential conficts of interest with respect to the research, authorship, and/or publication of this article.\u003c/p\u003e\n\u003cp\u003eFunding\u003c/p\u003e\n\u003cp\u003eThis project was supported by the National Natural Science Foundation of China (82203629 and 82372780), Shanghai Pujiang Program (22PJD054 and 22PJ1412000) and The Fundamental Research Funds for the Central University (22120240320).\u003c/p\u003e\n\u003cp\u003eAuthor contribution:\u003c/p\u003e\n\u003cp\u003ePW designed the study. PW, ZF, and WP performed the experiments. PW and ZF analyzed the data. PW drafted the manuscript. QW revised the manuscript. All authors read and approved the final manuscript.\u003c/p\u003e\n\u003cp\u003eAcknowledgements\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eWe are grateful to all the editors and reviewers for their hard work, and all the authors for their efforts.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\n\u003cli\u003eL\u0026oacute;pez-Ot\u0026iacute;n C, Blasco MA, Partridge L, Serrano M, Kroemer G. The Hallmarks of Aging. Cell. 2013;153:1194\u0026ndash;217.\u003c/li\u003e\n\u003cli\u003eNiccoli T, Partridge L. Ageing as a Risk Factor for Disease. Curr Biol. 2012;22:R741\u0026ndash;52.\u003c/li\u003e\n\u003cli\u003eVijg J, Suh Y. Genetics of Longevity and Aging. Annu Rev Med. 2005;56:193\u0026ndash;212.\u003c/li\u003e\n\u003cli\u003eCampisi J. Aging, Cellular Senescence, and Cancer. Annu Rev Physiol. 2013;75:685\u0026ndash;705.\u003c/li\u003e\n\u003cli\u003eHurria A, Jones L, Muss HB. Cancer Treatment as an Accelerated Aging Process: Assessment, Biomarkers, and Interventions. 2016;\u003c/li\u003e\n\u003cli\u003eHanahan D, Weinberg RA. Hallmarks of Cancer: The Next Generation. Cell. 2011;144:646\u0026ndash;74.\u003c/li\u003e\n\u003cli\u003eFranceschi C, Garagnani P, Parini P, Giuliani C, Santoro A. Inflammaging: a new immune\u0026ndash;metabolic viewpoint for age-related diseases. Nat Rev Endocrinol. 2018;14:576\u0026ndash;90.\u003c/li\u003e\n\u003cli\u003eLionaki E, Ploumi C, Tavernarakis N. One-Carbon Metabolism: Pulling the Strings behind Aging and Neurodegeneration. Cells. 2022;11:214.\u003c/li\u003e\n\u003cli\u003eTedeschi PM, Markert EK, Gounder M, Lin H, Dvorzhinski D, Dolfi SC, et al. Contribution of serine, folate and glycine metabolism to the ATP, NADPH and purine requirements of cancer cells. Cell Death Dis. 2013;4:e877\u0026ndash;e877.\u003c/li\u003e\n\u003cli\u003eYang M, Vousden KH. Serine and one-carbon metabolism in cancer. Nat Rev Cancer. 2016;16:650\u0026ndash;62.\u003c/li\u003e\n\u003cli\u003eDucker GS, Chen L, Morscher RJ, Ghergurovich JM, Esposito M, Teng X, et al. Reversal of Cytosolic One-Carbon Flux Compensates for Loss of the Mitochondrial Folate Pathway. Cell Metab. 2016;23:1140\u0026ndash;53.\u003c/li\u003e\n\u003cli\u003eNilsson R, Jain M, Madhusudhan N, Sheppard NG, Strittmatter L, Kampf C, et al. Metabolic enzyme expression highlights a key role for MTHFD2 and the mitochondrial folate pathway in cancer. Nat Commun. 2014;5:3128.\u003c/li\u003e\n\u003cli\u003eFan J, Ye J, Kamphorst JJ, Shlomi T, Thompson CB, Rabinowitz JD. Quantitative flux analysis reveals folate-dependent NADPH production. Nature. 2014;510:298\u0026ndash;302.\u003c/li\u003e\n\u003cli\u003eGreen NH, Galvan DL, Badal SS, Chang BH, LeBleu VS, Long J, et al. MTHFD2 links RNA methylation to metabolic reprogramming in renal cell carcinoma. Oncogene. 2019;38:6211\u0026ndash;25.\u003c/li\u003e\n\u003cli\u003eGalili T. dendextend: an R package for visualizing, adjusting and comparing trees of hierarchical clustering. Bioinformatics. 2015;31:3718\u0026ndash;20.\u003c/li\u003e\n\u003cli\u003eLawrence MS, Stojanov P, Polak P, Kryukov GV, Cibulskis K, Sivachenko A, et al. Mutational heterogeneity in cancer and the search for new cancer-associated genes. Nature. 2013;499:214\u0026ndash;8.\u003c/li\u003e\n\u003cli\u003eZhou M, Zhao H, Wang Z, Cheng L, Yang L, Shi H, et al. Identification and validation of potential prognostic lncRNA biomarkers for predicting survival in patients with multiple myeloma. J Exp Clin Cancer Res. 2015;34:102.\u003c/li\u003e\n\u003cli\u003eTang Z, Li C, Kang B, Gao G, Li C, Zhang Z. GEPIA: a web server for cancer and normal gene expression profiling and interactive analyses. Nucleic Acids Res. 2017;45:W98\u0026ndash;102.\u003c/li\u003e\n\u003cli\u003eNilsson R, Jain M, Madhusudhan N, Sheppard NG, Strittmatter L, Kampf C, et al. Metabolic enzyme expression highlights a key role for MTHFD2 and the mitochondrial folate pathway in cancer. Nat Commun. 2014;5:3128.\u003c/li\u003e\n\u003cli\u003eYang Y, Lane AN, Ricketts CJ, Sourbier C, Wei M-H, Shuch B, et al. Metabolic Reprogramming for Producing Energy and Reducing Power in Fumarate Hydratase Null Cells from Hereditary Leiomyomatosis Renal Cell Carcinoma. Tan M, editor. PLoS ONE. 2013;8:e72179.\u003c/li\u003e\n\u003cli\u003eZhang B, Zheng A, Hydbring P, Ambroise G, Ouchida AT, Goiny M, et al. PHGDH Defines a Metabolic Subtype in Lung Adenocarcinomas with Poor Prognosis. Cell Rep. 2017;19:2289\u0026ndash;303.\u003c/li\u003e\n\u003cli\u003eBalducci L, Ershler WB. Cancer and ageing: a nexus at several levels. Nat Rev Cancer. 2005;5:655\u0026ndash;62.\u003c/li\u003e\n\u003cli\u003eSaul D, Kosinsky RL, Atkinson EJ, Doolittle ML, Zhang X, LeBrasseur NK, et al. A new gene set identifies senescent cells and predicts senescence-associated pathways across tissues. Nat Commun. 2022;13:4827.\u003c/li\u003e\n\u003cli\u003eGorgoulis V, Adams PD, Alimonti A, Bennett DC, Bischof O, Bishop C, et al. Cellular Senescence: Defining a Path Forward. Cell. 2019;179:813\u0026ndash;27.\u003c/li\u003e\n\u003cli\u003eHurria A, Togawa K, Mohile SG, Owusu C, Klepin HD, Gross CP, et al. Predicting Chemotherapy Toxicity in Older Adults With Cancer: A Prospective Multicenter Study. J Clin Oncol. 2011;29:3457\u0026ndash;65.\u003c/li\u003e\n\u003cli\u003eMuller PAJ, Vousden KH. p53 mutations in cancer. Nat Cell Biol. 2013;15:2\u0026ndash;8.\u003c/li\u003e\n\u003cli\u003eFutreal PA, Coin L, Marshall M, Down T, Hubbard T, Wooster R, et al. A census of human cancer genes. Nat Rev Cancer. 2004;4:177\u0026ndash;83.\u003c/li\u003e\n\u003cli\u003eKandoth C, McLellan MD, Vandin F, Ye K, Niu B, Lu C, et al. Mutational landscape and significance across 12 major cancer types. Nature. 2013;502:333\u0026ndash;9.\u003c/li\u003e\n\u003cli\u003eSha D, Jin Z, Budzcies J, Kluck K, Stenzinger A, Sinicrope FA. Tumor Mutational Burden (TMB) as a Predictive Biomarker in Solid Tumors. Cancer Discov. 2020;10:1808\u0026ndash;25.\u003c/li\u003e\n\u003cli\u003eSamstein RM, Lee C-H, Shoushtari AN, Hellmann MD, Shen R, Janjigian YY, et al. Tumor mutational load predicts survival after immunotherapy across multiple cancer types. Nat Genet. 2019;51:202\u0026ndash;6.\u003c/li\u003e\n\u003cli\u003eChalmers ZR, Connelly CF, Fabrizio D, Gay L, Ali SM, Ennis R, et al. Analysis of 100,000 human cancer genomes reveals the landscape of tumor mutational burden. Genome Med. 2017;9:34.\u003c/li\u003e\n\u003cli\u003eChan TA, Yarchoan M, Jaffee E, Swanton C, Quezada SA, Stenzinger A, et al. Development of tumor mutation burden as an immunotherapy biomarker: utility for the oncology clinic. Ann Oncol. 2019;30:44\u0026ndash;56.\u003c/li\u003e\n\u003cli\u003eHellmann MD, Nathanson T, Rizvi H, Creelan BC, Sanchez-Vega F, Ahuja A, et al. Genomic Features of Response to Combination Immunotherapy in Patients with Advanced Non-Small-Cell Lung Cancer. Cancer Cell. 2018;33:843-852.e4.\u003c/li\u003e\n\u003cli\u003eHellmann MD, Ciuleanu T-E, Pluzanski A, Lee JS, Otterson GA, Audigier-Valette C, et al. Nivolumab plus Ipilimumab in Lung Cancer with a High Tumor Mutational Burden. N Engl J Med. 2018;378:2093\u0026ndash;104.\u003c/li\u003e\n\u003cli\u003eHanahan D, Weinberg RA. Hallmarks of Cancer: The Next Generation. Cell. 2011;144:646\u0026ndash;74.\u003c/li\u003e\n\u003cli\u003eVogelstein B, Papadopoulos N, Velculescu VE, Zhou S, Diaz LA, Kinzler KW. Cancer Genome Landscapes. Science. 2013;339:1546\u0026ndash;58.\u003c/li\u003e\n\u003cli\u003eMalumbres M, Barbacid M. Cell cycle, CDKs and cancer: a changing paradigm. Nat Rev Cancer. 2009;9:153\u0026ndash;66.\u003c/li\u003e\n\u003cli\u003ePalanichamy A, Apeltsin L, Kuo TC, Sirota M, Wang S, Pitts SJ, et al. Immunoglobulin class-switched B cells form an active immune axis between CNS and periphery in multiple sclerosis. Sci Transl Med [Internet]. 2014 [cited 2024 Jan 11];6. Available from: https://www.science.org/doi/10.1126/scitranslmed.3008930\u003c/li\u003e\n\u003cli\u003eFalini B, Fizzotti M, Pucciarini A, Bigerna B, Marafioti T, Gambacorta M, et al. A monoclonal antibody (MUM1p) detects expression of the MUM1/IRF4 protein in a subset of germinal center B cells, plasma cells, and activated T cells. Blood. 2000;95:2084\u0026ndash;92.\u003c/li\u003e\n\u003cli\u003eDe Visser KE, Joyce JA. The evolving tumor microenvironment: From cancer initiation to metastatic outgrowth. Cancer Cell. 2023;41:374\u0026ndash;403.\u003c/li\u003e\n\u003cli\u003eHartl FU. Cellular Homeostasis and Aging. Annu Rev Biochem. 2016;85:1\u0026ndash;4.\u003c/li\u003e\n\u003cli\u003eBoccaletto P, Machnicka MA, Purta E, Piątkowski P, Bagiński B, Wirecki TK, et al. MODOMICS: a database of RNA modification pathways. 2017 update. Nucleic Acids Res. 2018;46:D303\u0026ndash;7.\u003c/li\u003e\n\u003cli\u003eLee BY, Han JA, Im JS, Morrone A, Johung K, Goodwin EC, et al. Senescence‐associated \u0026beta;‐galactosidase is lysosomal \u0026beta;‐galactosidase. Aging Cell. 2006;5:187\u0026ndash;95.\u003c/li\u003e\n\u003cli\u003eDebacq-Chainiaux F, Erusalimsky JD, Campisi J, Toussaint O. Protocols to detect senescence-associated beta-galactosidase (SA-\u0026beta;gal) activity, a biomarker of senescent cells in culture and in vivo. Nat Protoc. 2009;4:1798\u0026ndash;806.\u003c/li\u003e\n\u003cli\u003eGoronzy JJ, Weyand CM. Mechanisms underlying T cell ageing. Nat Rev Immunol. 2019;19:573\u0026ndash;83.\u003c/li\u003e\n\u003cli\u003eMittelbrunn M, Kroemer G. Hallmarks of T cell aging. Nat Immunol. 2021;22:687\u0026ndash;98.\u003c/li\u003e\n\u003cli\u003eTakaoka A, Mitani Y, Sato M, Yokochi T, Noguchi S, Tanaka N, et al. Cross Talk Between Interferon-y and - d P Signaling Components in Caveolar Membrane Domains. 2000;288.\u003c/li\u003e\n\u003cli\u003eCoelho MA, Cooper S, Strauss ME, Karakoc E, Bhosle S, Gon\u0026ccedil;alves E, et al. Base editing screens map mutations affecting interferon-\u0026gamma; signaling in cancer. Cancer Cell. 2023;41:288-303.e6.\u003c/li\u003e\n\u003cli\u003eYang L, Canaveras JCG, Chen Z, Wang L, Liang L, Jang C, et al. Serine Catabolism Feeds NADH when Respiration Is Impaired. Cell Metab. 2020;31:809-821.e6.\u003c/li\u003e\n\u003cli\u003eShang M, Yang H, Yang R, Chen T, Fu Y, Li Y, et al. The folate cycle enzyme MTHFD2 induces cancer immune evasion through PD-L1 up-regulation. Nat Commun. 2021;12:1940.\u003c/li\u003e\n\u003cli\u003eGreen AC, Marttila P, Kiweler N, Chalkiadaki C, Wiita E, Cookson V, et al. Formate overflow drives toxic folate trapping in MTHFD1 inhibited cancer cells. Nat Metab. 2023;5:642\u0026ndash;59.\u003c/li\u003e\n\u003cli\u003eAchreja A, Yu T, Mittal A, Choppara S, Animasahun O, Nenwani M, et al. Metabolic collateral lethal target identification reveals MTHFD2 paralogue dependency in ovarian cancer. Nat Metab. 2022;4:1119\u0026ndash;37.\u003c/li\u003e\n\u003cli\u003eBonagas N, Gustafsson NMS, Henriksson M, Marttila P, Gustafsson R, Wiita E, et al. Pharmacological targeting of MTHFD2 suppresses acute myeloid leukemia by inducing thymidine depletion and replication stress. Nat Cancer. 2022;3:156\u0026ndash;72.\u003c/li\u003e\n\u003cli\u003eNilsson R, Jain M, Madhusudhan N, Sheppard NG, Strittmatter L, Kampf C, et al. Metabolic enzyme expression highlights a key role for MTHFD2 and the mitochondrial folate pathway in cancer. Nat Commun. 2014;5:3128.\u003c/li\u003e\n\u003cli\u003eSenovilla L, Vitale I, Martins I, Tailler M, Pailleret C, Michaud M, et al. An Immunosurveillance Mechanism Controls Cancer Cell Ploidy. Science. 2012;337:1678\u0026ndash;84.\u003c/li\u003e\n\u003cli\u003eAranda F, Chaba K, Bloy N, Garcia P, Bordenave C, Martins I, et al. Immune effectors responsible for the elimination of hyperploid cancer cells. OncoImmunology. 2018;7:e1463947.\u003c/li\u003e\n\u003cli\u003eMaggiorani D, Le O, Lisi V, Landais S, Moquin-Beaudry G, Lavall\u0026eacute;e VP, et al. Senescence drives immunotherapy resistance by inducing an immunosuppressive tumor microenvironment. Nat Commun. 2024;15:2435.\u003c/li\u003e\n\u003c/ol\u003e"},{"header":"Tables","content":"\u003cp\u003eTable1：Clinical information of breast cancer patients\u003c/p\u003e\n\u003ctable border=\"1\" cellspacing=\"0\" cellpadding=\"0\"\u003e\n \u003ctbody\u003e\n \u003ctr\u003e\n \u003ctd width=\"26.173285198555956%\" valign=\"top\"\u003e\n \u003cp\u003eVariable\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"24.36823104693141%\" valign=\"top\"\u003e\n \u003cp\u003e\u0026nbsp;\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"24.729241877256317%\" valign=\"top\"\u003e\n \u003cp\u003eNo recurrence\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"24.729241877256317%\" valign=\"top\"\u003e\n \u003cp\u003eRecurrence\u0026nbsp; \u0026nbsp; \u0026nbsp;\u0026nbsp;\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd width=\"26.173285198555956%\" rowspan=\"2\" valign=\"top\"\u003e\n \u003cp\u003eAge at diagnosis, years\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"24.36823104693141%\" valign=\"top\"\u003e\n \u003cp\u003e\u0026le;60\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"24.729241877256317%\" valign=\"top\"\u003e\n \u003cp\u003e879\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"24.729241877256317%\" valign=\"top\"\u003e\n \u003cp\u003e110\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd width=\"33.007334963325185%\" valign=\"top\"\u003e\n \u003cp\u003e\u0026gt;60\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"33.49633251833741%\" valign=\"top\"\u003e\n \u003cp\u003e209\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"33.49633251833741%\" valign=\"top\"\u003e\n \u003cp\u003e57\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd width=\"26.173285198555956%\" rowspan=\"3\" valign=\"top\"\u003e\n \u003cp\u003eHistopathological grade\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"24.36823104693141%\" valign=\"top\"\u003e\n \u003cp\u003e1\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"24.729241877256317%\" valign=\"top\"\u003e\n \u003cp\u003e55\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"24.729241877256317%\" valign=\"top\"\u003e\n \u003cp\u003e1\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd width=\"33.007334963325185%\" valign=\"top\"\u003e\n \u003cp\u003e2\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"33.49633251833741%\" valign=\"top\"\u003e\n \u003cp\u003e307\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"33.49633251833741%\" valign=\"top\"\u003e\n \u003cp\u003e56\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd width=\"33.007334963325185%\" valign=\"top\"\u003e\n \u003cp\u003e3\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"33.49633251833741%\" valign=\"top\"\u003e\n \u003cp\u003e165\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"33.49633251833741%\" valign=\"top\"\u003e\n \u003cp\u003e11\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd width=\"26.173285198555956%\" rowspan=\"2\" valign=\"top\"\u003e\n \u003cp\u003eT stage\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"24.36823104693141%\" valign=\"top\"\u003e\n \u003cp\u003eT1+T2\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"24.729241877256317%\" valign=\"top\"\u003e\n \u003cp\u003e693\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"24.729241877256317%\" valign=\"top\"\u003e\n \u003cp\u003e72\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd width=\"33.007334963325185%\" valign=\"top\"\u003e\n \u003cp\u003eT3+T4\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"33.49633251833741%\" valign=\"top\"\u003e\n \u003cp\u003e48\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"33.49633251833741%\" valign=\"top\"\u003e\n \u003cp\u003e7\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd width=\"26.173285198555956%\" rowspan=\"2\" valign=\"top\"\u003e\n \u003cp\u003eN stage\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"24.36823104693141%\" valign=\"top\"\u003e\n \u003cp\u003eN0\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"24.729241877256317%\" valign=\"top\"\u003e\n \u003cp\u003e477\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"24.729241877256317%\" valign=\"top\"\u003e\n \u003cp\u003e40\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd width=\"33.007334963325185%\" valign=\"top\"\u003e\n \u003cp\u003eN1+N2+N3\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"33.49633251833741%\" valign=\"top\"\u003e\n \u003cp\u003e57\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"33.49633251833741%\" valign=\"top\"\u003e\n \u003cp\u003e39\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003c/tbody\u003e\n\u003c/table\u003e\n\u003cp\u003e\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eTable2：Abbreviations\u003c/p\u003e\n\u003ctable border=\"1\" cellspacing=\"0\" cellpadding=\"0\"\u003e\n \u003ctbody\u003e\n \u003ctr\u003e\n \u003ctd width=\"25.49728752260398%\" valign=\"top\"\u003e\n \u003cp\u003eSASP\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"74.50271247739602%\" valign=\"top\"\u003e\n \u003cp\u003eSenescence-associated secretory phenotype\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd width=\"25.49728752260398%\" valign=\"top\"\u003e\n \u003cp\u003eMTHFD2\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"74.50271247739602%\" valign=\"top\"\u003e\n \u003cp\u003eMethylenetetrahydrofolate\u0026nbsp;dehydrogenase\u0026nbsp;2\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd width=\"25.49728752260398%\" valign=\"top\"\u003e\n \u003cp\u003eSCNAs\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"74.50271247739602%\" valign=\"top\"\u003e\n \u003cp\u003eSomatic copy\u0026nbsp;number alterations\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd width=\"25.49728752260398%\" valign=\"top\"\u003e\n \u003cp\u003eMAF\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"74.50271247739602%\" valign=\"top\"\u003e\n \u003cp\u003eMutation Annotation Format\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd width=\"25.49728752260398%\" valign=\"top\"\u003e\n \u003cp\u003eGDC\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"74.50271247739602%\" valign=\"top\"\u003e\n \u003cp\u003eGenomic Data Commons\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd width=\"25.49728752260398%\" valign=\"top\"\u003e\n \u003cp\u003eTMB\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"74.50271247739602%\" valign=\"top\"\u003e\n \u003cp\u003eTumor Mutation Burden\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd width=\"25.49728752260398%\" valign=\"top\"\u003e\n \u003cp\u003eDEGs\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"74.50271247739602%\" valign=\"top\"\u003e\n \u003cp\u003eDifferentially expressed genes\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd width=\"25.49728752260398%\" valign=\"top\"\u003e\n \u003cp\u003eFDR\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"74.50271247739602%\" valign=\"top\"\u003e\n \u003cp\u003eFalse discovery rate\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd width=\"25.49728752260398%\" valign=\"top\"\u003e\n \u003cp\u003eGO\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"74.50271247739602%\" valign=\"top\"\u003e\n \u003cp\u003eGene Ontology\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd width=\"25.49728752260398%\" valign=\"top\"\u003e\n \u003cp\u003eKEGG\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"74.50271247739602%\" valign=\"top\"\u003e\n \u003cp\u003eKyoto Encyclopedia of Genes and Genomes\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd width=\"25.49728752260398%\" valign=\"top\"\u003e\n \u003cp\u003eGSEA\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"74.50271247739602%\" valign=\"top\"\u003e\n \u003cp\u003eGene set enrichment analysis\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd width=\"25.49728752260398%\" valign=\"top\"\u003e\n \u003cp\u003escRNA-seq\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"74.50271247739602%\" valign=\"top\"\u003e\n \u003cp\u003eSingle-cell RNA sequencing\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd width=\"25.49728752260398%\" valign=\"top\"\u003e\n \u003cp\u003eMNN\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"74.50271247739602%\" valign=\"top\"\u003e\n \u003cp\u003eMutual nearest neighbor\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd width=\"25.49728752260398%\" valign=\"top\"\u003e\n \u003cp\u003eUMAP\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"74.50271247739602%\" valign=\"top\"\u003e\n \u003cp\u003eUniform manifold approximation and projection\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003ctr\u003e\n \u003ctd width=\"25.49728752260398%\" valign=\"top\"\u003e\n \u003cp\u003eCAFs\u003c/p\u003e\n \u003c/td\u003e\n \u003ctd width=\"74.50271247739602%\" valign=\"top\"\u003e\n \u003cp\u003eCancer-associated fibroblasts\u003c/p\u003e\n \u003c/td\u003e\n \u003c/tr\u003e\n \u003c/tbody\u003e\n\u003c/table\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":true,"highlight":"","institution":"","isAcceptedByJournal":false,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true},"keywords":"Senescence, Tumor, Cell cycle, MTHFD2","lastPublishedDoi":"10.21203/rs.3.rs-4353045/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-4353045/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003ch2\u003eBackground\u003c/h2\u003e \u003cp\u003eAge is a critical risk factor for cancer, as its incidence and mortality increase with age. However, there's limited understanding of the molecular changes aging induces in tumors.\u003c/p\u003e\u003ch2\u003eMethods\u003c/h2\u003e \u003cp\u003eWe explored demographic differences between young and old cancer patients and identified age sixty and above as pivotal in cancer prognosis. Subsequently, we developed an aging-related prognostic model based on genes to assess senescence's impact on aging-associated cancer.\u003c/p\u003e\u003ch2\u003eResults\u003c/h2\u003e \u003cp\u003eOur study revealed increased genomic instability and somatic mutations in tumors from older individuals. We also found alterations in carcinogenic signaling pathways, particularly immune responses, inflammatory pathways, and cell cycle arrest in susceptible populations. Single-cell RNA sequencing showed heightened frequencies of exhausted T cells, myeloid cells, and B cells in high-risk cohorts.\u003c/p\u003e\u003ch2\u003eConclusion\u003c/h2\u003e \u003cp\u003eMTHFD2 emerged as a crucial molecular switch regulating senescence in cancer. Its deletion promoted tumor growth by inducing cell senescence and stimulating the senescence-associated secretory phenotype (SASP) in senescent tumor cells. This highlights the need for tailored methodologies in effective cancer management.\u003c/p\u003e","manuscriptTitle":"Senescence Reprogramming by MTHFD2 Deficiency Facilitates Tumor Progression","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2024-05-14 17:46:03","doi":"10.21203/rs.3.rs-4353045/v1","editorialEvents":[{"type":"communityComments","content":0}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"793398ec-23c5-4d46-811d-b2fc7e936b80","owner":[],"postedDate":"May 14th, 2024","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"posted","subjectAreas":[],"tags":[],"updatedAt":"2024-05-29T22:53:18+00:00","versionOfRecord":[],"versionCreatedAt":"2024-05-14 17:46:03","video":"","vorDoi":"","vorDoiUrl":"","workflowStages":[]},"version":"v1","identity":"rs-4353045","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-4353045","identity":"rs-4353045","version":["v1"]},"buildId":"8U1c8b4HqxoKbykW_rLl7","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}