CPT1A mediates radiation sensitivity in colorectal cancer
eLife Assessment
This study reports a valuable finding for the treatment of colorectal cancer (CRC), as the authors demonstrated that the enzyme CPT1A plays an significant role in the response to radiotherapy in CRC patients. However, the reviewers found that the results presented are still incomplete.
https://doi.org/10.7554/eLife.97827.3.sa0Valuable: Findings that have theoretical or practical implications for a subfield
- Landmark
- Fundamental
- Important
- Valuable
- Useful
Incomplete: Main claims are only partially supported
- Exceptional
- Compelling
- Convincing
- Solid
- Incomplete
- Inadequate
During the peer-review process the editor and reviewers write an eLife Assessment that summarises the significance of the findings reported in the article (on a scale ranging from landmark to useful) and the strength of the evidence (on a scale ranging from exceptional to inadequate). Learn more about eLife Assessments
Abstract
The prevalence and mortality rates of colorectal cancer (CRC) are increasing worldwide. Radiation resistance hinders radiotherapy, a standard treatment for advanced CRC, leading to local recurrence and metastasis. Elucidating the molecular mechanisms underlying radioresistance in CRC is critical to enhance therapeutic efficacy and patient outcomes. Bioinformatic analysis and tumour tissue examination were conducted to investigate the CPT1A mRNA and protein levels in CRC and their correlation with radiotherapy efficacy. Furthermore, lentiviral overexpression and CRISPR/Cas9 lentiviral vectors, along with in vitro and in vivo radiation experiments, were used to explore the effect of CPT1A on radiosensitivity. Additionally, transcriptomic sequencing, molecular biology experiments, and bioinformatic analyses were employed to elucidate the molecular mechanisms by which CPT1A regulates radiosensitivity. CPT1A was significantly downregulated in CRC and negatively correlated with responsiveness to neoadjuvant radiotherapy. Functional studies suggested that CPT1A mediates radiosensitivity, influencing reactive oxygen species (ROS) scavenging and DNA damage response. Transcriptomic and molecular analyses highlighted the involvement of the peroxisomal pathway. Mechanistic exploration revealed that CPT1A downregulates the FOXM1-SOD1/SOD2/CAT axis, moderating cellular ROS levels after irradiation and enhancing radiosensitivity. CPT1A downregulation contributes to radioresistance in CRC by augmenting the FOXM1-mediated antioxidant response. Thus, CPT1A is a potential biomarker of radiosensitivity and a novel target for overcoming radioresistance, offering a future direction to enhance CRC radiotherapy.
Introduction
Colorectal cancer (CRC) is the second-highest cause of cancer-related mortality (Siegel et al., 2023). Radiotherapy is crucial for CRC management, especially in patients with locally advanced rectal cancer (cT3-4N+) (Glynne-Jones et al., 2017). Neoadjuvant therapies show clinical or pathological complete response in 16–30% of patients, realising downstaging in approximately 60% of patients, significantly enhancing local control, and facilitating curative surgery (Cercek et al., 2018). Moreover, radiotherapy is beneficial in initially unresectable and recurrent cases with limited metastasis to organs, such as the liver and lungs (Cervantes et al., 2023). Nonetheless, the effectiveness of radiotherapy is affected by radioresistance, which precipitates tumour relapse and metastasis and currently lacks an efficacious clinical resolution. Unlocking the molecular mechanisms underlying CRC radioresistance will enhance outcomes and improve patient prognoses.
Reactive oxygen species (ROS) are byproducts of normal cellular metabolism occurring in organelles, such as mitochondria, endoplasmic reticulum, and peroxisomes (Bhattacharyya et al., 2014). Within mitochondria, 90% of cellular ROS are generated by complexes I, II, and IV of the electron transport chain (Glasauer and Chandel, 2014). Exogenous stimuli, including radiation, can cause significant, acute elevations in ROS levels (Hecht et al., 2024). ROS function ambivalently in the intracellular signalling and redox homoeostasis of tumour cells (Shah et al., 2024). ROS amplify oncogenic phenotypes, such as proliferation and differentiation, hasten the accumulation of metastasis-inducing mutations, and aid tumour cell survival under hypoxic conditions (Palma et al., 2024). However, excess ROS precipitate apoptosis and other cell death types from oxidative stress (Palma et al., 2024). Thus, intracellular ROS generation is meticulously monitored and regulated by a comprehensive ROS-scavenging system encompassing antioxidants and antioxidative enzymes (Palma et al., 2024; Shah et al., 2024).
Radiation ionises water molecules, creating an intracellular surge of ROS, which indirectly cause two-thirds of DNA damage (Chio and Tuveson, 2017). Consequently, ROS scavenging inevitably influences cancer cell radiosensitivity (Skvortsova et al., 2015). Increased expression and activity of antioxidative enzymes, such as peroxidase (POD), catalase (CAT), glutathione peroxidase, and glutathione reductase, are correlated with radiosensitivity (Hecht et al., 2024).
Carnitine palmitoyltransferase 1 (CPT1) is an outer mitochondrial membrane that catalyses the rate-limiting step of fatty acid oxidation and is absent in several tumours (Melone et al., 2018). The CPT1 family contains three isoforms: CPT1A, CPT1B, and CPT1C. Research on CPT1A has been detailed (Schlaepfer and Joshi, 2020). CPT1A is critical to cancer cell growth, survival, and drug resistance, making it an attractive target (Qu et al., 2016). CPT1A also interacts with other key pathways and factors regulating gene expression and apoptosis in cancer cell (Qu et al., 2016). However, its role in CRC and radiotherapy resistance is unclear.
Previously, we found that various metabolic pathways, including fatty acid metabolism, are closely related to tumour radioresistance. Therefore, this study mainly focused on CPT1A, which affects CRC radiosensitivity, to reveal novel therapeutic strategies to mitigate radiotherapy resistance and improve clinical outcomes.
Results
Low expression of CPT1A in CRC tumours
Analysis of CRC mRNA sequencing arrays from GEO consistently indicated low CPT1A mRNA levels in CRC (Figure 1A). The transcript levels of CPT1A were lower in most colon cancers (19/24 pairs) than in the adjacent non-tumour tissues (Figure 1B). Meanwhile, CPT1A protein levels were lower in most CRC tissues (14/16 pairs) than in the adjacent non-tumour tissues (Figure 1C). Further exploration of CPT1A expression in TCGA revealed that CPT1A mRNA levels were significantly lower in colon adenocarcinoma (COAD) tissues than in the adjacent non-tumour tissues (Figure 1D). Similar results were observed for rectal adenocarcinoma (READ) (Figure 1E). Immunohistochemistry (IHC) staining of the cancer-adjacent borders of two patients showed that CPT1A protein levels were lower in CRC tissues than in the nearby non-tumour tissues (Figure 1F). These findings provide comprehensive evidence supporting the downregulation of CPT1A expression at both mRNA and protein levels in CRC.
Low CPT1A CRC exhibits radioresistance and poor overall survival
Based on the IHC scores, patients were divided into CPT1A high- and low-expression groups (Figure 2A). Kaplan-Meier survival analysis indicated that low CPT1A-expressing patients had low overall survival (OS) (Figure 2B). Survival analysis of READ patients with low CPT1A using GEPIA also suggested a low OS; the difference was statistically close but not significant (p=0.061, Figure 2C). To further investigate the correlation between CPT1A and radiotherapy efficacy, IHC staining and scoring for CPT1A were performed on samples from 43 patients with rectal cancer who received neoadjuvant radiochemotherapy. Among patients with a tumour regression grade (TRG) score of 1, which indicates a minimal number of residual tumour cells, a high IHC score for CPT1A is often observed. Conversely, patients with a TRG score of 2, indicating more residual tumour cells, exhibited relatively lower IHC staining intensity and overall scores (Figure 2D). A comparison of groups with TRG 2–3 and TRG 1 revealed significantly higher IHC scores in the TRG 1 group (Figure 2E). The TRG and IHC scores for CPT1A showed a negative correlation (R=−0.3430 and p=0.024) (Figure 2F). In summary, low CPT1A expression was associated with poor OS, high TRG scores, and a high probability of radioresistance.
Decreased CPT1A expression contributes to radiation resistance in CRC cells
Through a colony formation assay (CFA) and multi-target single-hit survival model, we found that SW480, Caco-2, SW620, HT-29, and cells were more resistant to radiation than HCT15, RKO, and HCT 116 cells (Figure 3—figure supplement 1A and B, Supplementary file 1). Furthermore, the background expression levels of CPT1A in the above cells revealed that CPT1A transcription and protein levels were higher in radiation-resistant cells than in radiosensitive cells (Figure 3—figure supplement 1C and D).
Accordingly, we constructed stable CRC cell lines with CPT1A knockout/overexpression. We transfected the CRISPR/Cas9 lentivirus into HCT 116 cells (highest CPT1A expression and radiosensitive) and used western blotting to verify that the knockout efficiency of the second site was the highest, while mRNA levels were significantly reduced (Figure 3A and B). We used the second knockout site for subsequent in vitro and in vivo experiments. We also transfected the CPT1A-overexpressing lentivirus into SW480 cells (lowest CPT1A expression and radioresistant) and verified that CPT1A protein and mRNA levels increased (Figure 3A and B).
The CFA and multi-target single-hit survival model suggested that radioresistance increased with CPT1A knockout (D0=1.526 vs 1.993, p<0.05, Figure 3C and D), whereas the radioresistance of cells decreased with CPT1A overexpression (D0=2.724 vs 1.963, p<0.01, Figure 3E and F). The comet assay suggested that the proportion of DNA in the tail of cells decreased with CPT1A knockout, indicating improved damage repair. With CPT1A overexpression, the proportion of DNA in the tail increased, indicating reduced damage repair capabilities (Figure 3G). We also detected γ-H2A.X expression at different time points after 6 Gy irradiation; with CPT1A knockout, γ-H2A.X disappeared from cells faster, indicating improved cell damage repair (Figure 3H). With CPT1A overexpression, the disappearance rate of γ-H2A.X in cells was slower, remaining detectable even after 24 hr, indicating diminished cellular repair capability (Figure 3I). Collectively, our data suggest that CPT1A radiosensitises intrinsically radioresistant cells.
We generated radioresistant cell lines (HCT-15-25F and HCT-15-5F) from the HCT-15 parent line by fractionated irradiation (Figure 3—figure supplement 2A). Using CFAs and the multi-target single-hit survival model, we observed increased radioresistance in these new cell lines, as indicated by higher D0 values than those of the parental cells (HCT-15 D0=2.957, HCT-15-5F D0=3.240, HCT-15-25F D0=3.822) (Figure 3—figure supplement 2B and C). We also found significantly decreased CPT1A protein expression in HCT-15-25F and HCT-15-5F cells compared with that in the original cells (Figure 3—figure supplement 2D). Thus, HCT-15-25F cells were selected for further analysis. Transfection with the CPT1A-overexpressing lentivirus led to stable overexpression in HCT-15-25F cells (Figure 3—figure supplement 2E). To assess the impact of CPT1A overexpression on DNA repair capacity, we monitored γ-H2A.X expression at different time points following 6 Gy irradiation. CPT1A overexpression resulted in a slow disappearance of intracellular γ-H2A.X, suggesting enhanced DNA damage repair capability (Figure 3—figure supplement 2F). CFAs and the multi-target single-hit survival model revealed that CPT1A overexpression increased the radiosensitivity of radioresistant cells (D0=2.871 vs 2.581, p<0.05; Figure 3—figure supplement 2G and H), indicating that CPT1A exerts a radiosensitising effect in inducible radioresistant cell lines.
Diminished CPT1A expression in vivo also leads to tumour radioresistance
To further validate the impact of CPT1A on radiation resistance in vivo, we established a xenograft model using HCT116-NC and HCT116-KO cell lines in nude mice (Figure 4A). CPT1A-stabilising knockout increased tumour weights in mice, which persisted after radiotherapy (Figure 4B), suggesting that CPT1A knockout promotes tumour growth and confers increased resistance to radiation. IHC staining demonstrated a significant increase in Ki-67 staining intensity and the percentage of positive cells in CPT1A knockout tumours, indicating enhanced proliferative capacity that was further pronounced after radiotherapy (Figure 4C and D). In the absence of CPT1A knockdown, radiotherapy reduced the percentage of Ki67-positive cells in the xenograft tumours by 32.9% (approximately 39.6% of the pre-irradiation baseline). In contrast, upon CPT1A knockdown, radiotherapy only led to a 14.5% reduction in the percentage of Ki67-positive cells (approximately 15.6% of the pre-irradiation baseline; Figure 4D). Similarly, we performed xenograft model experiments using the SW480-RFP and SW480-OE cell lines (Figure 4E). CPT1A overexpression resulted in reduced tumour weight in mice, a trend that persisted even after radiotherapy (Figure 4F). Furthermore, as illustrated in Figure 4E and F, in the absence of CPT1A overexpression, radiotherapy resulted in a 0.10 g decrease in tumour weight (approximately 52.5% of the pre-irradiation weight), whereas with CPT1A overexpression, radiotherapy induced a more pronounced 0.12 g reduction in tumour weight (approximately 89.7% of the pre-irradiation weight). These findings suggested that CPT1A overexpression inhibits tumour growth and sensitises tumour cells to radiation. IHC staining for Ki-67 revealed significantly decreased staining intensity and percentage of positive cells in CPT1A-overexpressing tumours, indicating weakened proliferative capacity, which was further accentuated after radiotherapy (Figure 4G and H).
Low CPT1A levels accelerate post-radiation ROS scavenging
To study the mechanism of low CPT1A expression in radiotherapy resistance, we conducted differential gene expression analysis between HCT116-KO and HCT116-NC cells. The gene expression heatmap demonstrated high consistency among technical replicates for both HCT 116-NC and HCT 116-KO cells (Figure 5—figure supplement 1A). With CPT1A knockdown, we found 363 upregulated and 1290 downregulated genes (|log2(fold change)|>1 and q<0.05) (Figure 5—figure supplement 1B). We conducted Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis and Gene Ontology (GO) annotation for all differentially expressed genes (DEGs; Figure 5A; Figure 5—figure supplement 1C–E), showing that the peroxisomes, cell cycle nucleotide excision repair, and fatty acid degradation pathways were among the enriched KEGG pathways (Figure 5A). Peroxisomes are important organelles that maintain cellular redox balance by clearing ROS. Therefore, we targeted peroxisomal pathways. ROS levels in the cells were dynamically balanced, including ROS production and scavenging (Figure 5B). To investigate the effect of CPT1A on ROS, we examined them in stable CPT1A knockout/overexpression cells following 6 Gy irradiation and 1 hr of incubation with 2',7'-dichlorodihydrofluorescein diacetate (DCFH-DA), a fluorescent redox probe used to detect ROS in cells (Figure 5C). The total ROS levels in CPT1A knockout cells decreased, whereas those in CPT1A-overexpressing cells increased (Figure 5D). The main mechanism by which cells produce ROS under irradiation is through the X-ray ionisation of water molecules, which far exceeds those from oxidative phosphorylation and NOX enzymes (Figure 5B). Therefore, we speculated that the regulation of CPT1A by intracellular ROS levels may be attributed to increased ROS scavenging. We examined the GSH/GSSG ratio (Figure 5E and F) and SOD (Figure 5G and H), CAT (Figure 5I), and POD enzyme activities (results were unchanged, data not shown) in stable CPT1A knockout/overexpression cells; in CPT1A knockout cells, the GSH/GSSG ratio and SOD and CAT enzyme activities increased, and these changes were also observed under 6 Gy irradiation (Figure 5I). In contrast, in CPT1A-overexpressing cells, the GSH/GSSG ratio, and SOD and CAT activities decreased, and these changes were also observed after 6 Gy irradiation (Figure 5I).
We further validated the effect of CPT1A on ROS scavenging in radioresistant cells. Total ROS significantly reduced in radioresistant cells compared to those in HCT-15 control cells following 6 Gy irradiation (Figure 5—figure supplement 2A and C). Furthermore, CPT1A overexpression increased total ROS levels in radioresistant cells (Figure 5—figure supplement 2B and D), suggesting that CPT1A enhances ROS accumulation in radioresistant cells. Additionally, we assessed the GSH/GSSG ratio (Figure 5—figure supplement 2E and F) and SOD (Figure 5—figure supplement 2G and H), POD (Figure 5—figure supplement 2I and J), and CAT enzyme activities (Figure 5—figure supplement 2K) in radioresistant cells, both at baseline and after 6 Gy irradiation, revealing an increased GSH/GSSG ratio and elevated activities of SOD, POD, and CAT enzymes in radioresistant cells compared to those in parental HCT-15 control cells. Furthermore, these changes persisted under 6 Gy irradiation.
To explore the influence of CPT1A overexpression on radiation resistance, we examined the GSH/GSSG ratio (Figure 5—figure supplement 2E and F) and SOD (Figure 5—figure supplement 2G and H), POD (Figure 5—figure supplement 2I and J), and CAT enzyme activities (Figure 5—figure supplement 2K) in HCT-15-25F-OE cells compared to those in control cells, showing that CPT1A overexpression restored the GSH/GSSG ratio and increased SOD and CAT enzyme activities but failed to restore POD enzyme activity. Additionally, under 6 Gy irradiation, the GSH/GSSG ratio and SOD enzyme activity were restored in CPT1A-overexpressing cells (Figure 5—figure supplement 2E–K).
Reduced CPT1A expression mediates radioresistance in CRC through increased expression of ROS-scavenging genes, facilitated by FOXM1
The reasons underlying the changes in SOD and CAT enzyme activities require further investigation. We examined the transcriptional and protein levels of SOD (SOD1, SOD2, and SOD3) and CAT (Figure 6A and B) and found that CPT1A knockout in cells increased both the mRNA and protein levels of SOD1, SOD2, and CAT, whereas the overexpression of CPT1A decreased them (Figure 6A and B). However, the function of CPT1A as a transcription factor is hitherto unreported, suggesting a potential regulatory mechanism mediated by other transcription factors. Thus, we employed a bioinformatics analysis by intersecting predicted or reported transcription factors known to regulate SOD1, SOD2, and CAT with upregulated DEGs. The Venn diagram highlights three transcription factors that met the criteria: FOXM1, LMNB1, and SAP30 (Figure 6C). FOXM1 transcription and protein levels were significantly increased in CPT1A knockout cells but decreased when CPT1A was overexpressed (Figure 6A and B). Other transcription factors showed no significant changes in protein levels (data not shown).
Additionally, we explored the correlation between FOXM1 and the downstream enzymes SOD1, SOD2, and CAT in READ and COAD. The results (Figure 6—figure supplement 1A and B) indicate a positive correlation between FOXM1 and SOD1 as well as SOD2. Furthermore, we used hTFtarget and JASPAR to predict FOXM1 binding sites in the promoters of SOD1, SOD2, and CAT (Table 1). Finally, a rescue experiment was conducted by overexpressing FOXM1 in HCT 116-NC and HCT 116-KO cells, demonstrating that downstream SOD1, SOD2, and CAT protein levels were restored by FOXM1 overexpression (Figure 6D). In summary, the downregulation of CPT1A increases FOXM1 mRNA and protein levels in CRC, promoting the transcription and translation of SOD1, SOD2, and CAT, thereby accelerating the scavenging of ROS produced after radiation exposure and ultimately leading to radiation resistance in CRC cells (Figure 6E).
Discussion
We elucidated the role of CPT1A in CRC and the molecular mechanisms involved in mediating radiosensitivity. CPT1A is often downregulated in CRC, and low CPT1A expression can worsen OS and increase the probability of radiochemotherapy resistance. Low CPT1A expression increases FOXM1 activity, promoting the transcription and translation of downstream SOD1, SOD2, and CAT, thereby facilitating the scavenging of radiation-induced ROS. Our study establishes CPT1A as an effective biomarker to predict CRC prognosis and radiotherapy sensitivity and proposes the molecular mechanism by which it mediates radiosensitivity through the FOXM1-SOD1/SOD2/CAT axis.
CPT1A is crucial to CRC initiation and progression (Mazzarelli et al., 2007). However, its role in CRC remains unclear. Using a weighted gene co-expression network analysis to explore hub genes in CRC development, we found that CPT1A is expressed at low levels in CRC and acts as a central anticancer gene, exhibiting excellent prognostic value (Wang et al., 2023). Sinomenine improves colitis-associated cancer by upregulating CPT1A (Zhang et al., 2022). In contrast, high CPT1A expression is associated with malignancy in CRC, and its inhibition ameliorates malignant phenotypes. CPT1A-mediated fatty acid oxidation promotes CRC metastasis (Wang et al., 2018). DHP-B, a CPT1A inhibitor, disrupts CPT1A-VDAC1 interaction in the mitochondria, increasing mitochondrial permeability and reducing oxygen consumption and energy metabolism in CRC cells (Hu et al., 2023). We observed significant downregulation of CPT1A expression in CRC, and low CPT1A expression was associated with worse prognosis and greater radiochemotherapy resistance, contrary to previous reports. This discrepancy may be related to the inherent heterogeneity of tumour tissue and differences in tumour stage.
The PGC1α/CEBPB/CPT1A axis, which enhances lipid β-oxidation, increases ATP and NADPH levels and promotes cellular radiation resistance in nasopharyngeal carcinoma (Tan et al., 2018). Tan et al. also discovered an interaction between CPT1A and Rab14, which transports fatty acids into the mitochondria, thereby facilitating lipid oxidation and cell survival under irradiation (Du et al., 2019). The HER1/2-MEK-ERK1/2-CPT1A/CPT2 axis reportedly enhances cell proliferation and confers radiation resistance in breast cancer (Han et al., 2019). However, no studies have investigated the association between CPT1A and radiosensitivity in CRC. Our research revealed that CPT1A is a radiation sensitivity gene, contradicting previous literature, possibly due to differences in cancer types.
Transcriptomic sequencing revealed that CPT1A regulates multiple pathways, including the peroxisomal pathway, which is responsible for ROS scavenging. CPT1A knockdown upregulated FOXM1, which, in turn, stimulated the transcription and translation of crucial antioxidative enzymes, SOD1, SOD2, and CAT, thereby expediting ROS clearance. This contradicts previous findings that CPT1A overexpression accelerates ROS production by increasing fatty acid β-oxidation, thereby promoting ageing phenotypes or augmenting cancer cells’ oxidative defences (Jiang et al., 2022; Joshi et al., 2020; Luo et al., 2021). Our findings diverge from those of other studies for two main reasons: first, previous research into the effects of CPT1A on ROS largely centred around mitochondria rather than the peroxisome; second, previous studies did not account for radiation, which significantly increases ROS production above cellular oxidative processes.
Forkhead box M1 (FOXM1) is a critical transcription factor for many cellular processes (Kalathil et al., 2020). Besides benefitting normal cell functions, it also regulates cancer processes, including growth, metastasis, and recurrence (Alimardan et al., 2023; Khan et al., 2023). FOXM1 also affects radiotherapy outcomes in many cancer types, including CRC (Kwon et al., 2021; Li et al., 2022; Liu et al., 2019; Pal et al., 2018; Takeshita et al., 2023; Xiu et al., 2018). In this study, we identified FOXM1 as a key regulator connecting CPT1A to ROS scavenging in CRC, exhibiting an inverse correlation with both CPT1A and ROS levels. FOXM1 is an essential transcription factor in intracellular redox, specifically in regulating the redox state of malignant mesothelioma cells (Cunniff et al., 2014). FOXM1-dependent fatty acid oxidation-mediated ROS modulation is a cell-intrinsic drug resistance mechanism in cancer stem cells (Choi et al., 2020). FOXM1 knockdown increases intracellular ROS levels and decreases the transcription levels of SOD2, CAT, PRDX, and GPX2 (Smirnov et al., 2016). Overall, our findings align with existing literature, highlighting the crucial role of FOXM1 in orchestrating the interplay between CPT1A and ROS homoeostasis. We further revealed that FOXM1 participates in the transcriptional regulation of SOD1.
Our study has some limitations. We only conducted radiosensitivity investigations of CPT1A in nude mice, which only demonstrated its regulatory role in CRC cell radiosensitivity. It remains unclear whether CPT1A can regulate the radiosensitivity of the entire tumour microenvironment in immunocompetent mice. The cell lines used in our study have different genetic backgrounds; specifically, HCT116 is microsatellite instable, whereas SW480 is not. Additionally, FOXM1 primarily localises to the nucleus, whereas CPT1A is a cytoplasmic protein; there is no known physiological basis for their co-localisation. Therefore, the specific mechanism through which CPT1A regulates FOXM1 expression requires further investigation. The reasons for the decreased expression of CPT1A in tumour cells remain unclear. Future studies should explore this in greater detail.
Materials and methods
Reagents and materials
Request a detailed protocolThe reagent suppliers are indicated in Key resources table.
Bioinformatic analyses
Request a detailed protocolRaw mRNA expression profiles and clinical features from GSE9348, GSE20916, GSE37364, GSE44076, GSE68468, and GSE110223 were downloaded from GEO (http://www.ncbi.nlm.nih.gov/geo/). The CPT1A mRNA expression in cancer and paired normal tissues was analysed using UALCAN (https://ualcan.path.uab.edu/) (Chandrashekar et al., 2022). Correlation analyses between genes were conducted using GEPIA2 (http://gepia2.cancer-pku.cn/) (Tang et al., 2019). Transcription factor prediction and promoter binding site analysis for SOD1, SOD2, and CAT were conducted using the hTF target database (http://bioinfo.life.hust.edu.cn/hTFtarget#) (Zhang et al., 2020). Additionally, JASPAR (http://jaspar.genereg.net/) was used to predict the binding sites of FOXM1 (Rauluseviciute et al., 2024).
Human
Human samples were collected from the Nanfang Hospital of Southern Medical University, Guangzhou, for CRC patients with or without radiochemotherapy with age over 18 and below 75 years. Two to four biopsies from each patient were collected to be used for qRT-PCR, western blotting, and sectioning. Informed consent forms were acquired from all the patients. Ethics approval for the use of human samples was obtained from Ethics Committee of Nanfang Hospital, Southern Medical University (No. NFEC-202304-K13).
RNA isolation and qRT-PCR
Request a detailed protocolA total of 48 rectal adenocarcinoma samples and their paired normal tissues were used for RNA isolation and qRT-PCR. Total RNA was extracted using TRIzol following the manufacturer’s instructions. cDNA was generated using the Evo M-MLV RT Mix Kit. mRNA expression was analysed using the SYBR Green Premix Pro Taq HS qPCR Kit on a QuantStudio6 Real-time PCR system, and β-actin was used for normalisation. Data were analysed using the 2-ΔΔCT method. The PCR primers are listed in Supplementary file 2.
Protein extraction and western blotting
Request a detailed protocolProteins of rectal adenocarcinoma and paired normal tissues were extracted from tissues and cells using a protein extraction kit, according to the manufacturer’s instructions. Briefly, NanoDrop was used to determine the protein concentration. The proteins were mixed in a 4:1 ratio with 5× loading buffer and denatured at 100°C in a water bath for 5 min. Subsequently, 20 μg of proteins were subjected to SDS-PAGE at 95 V. Subsequently, proteins were transferred onto PVDF membranes at a constant current. The membrane was blocked with 5% BSA and incubated with the corresponding primary antibody overnight at 4°C. Subsequently, the membrane was incubated with an HRP-conjugated secondary antibody (1:2000 dilution), and protein levels were detected using enhanced chemiluminescence.
Immunohistochemistry
Request a detailed protocolAll colon cancer (n=76) and rectal adenocarcinoma tissues (n=45) were collected before treatment and sectioned at 4 μm. Clinicopathological features of the patients were provided by the Department of Pathology (Table 2). Staining intensity was independently evaluated by two senior pathologists. The staining intensity was scored on a 4-point scale, where 0 represented no positive staining (negative), 1 represented light yellow (weak positive), 2 represented brownish yellow (positive), and 3 represented brown (strong positive) staining. The percentage of positive cells was also scored on a 4-point scale, where 1 point was assigned for ≤25% positive cells, 2 points for 26–50% positive cells, 3 points for 51–75% positive cells, and 4 points for >75% positive cells. The IHC score was calculated by multiplying the staining intensity score with the positive cell percentage score. Patients with an IHC score ≥6 were classified as the CPT1A-high group.
Cell culture and lentiviral infection
Request a detailed protocolThe CRC cell lines HCT-15, RKO, HCT 116, HT-29, Caco-2, SW480, and SW620 were purchased from ATCC (Manassas, VA, USA). All cell lines were authenticated with STR profiling, using cell line authentication services offered by ATCC, to avoid misidentification and tested negative for mycoplasma contamination. All cell lines were cultured in RPMI-1640 medium supplemented with 10% FBS. Cells were cultured at 37°C with 5% CO2. To establish radiation-resistant strains, HCT-15-25F and HCT-15-5F cells were generated using both conventional fractionated irradiation (2 Gy/fraction, 25 fractions, 5 fractions/week over 5 weeks) and large-fractionated irradiation (5 Gy/fraction, 5 fractions, for 1 week).
The full-length lentiviral expression of CPT1A with puromycin was constructed by VectorBuilder. CPT1A-targeting CRISPR/Cas9 lentiviral vectors (hCpt1a[gRNA#1]: AAATCTCTACTACACGGCCGATGTTACGACAGGTACCGTCCTT, hCpt1a[gRNA#2]: AGAAGGTAAGGACGGTACCTGTCGTAACATCGGCCGTGTAGTA, hCpt1a[gRNA#3]: CTGAACACTCCTGGGCAGATGCGCCGATCGTGGCCCACCTTTG) with RFP and Puro were constructed by VectorBuilder. Infection and in vitro transfection of cell lines were performed following the manufacturer’s protocol. The lentiviral full-length expression of FOXM1 with G418 was constructed using VectorBuilder.
CFA and multi-target single-hit survival model
Request a detailed protocolThe radiosensitivity of all human CRC cell lines was determined using a CFA with a multi-target single-hit model to the surviving fractions. Cells were plated in six-well plates and irradiated at doses of 0, 2, 4, 6, 8, and 10 Gy (6 MeV X-rays). The cells were then cultured for 10 days, stained with 1% crystal violet, and quantified using ImageJ version 1.8.0.
The surviving fraction for each dose was calculated using the following formula:
[(number of surviving colonies at dose X)/(number of cells seeded at dose X (average colonies arising from non-irradiated cells (0 Gy))/number of non-irradiated cells seeded)]. Survival curves were used to develop the multi-target single-hit model, SF = 1 − (1 − e−D/D0)×N, where SF is the surviving fraction, D is the radiation dose, and N is the extrapolation number. The multi-target single-hit model provides parameters related to radiation sensitivity and sensitisation ratios, thereby directly reflecting the cells’ radiation sensitivity.
Comet assay
Request a detailed protocolThe comet assay is a classical technique for assessing cellular radiation sensitivity. Essentially, cells that are more susceptible to radiation damage exhibit more DNA fragmentation post-irradiation, leading to a greater proportion of DNA in the tail (with less in the head). The comet assay was performed as previously described (Yu et al., 2022). Briefly, slides were covered with 100 μL of pre-warmed normal melting point agarose (2%) and placed on ice to solidify the first gel layer. Cells irradiated with 6 Gy were digested to obtain a cell suspension. Ten microlitres of the cell suspension were mixed with 80 μL of pre-warmed low melting point agarose (0.75%) and poured onto the slides. The slides were dipped in a cold lysis solution for 2 hr. After cell lysis, the slides were placed in a horizontal electrophoresis chamber filled with cold TAE solution, incubated for 25 min in the dark, and electrophoresed (1 V/cm). The slides were then neutralised in PBS for 5 min and stained with propidium iodide (PI) or DAPI. Comet images were captured using a fluorescence microscope. The percentage of DNA in the tail was analysed using the CASP 1.2.3 beta 1.
Animal model
Request a detailed protocolAll animal experiments were approved by the Institutional Animal Care and Use Committee of Nanfang Hospital (IACUC-PROJECT-20221128-003). Male BALB/c nude mice (5 weeks of age) were purchased from the Southern Medical University Laboratory Animal Center, China, and raised under specific pathogen-free conditions. All in vivo experiments were performed following institutional guidelines. To develop the xenograft tumour model, 5×106 cells were subcutaneously injected into the left flank of mice. On reaching a volume of 100 mm3, the tumours were irradiated twice at 8 Gy, and other parts of the mouse body were protected with a lead shield. The tumour volume was measured using Vernier callipers and calculated as 1/2×length×width×width. After the mice were euthanised with phenobarbital sodium, the tumours were excised, weighed, and embedded in paraffin for further experiments.
Transcriptomics
Request a detailed protocolHCT 116-NC and HCT 116-KO cells were collected and sent to RIBOBIO (Guangzhou, China) for polyA-seq transcriptome sequencing. RNA extraction, library preparation, and sequencing were performed according to the manufacturer’s instructions. To identify the DEGs between the two groups, the expression level of each transcript was calculated according to the transcript per million reads method. RSEM was used to quantify the gene abundance. Differential expression analysis was performed using DESeq2 (Liu et al., 2021). Genes with |log2(fold change)|>1 and q<0.05 were considered significantly differentially expressed. In addition, functional enrichment analysis, including using GO and KEGG, was performed to identify significantly enriched DEGs in GO terms and metabolic pathways at a Bonferroni-corrected p-value of 0.05 compared with the whole-transcriptome background. The heatmap and volcano of mRNA sequencing were conducted on ImageGP (https://www.bic.ac.cn/ImageGP) (Chen et al., 2022). GO functional enrichment and KEGG pathway analyses were performed using the ClusterProfiler package in R (Wu et al., 2021).
ROS detection
Request a detailed protocolROS were stained with a DCFH-DA probe and detected using flow cytometry. Briefly, the probe was diluted in cell culture media at a 1:1000 ratio to yield a final concentration of 10 μmol/L. The culture medium was replaced with 1 mL of DCFH-DA solution. The cells were incubated at 37°C for 20 min. After incubation, the cells were washed thrice with serum-free medium to remove any excess DCFH-DA.
One hour after 6 Gy irradiation, the cells were digested into single-cell suspensions and washed once with PBS. Finally, cells were resuspended in PBS, and the fluorescence intensity was analysed using flow cytometry at excitation and emission wavelengths of 488 and 525 nm, respectively.
GSH/oxidised GSH ratio
Request a detailed protocolCells (5×105 cells/well) were seeded into 10 cm wells until they reached 80% confluence. One hour after 6 Gy irradiation, the GSH/oxidised GSH (GSSG) ratio was determined using a GSH/GSSG ratio detection assay kit following the manufacturer’s protocols.
Enzyme activity
Request a detailed protocolCells (5×105 cells/well) were seeded into 10 cm wells until reaching 80% confluence before. One hour after irradiation with 6 Gy, crude enzyme extract was prepared. The cells were digested, collected, and centrifuged to remove the supernatant. Next, 1 mL of HPLC-grade water was added to each pellet. The cells were disrupted by ultrasonication (20% power, 3 s on and 10 s off, repeated 25–40 times). The resulting mixture was centrifuged at 8000×g and 4°C for 10 min, and the supernatant was collected as the crude enzyme extract. The enzyme activities of SOD1, SOD2, CAT, and POD were detected using respective enzyme activity kits, according to the manufacturer’s protocols.
Statistical analysis
Request a detailed protocolData are expressed as the mean ± standard deviation (SD), and p-values<0.05 were considered statistically significant in all experiments. Data were analysed using one-way analysis of variance (ANOVA), Spearman’s correlation, and Kaplan-Meier estimates. Statistical analyses were performed using SPSS software (version 20.0). All experiments were performed in triplicate.
Data availability
Raw mRNA expression profiles and clinical features of the GSE9348, GSE20916, GSE37364, GSE44076, GSE68468 and GSE110223 datasets are available in the GEO database (http://www.ncbi.nlm.nih.gov/geo/). The raw mRNA expression profiles of rectal and colon cancer patients in the TCGA database are available ualcan database (https://ualcan.path.uab.edu/). The mRNA sequencing data supporting the findings of this study are available in the China National Center for Bioinformation (CNCB) under accession number PRJCA030073.
-
NCBI BioProjectID PRJCA030073. CPT1A Mediates Radiation Sensitivity in Colorectal Cancer.
-
NCBI Gene Expression OmnibusID GSE9348. Expression data from healthy controls and early stage CRC patient's tumor.
-
NCBI Gene Expression OmnibusID GSE20916. Modeling oncogenic signaling in colon tumors by multidirectional analyses of microarray data.
-
NCBI Gene Expression OmnibusID GSE37364. Expression data from human colonic biopsy samples (adenoma-carcinoma).
-
NCBI Gene Expression OmnibusID GSE44076. Gene expression data from healthy, adjacent normal and tumor colon cells.
-
NCBI Gene Expression OmnibusID GSE68468. caArray_notte-00422: Molecular Dissection of Colon Cancer.
-
NCBI Gene Expression OmnibusID GSE110223. Expression data from 13 patients with colorectal cancer.
References
-
Heat shock proteins and cancer: The FoxM1 connectionBiochemical Pharmacology 211:115505.https://doi.org/10.1016/j.bcp.2023.115505
-
ROS in cancer: the burning questionTrends in Molecular Medicine 23:411–429.https://doi.org/10.1016/j.molmed.2017.03.004
-
Targeting antioxidants for cancer therapyBiochemical Pharmacology 92:90–101.https://doi.org/10.1016/j.bcp.2014.07.017
-
Rectal cancer: ESMO Clinical Practice Guidelines for diagnosis, treatment and follow-upAnnals of Oncology 28:iv22–iv40.https://doi.org/10.1093/annonc/mdx224
-
Regulation of antioxidants in cancerMolecular Cell 84:23–33.https://doi.org/10.1016/j.molcel.2023.11.001
-
Inhibition of Cpt1a alleviates oxidative stress-induced chondrocyte senescence via regulating mitochondrial dysfunction and activating mitophagyMechanisms of Ageing and Development 205:111688.https://doi.org/10.1016/j.mad.2022.111688
-
FOXM1 and cancer: faulty cellular signaling derails homeostasisFrontiers in Oncology 10:626836.https://doi.org/10.3389/fonc.2020.626836
-
FOXM1: A small fox that makes more tracks for cancer progression and metastasisSeminars in Cancer Biology 92:1–15.https://doi.org/10.1016/j.semcancer.2023.03.007
-
Acyl-CoA synthetase-4 mediates radioresistance of breast cancer cells by regulating FOXM1Biochemical Pharmacology 192:114718.https://doi.org/10.1016/j.bcp.2021.114718
-
The critical role of dysregulated RhoB signaling pathway in radioresistance of colorectal cancerInternational Journal of Radiation Oncology, Biology, Physics 104:1153–1164.https://doi.org/10.1016/j.ijrobp.2019.04.021
-
Three differential expression analysis methods for RNA sequencing: limma, EdgeR, DESeq2Journal of Visualized Experiments 1:62528.https://doi.org/10.3791/62528
-
Cpt1a promoted ROS-induced oxidative stress and inflammation in liver injury via the Nrf2/HO-1 and NLRP3 inflammasome signaling pathwayCanadian Journal of Physiology and Pharmacology 99:468–477.https://doi.org/10.1139/cjpp-2020-0165
-
Carnitine palmitoyltransferase I in human carcinomas: a novel role in histone deacetylation?Cancer Biology & Therapy 6:1606–1613.https://doi.org/10.4161/cbt.6.10.4742
-
The carnitine system and cancer metabolic plasticityCell Death & Disease 9:228.https://doi.org/10.1038/s41419-018-0313-7
-
JASPAR 2024: 20th anniversary of the open-access database of transcription factor binding profilesNucleic Acids Research 52:D174–D182.https://doi.org/10.1093/nar/gkad1059
-
CPT1A-mediated Fat OxidationMechanisms, and Therapeutic Potential. Endocrinology 161:bqz046.https://doi.org/10.1210/endocr/bqz046
-
Radiation resistance: Cancer stem cells (CSCs) and their enigmatic pro-survival signalingSeminars in Cancer Biology 35:39–44.https://doi.org/10.1016/j.semcancer.2015.09.009
-
FOXM1-mediated regulation of reactive oxygen species and radioresistance in oral squamous cell carcinoma cellsLaboratory Investigation; a Journal of Technical Methods and Pathology 103:100060.https://doi.org/10.1016/j.labinv.2022.100060
-
GEPIA2: an enhanced web server for large-scale expression profiling and interactive analysisNucleic Acids Research 47:W556–W560.https://doi.org/10.1093/nar/gkz430
-
The recovery of intestinal barrier function and changes in oral microbiota after radiation therapy injuryFrontiers in Cellular and Infection Microbiology 13:1288666.https://doi.org/10.3389/fcimb.2023.1288666
-
FOXM1 regulates radiosensitivity of lung cancer cell partly by upregulating KIF20AEuropean Journal of Pharmacology 833:79–85.https://doi.org/10.1016/j.ejphar.2018.04.021
-
hTFtarget: a comprehensive database for regulations of human transcription factors and their targetsGenomics, Proteomics & Bioinformatics 18:120–128.https://doi.org/10.1016/j.gpb.2019.09.006
Article and author information
Author details
Funding
National Natural Science Foundation of China (82273564)
- Yi Ding
National Natural Science Foundation of China (32370139)
- Hongying Fan
National Natural Science Foundation of China (32300085)
- Zhenhui Chen
National Natural Science Foundation of China (32070118)
- Hongying Fan
Key Science and Technology Brainstorm Project of Guangzhou (202206010045)
- Yi Ding
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
We would like to thank Editage (https://www.editage.cn) for English language editing. This research project was supported by the National Natural Science Foundation of China (No. 82273564, 32370139, 32300085, and 32070118), Key Science & Technology Brainstorm Project of Guangzhou (No. 202206010045).
Ethics
Human samples were collected from the Nanfang Hospital of Southern Medical University, Guangzhou, for colorectal cancer patients with or without radiochemotherapy with age over 18 and below 75 year. Two to four biopsies from each patient were collected to be used for qRT-PCR, western blotting, and sectioning. Informed consent forms were acquired from all the patients. Ethics approval for the use of human samples was obtained from Ethics Committee of Nanfang Hospital, Southern Medical University (No. NFEC-202304-K13).
All animal experiments were approved by the Institutional Animal Care and Use Committee of Nanfang Hospital (IACUC-PROJECT-20221128-003).
Version history
- Sent for peer review:
- Preprint posted:
- Reviewed Preprint version 1:
- Reviewed Preprint version 2:
- Version of Record published:
Cite all versions
You can cite all versions using the DOI https://doi.org/10.7554/eLife.97827. This DOI represents all versions, and will always resolve to the latest one.
Copyright
© 2024, Chen, Yu, Zheng et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Cancer Biology
- Evolutionary Biology
In growing cell populations such as tumours, mutations can serve as markers that allow tracking the past evolution from current samples. The genomic analyses of bulk samples and samples from multiple regions have shed light on the evolutionary forces acting on tumours. However, little is known empirically on the spatio-temporal dynamics of tumour evolution. Here, we leverage published data from resected hepatocellular carcinomas, each with several hundred samples taken in two and three dimensions. Using spatial metrics of evolution, we find that tumour cells grow predominantly uniformly within the tumour volume instead of at the surface. We determine how mutations and cells are dispersed throughout the tumour and how cell death contributes to the overall tumour growth. Our methods shed light on the early evolution of tumours in vivo and can be applied to high-resolution data in the emerging field of spatial biology.
-
- Cancer Biology
- Evolutionary Biology
In asexual populations that don’t undergo recombination, such as cancer, deleterious mutations are expected to accrue readily due to genome-wide linkage between mutations. Despite this mutational load of often thousands of deleterious mutations, many tumors thrive. How tumors survive the damaging consequences of this mutational load is not well understood. Here, we investigate the functional consequences of mutational load in 10,295 human tumors by quantifying their phenotypic response through changes in gene expression. Using a generalized linear mixed model (GLMM), we find that high mutational load tumors up-regulate proteostasis machinery related to the mitigation and prevention of protein misfolding. We replicate these expression responses in cancer cell lines and show that the viability in high mutational load cancer cells is strongly dependent on complexes that degrade and refold proteins. This indicates that the upregulation of proteostasis machinery is causally important for high mutational burden tumors and uncovers new therapeutic vulnerabilities.