Senescent tumor cells in colorectal cancer are characterized by elevated enzymatic activity of complexes 1 and 2 in oxidative phosphorylation
Article information
Abstract
Background
Cellular senescence is defined as an irreversible cell cycle arrest caused by various internal and external insults. While the metabolic dysfunction of senescent cells in normal tissue is relatively well-established, there is a lack of information regarding the metabolic features of senescent tumor cells.
Methods
Publicly available single-cell RNA-sequencing data from the GSE166555 and GSE178341 datasets were utilized to investigate the metabolic features of senescent tumor cells. To validate the single-cell RNA-sequencing data, we performed senescence-associated β-galactosidase (SA-β-Gal) staining to identify senescent tumor cells in fresh frozen colorectal cancer tissue. We also evaluated nicotinamide adenine dinucleotide dehydrogenase–tetrazolium reductase (NADH-TR) and succinate dehydrogenase (SDH) activity using enzyme histochemical methods and compared the staining with SA-β-Gal staining. MTT assay was performed to reveal the complex 1 activity of the respiratory chain in in-vitro senescence model.
Results
Single-cell RNA-sequencing data revealed an upregulation in the activity of complexes 1 and 2 in oxidative phosphorylation, despite overall mitochondrial dysfunction in senescent tumor cells. Both SA-β-Gal and enzyme histochemical staining using fresh frozen colorectal cancer tissues indicated a high correlation between SA-β-Gal positivity and NADH-TR/SDH staining positivity. MTT assay showed that senescent colorectal cancer cells exhibit higher absorbance in 600 nm wavelength.
Conclusions
Senescent tumor cells exhibit distinct metabolic features, characterized by upregulation of complexes 1 and 2 in the oxidative phosphorylation pathway. NADH-TR and SDH staining represent efficient methods for detecting senescent tumor cells in colorectal cancer.
Cellular senescence is an irreversible cell cycle arrest first reported in normal fibroblasts after serial cultivation in vitro by Hayflick and Moorhead in 1961 [1]. Cellular senescence can be induced by a variety of factors, including telomere shortening or dysfunction, DNA damage, oncogene activation or loss of tumor suppressor functions, epigenetic changes, and organelle damage [2,3]. Moreover, cellular senescence can be induced by oncogenic activation, a process known as oncogene-induced senescence (OIS) [4]. Since OIS was initially observed in premalignant lesions, it was initially considered to serve as a barrier to malignant transformation [4]. However, recent studies have demonstrated that senescent malignant cells can still be identified in various types of cancers [5-8]. Furthermore, these senescent cancer cells have been shown to promote the development and/or progression of cancer [5-7].
Despite being in a state of permanent proliferation arrest, senescent cells remain metabolically active [9,10]. Indeed, senescent cells are characterized by its vigorous synthesis of specific proteins: various cytokines and release factors, which are called senescence-associated secretory phenotype (SASP) [11]; the presence of senescent cells may massively affect adjacent microenvironment by inducing inflammatory responses [11]. Therefore, senescent cells may need sufficient energy to fit their demand to synthesize various metabolites. Since mitochondrial dysfunction is primarily involved in cellular senescence, the efficiency of oxidative phosphorylation (OXPHOS) in senescent cells from normal tissues is significantly decreased [12]. Therefore, it is widely accepted that senescent cells maintain their energy demand by increasing glycolysis and fermentation [13]. However, multiple studies have revealed that the OXPHOS status in malignant cells can vary. Some studies have shown that cancer cells exhibit upregulated glycolysis compared to normal cells, and OXPHOS is commonly downregulated in various cancers [14,15], while others have shown that OXPHOS can increase [16-18]. Nevertheless, the status of OXPHOS in senescent tumor cells remains unknown. From this background, we hypothesized that the main energy production pathway of senescent tumor cells might be different from non-senescent tumor cells. To check this hypothesis, we utilized single-cell RNA-sequencing to reveal the metabolic feature of senescent tumor cells. Subsequently, enzymatic activity related with OXPHOS was examined; enzyme histochemical methods, such as nicotinamide adenine dinucleotide dehydrogenase–tetrazolium reductase (NADH-TR) and succinate dehydrogenase (SDH) staining, were applied. By visualizing the activity of these enzymes, NADH-TR and SDH staining can provide valuable insights into the metabolic state of cells. The main purpose of the study is to check the correlation between senescent tumor cells and enzymatic activity related with OXPHOS in colorectal cancer (CRC) cells.
MATERIALS AND METHODS
CRC sample preparation
We obtained CRC tissue samples from pathological specimens of CRC patients who underwent surgery at Ajou University Hospital. An experienced pathologist sampled fresh tissues from both normal and tumor areas immediately after resection. We included patients who had not received chemotherapy or radiotherapy before surgery. Cases that failed to undergo any one of the three staining procedures for senescence and mitochondrial OXPHOS were previously excluded. Out of the 35 CRC cases, 28 met the inclusion criteria.
Single-cell RNA-sequencing
Publicly available CRC single-cell RNA-sequencing (scRNAseq) dataset, GSE166555 and GSE178341 were used for analysis [19,20]. Among patients, patient #8 and #C103 who markedly expressed CDKN2A was used for analysis to minimize the batch effect, respectively. CDKN2A– and CDKN2A+ cancer cells were analyzed and compared for further analysis. Detailed information about preprocessing steps and cancer cell extraction is described in our previous study [8].
Senescence-associated β-galactosidase staining
For senescence-associated β-galactosidase (SA-β-Gal) staining, representative fresh CRC tissue samples were immediately frozen at –20°C using a cryostat (Leica Biosystems, Nußloch, Germany) after being coated with OCT compound (Sakura Finetek, Torrance, CA, USA). Fresh frozen sections with a thickness of 6 μm were placed on slides. After fixation with phosphate buffer saline (PBS) buffer containing formaldehyde for 1 minute, the sections were incubated in a mixture containing 1 mg/mL of 5-bromo-4-chloro-3-indolyl-β-D-galactopyranoside (X-gal) (Bioneer, Daejeon, Korea), 40 mM citric acid/sodium phosphate buffer at pH 5.8, 5 mM potassium ferrocyanide, 5 mM potassium ferricyanide, 150 mM NaCl, and 2 mM MgCl2 at 37°C in an incubator for 14 hours. Following incubation, the slides were washed three times with PBS. Counterstaining was carried out with nuclear fast red solution at room temperature for approximately 5 minutes. After the counterstaining was complete, the slides were mounted using mounting solution (Sakura Finetek). All slides with SA-β-Gal staining were scanned at a minimum of 20× magnification using an Aperio AT2 slide scanner (Leica Biosystems).
NADH-TR staining
For NADH-TR staining, fresh frozen section slides were incubated in a mixture of nitro tetrazolium blue chloride solution (2 mg/mL NBT/0.05 M Tris buffer pH 7.6; Sigma-Aldrich, St. Louis, MO, USA) and reduced nicotinamide adenine dinucleotide solution (1.6 mg/mL; NADH/0.05 M Tris buffer pH 7.6; Sigma-Aldrich). The sections were then incubated at 37°C in an incubator for 2 hours. After incubation, the slides were washed three times with PBS and counterstained with nuclear fast red solution at room temperature for approximately 5 minutes. Once the counterstaining was complete, the slides were mounted using mounting solution (Sakura Finetek). All slides with NADH staining were scanned at a minimum of 20× magnification using an Aperio AT2 slide scanner (Leica Biosystems).
SDH staining
For SDH staining, fresh frozen section slides were incubated in a mixture containing 0.2 M phosphate buffer at pH 7.4, 0.1 M MgCl2, 0.2 M succinic acid (Sigma-Aldrich), and 2.4 mM nitro tetrazolium blue chloride (Sigma-Aldrich). The sections were incubated at 37°C in an incubator for 2 hours. After incubation, the slides were washed three times with PBS and counterstained with nuclear fast red solution at room temperature for approximately 5 minutes. Once the counterstaining was complete, the slides were mounted using mounting solution (Sakura Finetek). All slides with SDH staining were scanned at a minimum of 20× magnification using an Aperio AT2 slide scanner (Leica Biosystems).
Evaluation of OXPHOS in CRC tissues
To investigate senescent cells in CRC, we analyzed virtual slides of SA-β-Gal staining in CRC tissue sections. SA-β-Gal positivity was identified in cancer cells, stromal cells, and/or macrophages. We specifically assessed the percentages of SA-β-Gal positive cancer cell areas relative to the total cancer cell areas. For the analysis of OXPHOS in CRC, we examined virtual slides of NADH-TR and SDH staining in CRC tissue sections. We compared the intensity of purple-blue formazan pigments in cancer cells to that in normal epithelial cells or stromal cells and considered it positive if the intensity was similar or more intense. We also evaluated the percentages of positive cancer cell areas within the total cancer cell areas. Furthermore, we compared the patterns of positivity in SA-β-Gal, NADH-TR, and SDH staining in CRC tissue sections. We conducted further analysis to determine the association between the positivity of SA-β-Gal, NADH-TR, and SDH staining and clinicopathologic parameters of CRC patients.
Cell culture and hydrogen peroxide treatment
CRC cell lines SNU254 and SNU1544 were purchased from the Korean Cell Line Bank (KCLB, Seoul, Korea). The cells were maintained in RPMI supplemented with 10% fetal bovine serum. For the hydrogen peroxide treatment experiment, cells were treated with 100 μM and 200 μM H2O2/dimethyl sulfoxide solution. After treatment for the indicated times, cells were incubated at 37°C for 2 days. The cells were then washed with PBS buffer and subjected to measurement of hydrogen peroxide using live cell imaging microscopy.
3-(4,5-Dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) assay
Control and hydrogen peroxide-treated SNU254 and SNU1544 cells were preincubated at a concentration of 1 × 106 cells/mL in culture medium with 1 μg/mL actinomycin C1 for 3 hours at 37°C and 5% CO2. Cells were seeded at a concentration of 5 × 104 cells/well in 100 μL culture medium containing 1 μg/mL actinomycin C1 into microplates. Cell cultures were incubated for 24 hours at 37°C with 5% CO2. After the incubation period, MTT labeling reagent (final concentration 0.5 mg/mL) was added to each well. The microplate was incubated for 4 hours in a humidified atmosphere. Solubilization solution was added into each well followed by overnight incubation in the incubator in a humidified atmosphere. Complete solubilization of the purple formazan crystals was checked and the absorbance of the samples was measured using a microplate reader. The wavelength to measure absorbance of the formazan product was 600 nm.
Statistical analysis
The statistical analysis was conducted using IBM SPSS ver. 25.0 (IBM Corp., Armonk, NY, USA). The Student independent t test or Fisher’s exact tests were used to analyze the parameters. A linear regression was performed in Prism 9 software (GraphPad Software, San Diego, CA, USA). A p-value of less than or equal to .05 was considered statistically significant. For adjusted p-values in gene set enrichment analysis (GSEA), the raw p-values were ranked in ascending order. The adjusted p-values are calculated for each gene set using the ranked raw p-values and the formula of the Benjamin-Hochberg procedure.
RESULTS
Metabolic features of senescent tumor cells
To elucidate the metabolic features of senescent tumor cells, we employed scRNA-seq using the publicly available CRC dataset, GSE166555 [19]. To minimize the batch effect, we exclusively analyzed patients who exhibited marked expression of p16INK4A (CDKN2A), as described in our previous study. Among the 274 CDX2 and EPCAM double-positive cancer cells, 198 cells tested positive for CDKN2A, while 76 cells tested negative (Fig. 1A). To validate the senescence-related phenotype of CDKN2A-positive CRC cells, we conducted GSEA. GSEA clearly demonstrated that CDKN2A-positive cells exhibit senescent phenotypes (Fig. 1B). Numerous prior studies have already reported that mitochondrial dysfunction is a hallmark of senescence in primary cells from normal tissues [21-23]. However, it remains unknown in the context of senescent tumor cells. Similar to senescent cells from normal tissues, senescent tumor cells exhibit OXPHOS dysfunction (Fig. 1C). GSEA also revealed a significant upregulation of the electron transport from NADH to ubiquinone in senescent tumor cells (Fig. 1D). However, the electron transport from ubiquinol to cytochrome c was inhibited in senescent tumor cells (Fig. 1E). Furthermore, ATP synthesis was downregulated in senescent tumor cells (Fig. 1F). These findings were validated in another scRNA-seq dataset, GSE178341 (Supplementary Fig. S1) [20]. Therefore, it is highly suggested that the upregulation of complexes 1 and 2 of OXPHOS in senescent tumor cells is a result of a potential feedback loop caused by dysfunctions in the subsequent electron transport pathways (Fig. 1G).
Single-cell RNA-sequencing (scRNA-seq) of human colorectal cancer tissues. (A) Uniform Manifold Approximation and Projection (UMAP) of public scRNA-seq dataset GSE166555 patient #8 who markedly expressed CDKN2A is shown. A total of 274 cells were analyzed. (B) Gene set enrichment analysis (GSEA) of “Fridman Senescence Up” (left panel) and “GOBP Cellular Senescence” gene sets was performed in CDKN2A+ vs. CDKN2A– cancer cells. (C) GSEA of a “GOBP Oxidative Phosphorylation” gene set was performed in CDKN2A+ vs. CDKN2A– cancer cells. (D) GSEA of a “GOBP Mitochondrial Electron Transport NADH to Ubiquinone” gene set was performed in CDKN2A+ vs. CDKN2A– cancer cells. (E) GSEA of a “GOBP Mitochondrial Electron Transport Ubiquinol to Cytochrome c” gene set was performed in CDKN2A+ vs. CDKN2A– cancer cells. (F) GSEA of a “GOBP ATP Synthesis Coupled Electron Transport” gene set was performed in CDKN2A+ vs. CDKN2A– cancer cells. (G) Schematic image of expected electron transport efficiency in senescent and non-senescent tumor cells. respectively. The thickness of red arrows indicates the efficiency of electron flow. ‘padj’ indicates the adjusted p-value. CRC, colorectal cancer; ES, enrichment score; NES, normalized enrichment score.
Patient characteristics according to the SA-β-Gal positivity
To confirm the metabolic phenotype suggested by the scRNAseq results, we performed enzymatic histochemistry analysis using NADH-TR and SDH. NADH is a component of complex 1 of the electron transport chain, which is part of OXPHOS [24]. In contrast, SDH is a component of complex 2 of OXPHOS [24]. However, there was a challenge in performing p16INK4A immunostaining because it is typically performed on formalin-fixed paraffin-embedded (FFPE) CRC samples, while enzymatic histochemistry can only be conducted on fresh or frozen tissues. To address this issue, we employed SA-β-Gal staining, another standard senescence marker that can be utilized in frozen tissues. Previously, we had confirmed the correlation between p16INK4A and SA-β-Gal staining in CRC cells [7,8]. Therefore, following analysis was performed based on the SA-β-Gal staining rather than immunostaining. A wide variation in the number of SA-β-Gal– positive cells were observed among CRC samples. These cells were distributed not only in the cancerous epithelium but occasionally in the normal epithelium (Supplementary Fig. S2), but also in stromal fibroblasts, and macrophages [2,7]. In this study, SA-β-Gal–positive malignant epithelial cells were exclusively examined to analyze the traits of senescent tumor cells in each CRC sample. The proportion of SA-β-Gal staining ranged from less than 1% to more than 90% across the CRC cases. We divided the CRC samples into two groups based on the proportion of SA-β-Gal–positive cancer cells: a low SA-β-Gal group (less than 30%) and a high SA-β-Gal group (more than 30%). The clinicopathologic characteristics of patients according to the SA-β-Gal group were shown in Table 1. We did not find any significant associations between clinicopathologic characteristics and the low or high SA-β-Gal–positive groups.
Clinicopathologic characteristics of CRC according to the SA-β-Gal staining intensity
NADH-TR and SDH staining in senescent tumor cells in vivo
To assess the OXPHOS activity in senescent tumor cells, we conducted NADH-TR staining (Supplementary Fig. S3) and SDH staining (Supplementary Fig. S4) on consecutive fresh frozen CRC tissue sections. Additionally, we performed SA-β-Gal staining and enzymatic histochemical analysis on serial sections of frozen tissues. Our findings revealed a consistent correlation between SA-β-Gal staining and both NADH-TR and SDH staining. The intensity of NADH-TR and SDH staining was found to be higher in cancer tissues with high SA-β-Gal activity (Fig. 2). Conversely, cancer tissues with low SA-β-Gal activity exhibited lower NADH-TR and SDH staining intensity (Fig. 3). A linear regression analysis clearly showed that there is a positive correlation between SA-β-Gal and NADH-TR/SDH staining (Fig. 4), suggesting that complex 1 and 2 activities of OXPHOS in senescent tumor cells are upregulated despite of the mitochondrial dysfunctions.
Nicotinamide adenine dinucleotide dehydrogenase–tetrazolium reductase (NADH-TR) and succinate dehydrogenase (SDH) staining results in the senescence-associated β-galactosidase (SA-β-Gal) high patient. The staining results of SA-β-Gal (A), NADH-TR (B), and SDH staining (C) are shown in the SA-β-Gal high patient.
Nicotinamide adenine dinucleotide dehydrogenase–tetrazolium reductase (NADH-TR) and succinate dehydrogenase (SDH) staining results in the senescence-associated β-galactosidase (SA-β-Gal) low patient. The staining results of SA-β-Gal (A), NADH-TR (B), and SDH staining (C) are shown in the SA-β-Gal low patient.
Scatter plots according to the nicotinamide adenine dinucleotide dehydrogenase–tetrazolium reductase (NADH-TR) (A) and succinate dehydrogenase (SDH) (B) staining results compared with senescence-associated β-galactosidase (SA-β-Gal) staining. p-values are obtained using the hypothesis test for regression slope.
MTT assay in control and senescent tumor cells in vitro
MTT is a widely used method for assessing the viability of cultured cells [25]. Additionally, MTT can be utilized to evaluate mitochondrial function. MTT is reduced to formazan, which strongly absorbs light in the visible range and has low solubility in water [26]. Notably, MTT can be directly reduced by complex I of the respiratory chain, providing a measure of complex I activity in the mitochondria [26]. To assess the complex I activity of senescent tumor cells in vitro, senescence in CRC cell lines SNU1544 and SNU254 was induced by hydrogen peroxide treatment, confirmed by SA-β-Gal staining (Supplementary Fig. S5A). The mRNA levels of representative SASPs were assessed, and SASPs were found to be increased in hydrogen peroxidetreated cancer cells (Supplementary Fig. S5B). MTT assays revealed an upregulation of complex I activity in senescent tumor cells in vitro (Supplementary Fig. S5C), which is consistent with in vivo findings.
DISCUSSION
Cancer cells, owing to their rapid proliferation rate and heightened metabolic activity, employ distinctive energy production mechanisms to meet their energy requirements [27]. Warburg originally observed the phenomenon of markedly increased glucose uptake in cancer cells, now widely recognized as the ‘Warburg effect’ [28,29]. Although Warburg suggested that elevated glycolysis in cancer cells resulted from mitochondrial defects [29], recent studies have shown that the mitochondrial pathway in cancer cells may still be partially functional and contribute to energy production [30]. Consequently, mitochondrial metabolism remains an attractive target for cancer therapy, with numerous investigations exploring this avenue [31,32]. Nevertheless, a notable challenge in the development of agents targeting mitochondrial metabolism is the resemblance between the metabolic processes of cancer cells and cytotoxic T cells [32-34]. The use of these agents to target cancer cells may inadvertently lead to the elimination of cancer-fighting cytotoxic T cells, potentially exacerbating a prognosis [30,32]. Therefore, the pursuit of metabolism-targeting agents specific to cancer cells represents an emerging frontier in cancer treatment.
Research into cancer metabolism has made significant progress in understanding non-senescent tumor cells [30]. However, the metabolic characteristics of senescent tumor cells remain less clear. Cellular senescence is widely recognized to be closely associated with mitochondrial dysfunction [35]. Although there is ongoing debate about whether cellular senescence induces mitochondrial dysfunction [36] or vice versa [37], it is evident that there is a positive correlation between cellular senescence and mitochondrial dysfunction [38]. Indeed, mitochondrial dysfunction has been observed in various cellular senescence models. Nelson and colleagues demonstrated a link between replicative senescence and mitochondrial dysfunction [39]. Moiseeva and colleagues reported a correlation between OIS and mitochondrial dysfunction [23]. Furthermore, Yoon et al. [21] revealed that the induction of senescence through iron chelation resulted in a reduction in complex 2 activity, preceding the onset of p27kip1-mediated cell cycle arrest in both Chang and normal liver cells. Likewise, senescence induced by transforming growth factor β1 in Mv1Lu lung epithelial cells involved the inhibition of complex 4, leading to the generation of mitochondrial reactive oxygen species [22]. To the best of our knowledge, this is the first study to analyze the metabolic features of senescent cancer cells in CRC. In our results, senescent tumor cells exhibited relatively more significant defects in OXPHOS compared to non-senescent tumor cells, as indicated by impaired electron transport to cytochrome c and reduced ATP production. This is the same line with the various studies based on multiple senescence models using normal cells [21-23]. However, intriguingly, there was a notable upregulation in the activity of electron transport from NADH to ubiquinone, a function typically associated with complex 1 [24]. While further direct evidence is required, these data suggest that electron transport from ubiquinone to complex 3 may be compromised in senescent tumor cells, leading to an upregulation in the gene expression of complex 1 and 2 as part of a potential feedback loop.
NADH-TR and SDH staining are enzyme histochemical methods that are widely used to examine the oxidative activity of cells, particularly in muscle biopsy. These stains can reveal the oxidative capacity of muscle fibers and their architectural changes in some muscle disorders. NADH and SDH are components of complex 1 and 2 of the electron transport chain, respectively. By visualizing the activity of these enzymes, NADH-TR and SDH staining can provide valuable insights into the metabolic state of cells [40]. In addition to diagnosing muscle disorders, these methods can be applied to screen the OXPHOS status of cancer cells as well as senescent cells. OXPHOS is an important metabolic process not only in cells that drive cancer drug resistance [41], but also in senescent cells associated with the aging process [9]. Therefore, NADH-TR and SDH staining can be used to identify the metabolic state of these cells and provide insights into their response to anticancer therapy or anti-aging therapy.
While cellular senescence in cancer cells has traditionally been considered a defense mechanism against cancer progression [42], recent studies have revealed that senescent tumor cells can actually promote cancer progression by influencing the tumor microenvironment [5-7]. In fact, senescent tumor cells themselves exhibit a higher invasion ability due to the secretion of various extracellular matrix-modulating proteins, such as matrix metalloproteinases [2]. As a result, senescent tumor cells are often located at the invasive front of the cancer tissue rather than its center [5,7]. Moreover, the reversal of the halted proliferation of senescent tumor cells may be associated with cancer relapse in distant metastatic organs. Therefore, the detection of senescent tumor cells is crucial for predicting patient prognosis and potential cancer relapse in the near future. However, there are limited methods available for the detection of senescent tumor cells. Furthermore, many markers for cellular senescence are tumor suppressor genes, such as p53, p21Waf1, and p16INK4A [2,43]. Consequently, they have a higher mutation rate, which can result in the overexpression of malfunctioning proteins due to potential feedback loops [2,44,45]. To overcome this challenge, additional senescence markers, such as a decrease in lamin B1 or detection of H3K9-trimethylation, are utilized [46]. In our experience, the staining of multiple senescence markers using serial sections, combined with CDK inhibitors, and accessory senescence markers are essential to minimize false-positive results when detecting senescent tumor cells in FFPE tissues [2]. While SA-β-Gal is considered a reliable marker for detecting senescent tumor cells in fresh or frozen CRC tissues [47], it often produces false positives due to its inherent β-galactosidase activity [48]. In our study, we found that enzymatic histochemistry using NADH-TR and SDH staining can effectively identify senescent tumor cells in CRC. Therefore, utilizing NADH-TR and SDH staining in conjunction with SA-β-Gal staining, using serial sections, represents a valuable approach to reduce the risk of false positivity associated with SA-β-Gal staining alone.
There are several limitations in this study. Firstly, the results mainly rely on scRNA-seq, enzyme histochemistry analysis, and in vitro MTT assay, which provide relatively indirect evidence compared to various mitochondrial supercomplex assays [39]. Therefore, to confirm the specific dysfunction of the mitochondrial electron transport system, various assays from complex 1 to 4 should be conducted, which will be part of our future studies. Enzyme histochemistry along with functional assays for electron transport complexes may address overall efficiency of the electron transport and OXPHOS system. However, despite similar results in different scRNA-seq datasets in our study consistently indicating that complex 1 and 2 are functionally upregulated in senescent tumor cells [19,20], the reasons for this phenomenon remain unknown. Secondly, the number of patients analyzed in this study was relatively small. Therefore, a larger number of patients will be needed to strengthen the results.
Supplementary Information
The Data Supplement is available with this article at https://doi.org/10.4132/jptm.2023.10.09.
Single-cell RNA-sequencing (scRNA-seq) of human colorectal cancer (CRC) tissue, GSE178341. UMAP of Public scRNA-seq dataset GSE178341 Patient #C103 who markedly expressed CDKN2A is shown. Total 1,309 cells were analyzed (left upper panel). GSEA of “Fridman Senescence Up” was performed in CDKN2A+ vs. CDKN2A- cancer cells (right upper panel). GSEA of a “GOBP Mitochondrial Electron Transport NADH to Ubiquinone”, “GOBP Mitochondrial Electron Transport Ubiquinol to Cytochrome c”, “GOBP ATP Synthesis Coupled Electron Transport” gene sets were performed in CDKN2A+ vs. CDKN2A- cancer cells (lower panels). ‘padj’ indicates the adjusted p-value.
Senescence-associated β-galactosidase (SA-b-Gal) staining staining results in human colorectal cancer tissues. (A) The staining result of SA-b-Gal–low patient. (B) The staining result of SA-b-Gal–high patient.
Nicotinamide adenine dinucleotide dehydrogenase–tetrazolium reductase (NADH-TR) staining results in human colorectal cancer tissues. (A) The staining result of NADH-TR–low patient. (B) The staining result of LADH-TR–high patient.
Succinate dehydrogenase (SDH) staining results in human colorectal cancer tissues. (A) The staining result of SDH-low patient. (B) The staining result of SDH-high patient.
In vitro MTT assays in H2O2-induced senescence model in colorectal cancer (CRC) cell lines. (A) Senescenceassociated β-galactosidase (SA-b-Gal) staining in control and H2O2-treated CRC cell line, SNU1544 (left panel). Quantification data for SA-b-Gal staining of SNU1544 and SNU254 is shown in the right panel. (B) The mRNA expression of representative senescence-associated secretory phenotypes in control and H2O2-treated CRC cell lines. (C) Relative absorbance in 600 nm of blank, control, and H2O2-treated CRC cells in MTT assay.
Notes
Ethics Statement
This study was approved by the Institutional Review Board of Ajou University Hospital (AJIRB-BMR-OBS-16-218). All patients provided informed consent.
Availability of Data and Material
The datasets generated or analyzed during the current study are available in the Gene expression Omnibus (GEO) repository from National Center for Biotechnology Information (NCBI) with accesion number GSE166555 and GSE178341. Persistent weblinks are https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE166555 and https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE178341.
Code Availability
Not applicable.
Author contributions
Conceptualization: JSS, JHK. Data curation: SSP, SHP, JHK. Formal analysis: JSS, SSP, YHK, TGK, SYE, DHL, SHP, TJP, JHK. Funding acquisition: SSP, JHK. Investigation: JSS, SSP, JHK. Methodology: JSS, SSP, YHK, SHP, TGK, JHK. Project administration: JHK. Resources: JSS, JHK. Software: SSP. Supervision: TJP, JHK. Validation: SSP, JHK. Visualization: SSP, SHP, JHK. Writing—original draft: JSS, TGK, SSP, JHK. Approval of final manuscript: all authors.
Conflicts of Interest
J.H.K., a contributing editor of the Journal of Pathology and Translational Medicine, was not involved in the editorial evaluation or decision to publish this article. All remaining authors have declared no conflicts of interest.
Funding Statement
National Research Foundation of Korea, NRF-2016R1D1A1B02010452 (J-HK).
Korea Health Industry Development Institute (KHIDI), a grant of the MD-PhD/Medical Scientist Training Program; Funding number is not applicable (SSP).
Ajou University Hospital, a grant of the MD-PhD/Medical Scientist Training Program; Funding number is not applicable (SSP).