Novel biomarkers to predict treatment response and prognosis in locally advanced rectal cancer undergoing neoadjuvant chemoradiotherapy
BMC Cancer volume 23, Article number: 1099 (2023)
To identify genes associated with treatment response and prognosis for locally advanced rectal cancer (LARC) patients receiving neoadjuvant chemoradiotherapy (NCRT).
In our cohort, gene expression profiles of 64 tumor biopsy samples before NCRT were examined and generated. Weighted gene co-expression network analysis was performed to identify gene modules. External validation datasets included GSE3493, GSE119409, and GSE133057. The expression of candidate genes was evaluated using immunohistochemistry (IHC). TIMER was used to assess immune infiltration.
We identified and validated the capability to predict the treatment response of CCT5 and ELF1 using our data and external validation datasets. The trends of survival differences of candidate genes in the GSE133057 dataset were similar to our cohort. High levels of CCT5 and ELF1 expression were associated with NCRT resistance and poor prognosis. Furthermore, the expression of CCT5 and ELF1 were also assessed in 117 LARC patients’ samples by the IHC method. Based on IHC results and Cox analysis, the risk score model with CCT5 and ELF1 was constructed and performed well. The risk score was an independent prognostic factor for progression-free survival and overall survival in LARC patients and was then used to build nomogram models. The underlying mechanisms of CCT5 and ELF1 were explored using gene set enrichment analysis. The underlying pathway including apoptosis, cell cycle, and other processes. CCT5 and ELF1 expressions were significantly correlated with immune cell infiltration.
CCT5 and ELF1 were determined as biomarkers for treatment response and prognosis in LARC patients. The risk score model and nomograms helped predict treatment response and survival outcomes for LARC patients undergoing NCRT.
Colorectal cancer remains a leading disease burden and ranks second in cancer-related mortality worldwide [1, 2]. Preoperative neoadjuvant chemoradiotherapy (NCRT) presents advantageous features in tumor downstaging, surgical resection rates, and anus conservation rates [3,4,5]. Currently, NCRT followed by radical surgery is a mainstay of therapeutic strategy for patients with locally advanced rectal cancer (LARC). However, roughly 15-45% of patients exhibit therapy resistance and are likely to suffer from potential complications and toxicity of NCRT, which must not be neglected . Thus, identifying the regulator genes of NCRT is crucial to improving the treatment efficacy.
Previous studies have reported differentially expressed genes (DEGs) between NCRT-sensitive patients and NCRT-resistant patients and potential markers to predict tumor regression grading (TRG) [7,8,9,10]. Whereas, the prognostic value of the predictor was not explored in several studies. The neoadjuvant rectal (NAR) score has recently been proposed as a composite endpoint for LARC patients to predict clinical outcomes [11, 12]. A low NAR score indicates a positive treatment response and a better prognosis. In view of that, joint analysis of NAR score, TRG, and prognosis may contribute to a better understanding of NCRT regulatory factors and clinical outcomes.
In this study, we performed gene expression profiles on tumor biopsy samples from patients with LARC undergoing NCRT. Weighted gene co-expression network analysis (WGCNA) was used to determine NAR score-related modules and identify candidate genes with predictive and prognostic significance. Patient tissue samples and external datasets were used for validation.
Patients and clinical data collection
A total of 64 patients with LARC undergoing NRCT between 2015 and 2018 in Fujian Medical University Union Hospital were enrolled. Specifically, as we previously reported , the radiotherapy consisted of a 45 Gy dose in 25 fractions over five weeks and a boost dose of 5.4 Gy for the tumor. And concurrent chemotherapy was as follows: oral capecitabine 825 mg/m2 twice per day for two weeks. Patients received radical surgery after 6 to 8 weeks from the last dose of radiotherapy. All patients were recommended to receive adjuvant chemotherapy after surgery. Tumor biopsy samples, used for gene expression profiles analysis, were obtained from colonoscopy before NCRT. In addition, LARC patients undergoing NRCT and radical surgery between 2012 and 2014 were also included. Their colonoscopy samples before neoadjuvant treatment were collected to validate and further identify the protein expression of candidate genes. This study was approved by the Institutional Review Board of Fujian Medical University Union Hospital (2019KY006).
TRG was used to assess pathological response to NCRT. In detail, TRG 0, no residual tumor cells; TRG 1, near-complete regression with tumor cells individually or in small groups; TRG2, residual tumor cells with a desmoplastic response; and TRG 3, minimal or no regression. Patients with TRG 0 and TRG 1 were classified as NCRT-sensitive groups, while TRG 2 and TRG 3 as NCRT-resistant groups. NAR score was calculated based on the equation: [5ypN–3 (cT–ypT) + 12]² /9.61 . Hereinto, cT refers to clinical T stage (value: 1, 2, 3, 4), ypT refers to pathological T stage (value: 0, 1, 2, 3, 4), and ypN refers to pathological nodal status (value: 0, 1, 2).
Gene expression analysis
The total RNA was extracted from tumor biopsy samples using Trizol reagent (Invitrogen) according to the manufacturer’s instructions. NanoDrop ND-1000 monitored RNA quality control and quantification. Labeling, hybridization, and scanning were carried out according to standard protocols. Then, data quality control and normalization were performed.
DEGs identification and enrichment analysis
The DEGs were identified using the R “limma” package. Genes with |log fold change (logFC)|> 0.5 and p-value < 0.05 were determined as DEGs between the NCRT-sensitive and the NCRT-resistant group. Visualization and comparison of DEGs use volcano plots and heatmaps. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) were used to explore the biological functions of DEGs by the R “clusterProfiler” package.
WGCNA was carried out using the R “WGCNA” package . That is, the soft threshold power was set, and the topological overlap matrix (TOM) was calculated. Modules were determined by the Dynamic Tree Cut method. Then, high similarity modules were merged by clustering analysis. The correlation between modules and clinicopathological features is calculated.
The protein expression of the candidate genes was assessed by immunohistochemistry. Immunohistochemical staining for candidate genes was conducted as described earlier . The following antibodies were used: anti-FBXO7 (203,049-T40, Sino Biological, China), anti-GSTT4 (bs-16345R, Bioss, China), anti-CCT5 (11603-1-AP, Proteintech, China), and anti-ELF1 (22565-1-AP, Proteintech, China), and anti-SLC44A1 antibody (14687-1-AP, Proteintech, China). Immunohistochemical results were scored using a semi-quantitative scoring method. Specifically, data were collected from random visual fields of five different areas. The intensity of staining was scored as 0 (negative), 1 (light yellow), 2 (brown), and 3 (deep brown). The percentage of positive cells was scored as 0 (< 5%), 1 (5–25%), 2 (25–50%), 3 (50–75%), and 4 (> 75%). The two values were multiplied and calculated as immunohistochemical scores. Scores of 0–4 were considered low expression, and those with scores above 4 were classified as high expression.
External validation datasets
External validation datasets were obtained from Gene Expression Omnibus (GEO) database. The GSE3493, including 46 LARC patients undergoing NRCT, was identified to assess the gene expression of the candidate genes between NCRT-sensitive groups and NCRT-resistant groups. After excluding samples with unknown TRG information, a total of 56 LARC patients undergoing NRCT in GSE119409 were also included in the external validation. The GSE133057 (33 LARC patients) was adopted for the validation of survival outcomes.
Gene set enrichment analysis
Gene set enrichment analysis (GSEA) was employed to investigate the potential biological pathways of candidate genes. LARC patients were enrolled in high and low-expression groups of candidate genes based on the median expression. The false discovery rate (FDR) < 0.25 and P < 0.05 were accepted as statistically significant.
Immune infiltration analysis
The relationship between gene expression and immune cells infiltration (including B cell, CD8 + T cell, CD4 + T cell, macrophage, neutrophil, and dendritic cell) was explored using Tumor IMmune Estimation Resource (TIMER).
The drug sensitivity dataset was obtained from the CellMiner website . Pearson’s correlation test determined the correlation between target genes and drug sensitivity. These correlation results were visualized using the R package ggplot2.
Statistical analysis was carried out with R software (version 4.1.2), GraphPad Prism 8, and SPSS (version 22). X-tile software was adopted to determine the cut-off points for the gene expression of the candidate genes. T-tests or non-parametric tests were done to test differences between continuous variables. The categorical data were analyzed by Fisher’s exact or Chi-square tests. Kaplan–Meier (KM) analysis was used to estimate survival outcomes using a log-rank test. Correlations were evaluated using Pearson correlation test. Receiver operating characteristics (ROC) analysis was performed, and the area under the curve (AUC) was calculated to assess the TRG and predictive survival capacities of the candidate genes. Cox proportional hazards regression analysis identified independent risk factors for progression-free survival (PFS) and overall survival (OS). A nomogram was created based on the above factors by the R “rms” package. The calibration curve was used to evaluate the performance of the model. P-value < 0.05 was considered statistically significant.
DEGs and functional enrichment
The gene expression profile based on gene chips was obtained from 64 LARC patients (Supplementary Table 1) before NCRT in our cohort. A total of 333 DEGs were identified between NCRT-sensitive patients and NCRT-resistant patients. Among these, compared to the NCRT-sensitive group, 94 genes are upregulated, and 239 genes are downregulated in the NCRT-resistant group. Heatmap and volcano plot were shown in Fig. 1A and B. We performed GO and KEGG analyses to explore the potential biological significance of DEGs. GO analyses revealed that the DEGs were significantly enriched in serine phosphorylation of STAT protein, natural killer cell activation involved in immune response, response to exogenous dsRNA, and so forth (Fig. 1C). Furthermore, the top 10 enriched KEGG pathways are presented in Fig. 1D, including cytokine-cytokine receptor interaction, the Toll-like receptor signaling pathway, and cytosolic DNA-sensing pathway, respectively.
Identification of candidate genes
A weighted gene co-expression network was constructed to identify NCRT-regulated genes (Fig. 2A). A total of 3 modules were identified by merging the high similarity modules. Here, as seen in Fig. 2B, a negative correlation of the brown module with TRG was present (r=-0.24, P = 0.06), and the blue module positively correlated with carcinoembryonic antigen (CEA, r = 0.25, P = 0.05). The turquoise module (r = 0.28, P = 0.02) had the highest positive correlation with the NAR score, while the brown module (r=-0.47, P < 0.01) with the highest negative correlation. Next, we also identified 1147 TRG-associated genes and 692 cancer progression genes(PFS-associated genes), respectively (all P < 0.05). The intersection among the three sets was estimated to determine the candidate genes with the predictive capabilities of both pathological response and prognosis. Finally, there were five genes overlapped among three sets, including FBXO7, CCT5, ELF1, GSTT4, and SLC44A1, which were identified as the candidate genes (Fig. 2C). In addition, GO and KEGG pathway analyses were performed for the turquoise module (Fig. 2D and E) and the brown module (Fig. 2F and G) genes to gain a more comprehensive understanding of biological effects.
Internal validation of predictive capabilities of candidate genes
Next, to test the discriminant power of candidate genes in NCRT sensitivity, the correlations between candidate genes’ expression and TRG grade were calculated. The R-values of the correlation for FBXO7, CCT5, ELF1, GSTT4, and SLC44A1 were 0.27 (P = 0.032), 0.27 (P = 0.028), 0.26 (P = 0.041), -0.25 (P = 0.047), and 0.21 (P = 0.092), respectively (Fig. 3A and E). We also evaluated the association between candidate genes and NAR score, and the results revealed a significant correlation between these (Fig. 3F J). Next, we employed the gene expression of five candidate genes to predict the pathologic response of NCRT. As shown in Fig. 3K and O, the AUC values of FBXO7, CCT5, ELF1, GSTT4, and SLC44A1 to predict NCRT response were 0.637 (P = 0.059), 0.654 (P = 0.034), 0.638 (P = 0.057), 0.632 (P = 0.071), and 0.672 (P = 0.018), respectively. Furthermore, the ability to assess the prognosis of five candidate genes was further explored. The cut-off points of gene expression were determined by X-tile software for survival analysis (Supplementary Fig. 1). KM analysis revealed the candidate genes expression at diagnosis could predict the survival of LARC patients undergoing NRCT. Specifically, high expression of FBXO7, CCT5, ELF1, and SLC44A1 were correlated with poor survival, while high expression of GSTT4 was associated with a better prognosis (Fig. 3P and Y). Notably, no significant overall survival differences were observed in high ELF1 expression (P = 0.242).
External validation analysis
The GSE3493 and GSE119409 datasets were collected from the GEO database for external validation. Gene expression was applied to predict the treatment response (NCRT-sensitive). As depicted in Supplementary Fig. 2A-2D, in GSE3943 dataset, the AUC of FBXO7, CCT5, ELF1 and GSTT4 were 0.636 (P = 0.18), 0.512 (P = 0.91), 0.697(P = 0.05), and 0.575 (P = 0.46), respectively. As displayed in Supplementary Fig. 2E-2I, in GSE119409 dataset, CCT5 showed the highest power (AUC = 0.689, P = 0.03), followed by FBXO7 (AUC = 0.593, P = 0.29), SLC44A1 (AUC = 0.586, P = 0.33) GSTT4 (AUC = 0.527, P = 0.76) and ELF1 (AUC = 0.515, P = 0.86). Furthermore, we also evaluated the effects of candidate genes on clinical outcomes. Based on X-tile, candidate genes were separated as the high and low expression groups in GSE133057. We observed that patients of the CCT5 high expression groups had significantly worse prognoses (P = 0.012) and that a trend toward poor survival in the high expression of FBXO7, ELF1, GSTT4, SLC44A1 (all P > 0.05,Supplementary Fig. 2J-2 N). However, due to SLC44A1 of GSE3493 being missing, we were unable to verify its performance.
Immunohistochemistry validation of candidate genes
To further validate the gene expression of five candidate genes in LARC patients undergoing NCRT, a total of 117 patients were enrolled in the validation cohort. Clinicopathological features of the validation cohort were shown in Supplementary Table 2. We assessed the protein expression of five candidate genes by immunohistochemical staining in tumor biopsy samples from these patients before NCRT. Supplementary Fig. 3 illustrated the different expression levels of five genes. The result demonstrated a higher immunohistochemistry score of CCT5 and ELF1 in NCRT-resistant patients compared with NCRT-sensitive patients (all P < 0.05,Fig. 4B C), while no difference was observed in the expression of FBXO7, SLC44A1, and GSTT4 (Fig. 4A, D and E). The ROC curve indicated the AUC of CCT5 and ELF1 to predict the treatment response to NCRT were 0.727 and 0.717 (Fig. 4G H, all P < 0.05), and FBXO7, GSTT4, and SLC44A1 cannot differentiate between two groups (Fig. 4F J, all P > 0.05). Moreover, CCT5 and ELF1 expression were significantly related to the NAR score (Fig. 4L M), while the trend was not observed in FBXO7, GSTT4, and SLC44A1 (Fig. 4K N-4O). Survival analysis demonstrated that high expression of FBXO7, CCT5, and ELF1 had worse PFS (P = 0.010, P = 0.058, P = 0.018, respectively, Fig. 4P and R) and OS (P = 0.043, P = 0.069, P = 0.003, respectively,Fig. 4U W). At the same time, there was no difference in survival between high and low GSTT4 and SLC44A1 expression, (Fig. 4S, T, X and Y).
Construction of a risk score model
We further performed Cox analysis to determine the prognostic value of genes. On univariate analysis, a high CCT5 and ELF1 expression were associated with shorter OS (all P < 0.05). Multivariate analysis revealed that CCT5 (hazard ratio [HR], 1.141, 95% confidence interval [CI]: 1.007–1.294, P = 0.039) and ELF1 (HR,1.179, 95%CI, 1.067–1.304, P = 0.001) were independently associated with prognosis (Table 1). Hence, based on coefficients obtained from Cox analysis, the risk score was computed as: 0.132×CCT5 expression + 0.165×ELF1 expression. We divided the validation cohort into two groups based on the median risk score, and the clinicopathologic characteristics of the two groups were presented in Supplementary Table 3. The risk score of the pathological complete response (pCR) group was significantly lower than the non-pCR group (Fig. 5A), with excellent predictive capacity (AUC = 0.780, P < 0.01, Fig. 5B). A significant correlation was observed between the risk score and NAR score (R = 0.43, P < 0.01, Fig. 5C). Survival analysis showed the risk score has an excellent ability to discriminate clinical outcomes (Fig. 5D H).
Prognosis analyses and development nomogram
We after that evaluated the prognostic value of the risk score. Univariate analysis revealed that the ypTMN stage (P < 0.001), TRG (P < 0.001), pathological type (P = 0.036), risk score (P < 0.001), and NAR score (P = 0.006) were significantly associated with PFS. Multivariate analysis manifested that the ypTMN stage ( P = 0.040) and risk score ( P = 0.025) were independent risk factors of PFS (Table 2). Furthermore, as shown in Table 3, the ypTMN stage ( P = 0.031) and risk score ( P = 0.008) were significantly correlated with OS by multivariate analysis. Results of the analysis based on Cox analysis, nomograms were constructed to predict 1-year, 3-year-, and 5-year PFS (Fig. 6A) and OS (Fig. 6C), and good calibration was also confirmed (Fig. 6B and D).
External validation for risk score
The risk score of the NCRT-resistant group was considerably higher in the GSE3493 cohort than the NCRT-sensitive group ( Supplementary Fig. 4A) with an adequate capacity to discriminate NCRT-sensitive patients (AUC = 0.701, P = 0.046, Supplementary Fig. 4C). However, we found no statistical difference between the gene expression of the GSE119409 cohort ( Supplementary Fig. 4B, 4D) and the prognosis of the GSE133057 cohort (all P > 0.05, Supplementary Fig. 4E).
Following that, we look into the potential mechanisms by which CCT5 and ELF1 might impact NCRT sensitivity. The median expression level defined the groups with high and low CCT5 and ELF1 expression. Colorectal cancer, apoptosis, and DNA replication were enriched in the high CCT5 expression group, according to GSEA analysis (Fig. 7A). Furthermore, the high ELF1 expression group was significantly enriched with apoptosis, cell cycle, mTOR signaling pathway, and cancer pathway (Fig. 7B). These pathways contribute to a better understanding of the mechanism of NCRT resistance.
Drug sensitivity analysis
We then examined the association of CCT5 and ELF1 with drug sensitivity using the CellMiner database. Figure 8 demonstrated the results of the drug sensitivity analysis. For instance, the ELF1 expression was positively associated with the drug sensitivity of vorinostat, artesunate, nilotinib, and selumetinib (All P < 0.05).
Analysis of immune infiltration
Next, we also investigated immune cell infiltration. As shown in Fig. 9A, positive correlations were observed between CCT5 expression and the infiltration of CD8 + T cell (P < 0.01) and neutrophil (P < 0.05). B cell, CD8 + T cell, macrophage, and neutrophil infiltration were significantly associated with ELF1 expression (All P < 0.05, Fig. 9B).
Previous studies have demonstrated the superiority of NCRT in LARC treatment. At the moment, the mechanism of NCRT resistance remains unclear. Given treatment resistance in certain patients and potential complications and toxicity, further research is warranted. The expression profile of the pre-treatment tumor sample from our cohort was used for the analyses in this study. WGCNA was utilized to determine NAR score-related modules and genes. To determine the potential genes with a high ability to predict sensitivity to NCRT and prognosis, TRG-related and cancer progression (also PFS)-related genes were also identified. The intersection of three gene sets was regarded as candidate genes for additional validation at the protein level. The protein expression of CCT5 and ELF1 were significantly different in the NCRT-sensitive and NCRT-resistant groups according to the results of immunohistochemistry staining of candidate genes. Furthermore, high-level expression of CCT5 and ELF1 were linked to a poor prognosis.
The NAR score has recently been considered a potential surrogate endpoint for the prognosis of LARC patients undergoing NCRT . This score formula was established on the basis of Valentini et al. nomogram for predicting the overall survival in LARC patients, and it considered both pretreatment tumor burden and posttreatment tumor regression . Previous researches have demonstrated that the NAR score can accurately predict clinical outcomes and assist in assessing the necessity of adjunctive therapy [18,19,20]. Thus, we found NAR score-related genes using WCGNA, a more biological approach to correlating modules with phenotypic features. However, there is some debate as to the prognostic value of the NAR score. A study from the Netherlands demonstrated that a combination model based on clinical data and pathological data exceeded the NAR score in assessing survival outcomes . As a result, cancer progression (also PFS)-associated genes were also discovered. Therefore, discerning therapy response before the surgical procedure aided clinical decision-making. After adequate evaluation, Watch & Wait and local excision might become a suitable therapy for NCRT-sensitive patients, maximizing the patient and enhancingn appropriat the quality of life . Hence, we identify TRG-associated genes, and the intersection of three gene sets contributes to determining the genes with the value of predicting prognosis and treatment response.
FBXO7, CCT5, ELF1, GSTT4, and SLC44A1 were discovered as biomarkers for treatment response and prognosis in our cohort. The value of CCT5 and ELF1 in predicting response to therapy was validated in the GSE3493 and GSE119409 cohorts. The survival analysis of the GSE133057 cohort showed that increased expression of FBXO7, CCT5, ELF1, and SLC44A1 implies poor survival, whereas high expression of GSTT4 predicts the opposite outcomes. These survival trends match those obtained in our cohort. CCT5 and ELF1 were eventually determined after the IHC validation of protein levels.
CCT5 is a protein folding subunit of the chaperonin containing TCP1 complex (also known as TCP1-ring complex). Recently, accumulating evidence reveals that CCT5 was implicated in tumor progression, and high expression of CCT5 has been discovered in a range of malignancies and associated with worse survival [23,24,25]. Notably, CCT5 is involved in chemotherapy resistance. In breast cancer, the knockdown of CCT5 leads to increased apoptosis after docetaxel treatment . CCT5 expression was elevated in multidrug-resistant gastric carcinoma cells . Nevertheless, the functions of CCT5 in NCRT of LARC patients remain elusive. Our studies demonstrate that high CCT5 expression correlated with NAR score, NCRT resistance, and poor prognosis in LARC patients.
ELF1 is a transcription factor that belongs to the ETS family and plays contrasting roles in different tumors. Some reports pointed out the part of ELF1 has tumor-promoting effects in glioma [28, 29], oral squamous cell carcinoma , choroidal melanoma , endometrial carcinoma , acute myeloid leukemia , nasopharyngeal carcinoma , and osteosarcoma , while others have reported a tumor-suppressive function in Hodgkin lymphoma  and prostate cancer . Starr et al. found that ELF1 regulates the expression of TM9SF2, an oncogene in colorectal cancer . Herein, our findings indicated that ELF1 expression was more strongly raised in NCRT-resistant patients, and ELF1 was an independent predictor of survival. Going further, GSEA analysis revealed apoptosis, cell cycle, and mTOR signaling pathway are significantly enriched in high ELF1 expression, and those pathways were implicated in treatment resistance [39,40,41].
Compared to a single prognostic factor, the advantages of the risk model and nomogram have previously been reported [42, 43]. In our study, we build a risk score model based on CCT5 and ELF1 IHC expression that effectively discriminated between non-pCR and pCR patients and survival. Nomograms were created by taking together the result of Cox regression analysis. These models are helpful in predicting the prognosis of LARC patients undergoing NCRT.
Furthermore, a growing body of studies has revealed the impact of tumor microenvironment on chemoradiotherapy sensitivity. One study reported that increased CD163 + tumor-associated macrophages (TAM) were observed in non-pCR group in LARC patients . Meanwhile, TAM infiltration was confirmed with chemoradiotherapy resistance in oral squamous cell carcinoma and cervical cancer [45, 46]. In addition, several previous studies have revealed that increased neutrophils were associated with poor chemoradiotherapy response and worse clinical outcomes [47,48,49,50]. Here, our study also explored the potential role of CCT5 and ELF1 in the tumor microenvironment and found the infiltration level of CD8 + T cell and neutrophil were positively associated with CCT5 and ELF1 expression, and the same trend was observed in the correlation between B cell and macrophage and ELF1 expression. However, there are limited studies about CCT5 or ELF1 leading to immune infiltration, and the specific mechanism remains to be elucidated. Nonetheless, these findings provided valuable information and potential directions for future research.
There are limitations to this study. More external datasets are required to confirm and validate our findings. Meanwhile, bioinformatics has deconstructed the probable processes of CCT5 and ELF1, but further investigation will be needed to clarify their biological role and impact on the cancer microenvironment in future studies.
In summary, CCT5 and ELF1 were determined as biomarkers for NCRT treatment response and prognosis by internal and external validation. The risk score model and nomogram were constructed to predict survival. These findings contributed to personalized clinical decision-making for LARC patients undergoing NCRT.
External validation datasets (GSE3493, GSE119409, and GSE133057) were obtained from Gene Expression Omnibus (http://www.ncbi.nlm.nih.gov/geo/). The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.
Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, Bray F. Global Cancer Statistics 2020: GLOBOCAN estimates of incidence and Mortality Worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2021;71(3):209–49.
Siegel RL, Miller KD, Goding Sauer A, Fedewa SA, Butterly LF, Anderson JC, Cercek A, Smith RA, Jemal A. Colorectal cancer statistics, 2020. CA Cancer J Clin. 2020;70(3):145–64.
Li Y, Wang J, Ma X, Tan L, Yan Y, Xue C, Hui B, Liu R, Ma H, Ren J. A review of Neoadjuvant Chemoradiotherapy for locally advanced rectal Cancer. Int J Biol Sci. 2016;12(8):1022–31.
Ominelli J, Valadão M, Araujo ROC, Cristina de Melo A, Araujo LH. The Evolving Field of Neoadjuvant Therapy in locally-advanced rectal Cancer: evidence and prospects. Clin Colorectal Cancer. 2021;20(4):288–98.
Gollins S, Sebag-Montefiore D. Neoadjuvant treatment strategies for locally advanced rectal Cancer. Clin Oncol (R Coll Radiol (G B)). 2016;28(2):146–51.
Ha YJ, Tak KH, Kim CW, Roh SA, Choi EK, Cho DH, Kim JH, Kim SK, Kim SY, Kim YS, et al. PSMB8 as a candidate marker of responsiveness to Preoperative Radiation Therapy in rectal Cancer patients. Int J Radiat Oncol Biol Phys. 2017;98(5):1164–73.
Seo I, Lee HW, Byun SJ, Park JY, Min H, Lee SH, Lee JS, Kim S, Bae SU. Neoadjuvant chemoradiation alters biomarkers of anticancer immunotherapy responses in locally advanced rectal cancer. J Immunother Cancer 2021, 9(3).
Momma T, Okayama H, Kanke Y, Fukai S, Onozawa H, Fujita S, Sakamoto W, Saito M, Ohki S, Kono K. Validation of Gene Expression-Based Predictive Biomarkers for Response to Neoadjuvant Chemoradiotherapy in Locally Advanced Rectal Cancer. Cancers (Basel) 2021, 13(18).
Agostini M, Zangrando A, Pastrello C, D’Angelo E, Romano G, Giovannoni R, Giordan M, Maretto I, Bedin C, Zanon C, et al. A functional biological network centered on XRCC3: a new possible marker of chemoradiotherapy resistance in rectal cancer patients. Cancer Biol Ther. 2015;16(8):1160–71.
Millino C, Maretto I, Pacchioni B, Digito M, De Paoli A, Canzonieri V, D’Angelo E, Agostini M, Rizzolio F, Giordano A, et al. Gene and MicroRNA expression are Predictive of Tumor response in rectal adenocarcinoma patients treated with preoperative chemoradiotherapy. J Cell Physiol. 2017;232(2):426–35.
Glynne-Jones R, Glynne-Jones S. The concept and use of the neoadjuvant rectal score as a composite endpoint in rectal cancer. Lancet Oncol. 2021;22(7):e314–26.
Fokas E, Fietkau R, Hartmann A, Hohenberger W, Grützmann R, Ghadimi M, Liersch T, Ströbel P, Grabenbauer GG, Graeven U, et al. Neoadjuvant rectal score as individual-level surrogate for disease-free survival in rectal cancer in the CAO/ARO/AIO-04 randomized phase III trial. Ann Oncol. 2018;29(7):1521–7.
Li A, Huang T, Zheng R, Chi P, Li Z, Wang X, Xu B. Preoperative chemoradiotherapy with capecitabine and triweekly oxaliplatin versus capecitabine monotherapy for locally advanced rectal cancer: a propensity-score matched study. BMC Cancer. 2022;22(1):789.
Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. 2008;9:559.
Feng C, Zhang Y, Huang J, Zheng Q, Yang Y, Xu B. The Prognostic significance of APOBEC3B and PD-L1/PD-1 in nasopharyngeal carcinoma. Appl Immunohistochem Mol Morphology: AIMM. 2021;29(3):239–44.
Shankavaram UT, Varma S, Kane D, Sunshine M, Chary KK, Reinhold WC, Pommier Y, Weinstein JN. CellMiner: a relational database and query tool for the NCI-60 cancer cell lines. BMC Genomics. 2009;10:277.
Valentini V, van Stiphout RG, Lammering G, Gambacorta MA, Barba MC, Bebenek M, Bonnetain F, Bosset JF, Bujko K, Cionini L, et al. Nomograms for predicting local recurrence, distant metastases, and overall survival for patients with locally advanced rectal cancer on the basis of european randomized clinical trials. J Clin Oncol. 2011;29(23):3163–72.
Lim YJ, Song C, Jeon SH, Kim K, Chie EK. Risk stratification using neoadjuvant rectal score in the era of Neoadjuvant Chemoradiotherapy: Validation with Long-term Outcome Data. Dis Colon Rectum. 2021;64(1):60–70.
Huang WS, Kuan FC, Lin MH, Chen MF, Chen WC. Prognostic significance of Neoadjuvant rectal scores in Preoperative Short-Course Radiotherapy and Long-Course Concurrent Chemoradiotherapy for patients with locally advanced rectal Cancer. Ann Surg Oncol. 2020;27(11):4309–18.
Roselló S, Frasson M, García-Granero E, Roda D, Jordá E, Navarro S, Campos S, Esclápez P, García-Botello S, Flor B, et al. Integrating downstaging in the Risk Assessment of patients with locally advanced rectal Cancer treated with Neoadjuvant Chemoradiotherapy: validation of Valentini’s Nomograms and the neoadjuvant rectal score. Clin Colorectal Cancer. 2018;17(2):104–112e102.
van der Valk MJM, Vuijk FA, Putter H, van de Velde CJH, Beets GL, Hilling DE. Disqualification of neoadjuvant rectal score based on data of 6596 patients from the Netherlands Cancer Registry. Clin Colorectal Cancer. 2019;18(2):e231–6.
Smith JJ, Strombom P, Chow OS, Roxburgh CS, Lynn P, Eaton A, Widmar M, Ganesh K, Yaeger R, Cercek A, et al. Assessment of a Watch-and-wait strategy for rectal Cancer in patients with a complete response after Neoadjuvant Therapy. JAMA Oncol. 2019;5(4):e185896.
Meng Y, Yang L, Wei X, Luo H, Hu Y, Tao X, He J, Zheng X, Xu Q, Luo K, et al. CCT5 interacts with cyclin D1 promoting lung adenocarcinoma cell migration and invasion. Biochem Biophys Res Commun. 2021;567:222–9.
He J, McLaughlin RP, van der Beek L, Canisius S, Wessels L, Smid M, Martens JWM, Foekens JA, Zhang Y, van de Water B. Integrative analysis of genomic amplification-dependent expression and loss-of-function screen identifies ASAP1 as a driver gene in triple-negative breast cancer progression. Oncogene. 2020;39(20):4118–31.
Li Y, Liu C, Zhang X, Huang X, Liang S, Xing F, Tian H. CCT5 induces epithelial-mesenchymal transition to promote gastric cancer lymph node metastasis by activating the Wnt/β-catenin signalling pathway. Br J Cancer. 2022;126(12):1684–94.
Ooe A, Kato K, Noguchi S. Possible involvement of CCT5, RGS3, and YKT6 genes up-regulated in p53-mutated tumors in resistance to docetaxel in human breast cancers. Breast Cancer Res Treat. 2007;101(3):305–15.
Ludwig A, Dietel M, Lage H. Identification of differentially expressed genes in classical and atypical multidrug-resistant gastric carcinoma cells. Anticancer Res. 2002;22(6a):3213–21.
Cheng M, Zeng Y, Zhang T, Xu M, Li Z, Wu Y. Transcription factor ELF1 activates MEIS1 transcription and then regulates the GFI1/FBW7 Axis to promote the development of Glioma. Mol Therapy Nucleic Acids. 2021;23:418–30.
Hu M, Li H, Xie H, Fan M, Wang J, Zhang N, Ma J, Che S. ELF1 transcription factor enhances the progression of Glioma via ATF5 promoter. ACS Chem Neurosci. 2021;12(7):1252–61.
Qiao C, Qiao T, Yang S, Liu L, Zheng M. SNHG17/miR-384/ELF1 axis promotes cell growth by transcriptional regulation of CTNNB1 to activate Wnt/β-catenin pathway in oral squamous cell carcinoma. Cancer Gene Ther. 2022;29(1):122–32.
Wang L, Tang D, Wu T, Sun F. ELF1-mediated LUCAT1 promotes choroidal melanoma by modulating RBX1 expression. Cancer Med. 2020;9(6):2160–70.
Takai N, Miyazaki T, Nishida M, Shang S, Nasu K, Miyakawa I. Clinical relevance of Elf-1 overexpression in endometrial carcinoma. Gynecol Oncol. 2003;89(3):408–13.
Pang Y, Zhao Y, Wang Y, Wang X, Wang R, Liu N, Li P, Ji M, Ye J, Sun T, et al. TNFAIP8 promotes AML chemoresistance by activating ERK signaling pathway through interaction with Rac1. J Exp Clin Cancer Res. 2020;39(1):158.
Chen CH, Su LJ, Tsai HT, Hwang CF. ELF-1 expression in nasopharyngeal carcinoma facilitates proliferation and metastasis of cancer cells via modulation of CCL2/CCR2 signaling. Cancer Manage Res. 2019;11:5243–54.
Wang L. ELF1-activated FOXD3-AS1 promotes the migration, invasion and EMT of osteosarcoma cells via sponging mir-296-5p to upregulate ZCCHC3. J bone Oncol. 2021;26:100335.
Paczkowska J, Soloch N, Bodnar M, Kiwerska K, Janiszewska J, Vogt J, Domanowska E, Martin-Subero JI, Ammerpohl O, Klapper W, et al. Expression of ELF1, a lymphoid ETS domain-containing transcription factor, is recurrently lost in classical Hodgkin lymphoma. Br J Haematol. 2019;185(1):79–88.
Budka JA, Ferris MW, Capone MJ, Hollenhorst PC. Common ELF1 deletion in prostate cancer bolsters oncogenic ETS function, inhibits senescence and promotes docetaxel resistance. Genes & cancer. 2018;9(5–6):198–214.
Clark CR, Maile M, Blaney P, Hellweg SR, Strauss A, Durose W, Priya S, Habicht J, Burns MB, Blekhman R, et al. Transposon mutagenesis screen in mice identifies TM9SF2 as a novel colorectal cancer oncogene. Sci Rep. 2018;8(1):15327.
Adamsen BL, Kravik L, De Angelis PM. Cellular response to chemoradiotherapy, radiotherapy and chemotherapy in two colorectal cancer cell lines. Radiat Res. 2009;171(5):562–71.
Suzuki T, Sadahiro S, Tanaka A, Okada K, Saito G, Kamijo A, Akiba T, Kawada S. Predictive markers of chemoradiotherapy for rectal cancer: comparison of biopsy specimens taken before and about 1 week after the start of chemoradiotherapy. Int J Clin Oncol. 2015;20(6):1130–9.
Wanigasooriya K, Tyler R, Barros-Silva JD, Sinha Y, Ismail T, Beggs AD. Radiosensitising Cancer using Phosphatidylinositol-3-Kinase (PI3K), protein kinase B (AKT) or mammalian target of Rapamycin (mTOR) inhibitors. Cancers (Basel) 2020, 12(5).
Zheng W, Lu Y, Feng X, Yang C, Qiu L, Deng H, Xue Q, Sun K. Improving the overall survival prognosis prediction accuracy: a 9-gene signature in CRC patients. Cancer Med. 2021;10(17):5998–6009.
Hong T, Cai D, Jin L, Zhang Y, Lu T, Hua D, Wu X. Development and validation of a nomogram to predict survival after curative resection of nonmetastatic colorectal cancer. Cancer Med. 2020;9(12):4126–36.
Yang Y, Tian W, Su L, Li P, Gong X, Shi L, Zhang Q, Zhao B, Zhao H. Tumor-infiltrating cytotoxic T cells and Tumor-Associated Macrophages Correlate with the Outcomes of Neoadjuvant Chemoradiotherapy for locally advanced rectal Cancer. Front Oncol. 2021;11:743540.
Matsuoka Y, Yoshida R, Nakayama H, Nagata M, Hirosue A, Tanaka T, Kawahara K, Nakagawa Y, Sakata J, Arita H, et al. The tumour stromal features are associated with resistance to 5-FU-based chemoradiotherapy and a poor prognosis in patients with oral squamous cell carcinoma. APMIS: Acta Pathologica Microbiologica et immunologica Scandinavica. 2015;123(3):205–14.
Lippens L, Van Bockstal M, De Jaeghere EA, Tummers P, Makar A, De Geyter S, Van de Vijver K, Hendrix A, Vandecasteele K, Denys H. Immunologic impact of chemoradiation in cervical cancer and how immune cell infiltration could lead toward personalized treatment. Int J Cancer. 2020;147(2):554–64.
Hu P, Liu Q, Deng G, Zhang J, Liang N, Xie J, Zhang J. Radiosensitivity nomogram based on circulating neutrophils in thoracic cancer. Future Oncol. 2019;15(7):727–37.
Wisdom AJ, Hong CS, Lin AJ, Xiang Y, Cooper DE, Zhang J, Xu ES, Kuo HC, Mowery YM, Carpenter DJ, et al. Neutrophils promote tumor resistance to radiation therapy. Proc Natl Acad Sci U S A. 2019;116(37):18584–9.
Schernberg A, Nivet A, Dhermain F, Ammari S, Escande A, Pallud J, Louvel G, Deutsch E. Neutrophilia as a biomarker for overall survival in newly diagnosed high-grade glioma patients undergoing chemoradiation. Clin Translational Radiation Oncol. 2018;10:47–52.
Schernberg A, Vernerey D, Goldstein D, Van Laethem JL, Glimelius B, van Houtte P, Bonnetain F, Louvet C, Hammel P, Huguet F. Predictive value of neutrophils count for local Tumor Control after Chemoradiotherapy in patients with locally advanced pancreatic carcinoma. Int J Radiat Oncol Biol Phys. 2021;110(4):1022–31.
The authors appreciate the GEO database for open access to the data. We thank all the Department of Radiation Oncology staff, Fujian Medical University Union Hospital, Fuzhou, Fujian Province, People’s Republic of China.
This study was sponsored by Fujian provincial health technology project (2021CXA011), National Natural Science Foundation of Fujian Province, China (No.2022J02037), National Natural Science Foundation of China (No. 82172800), Special Financial Foundation of Fujian Provincial (No. 2020B019), Joint Funds for the Innovation of Science and Technology, Fujian Province (No. 2020Y9125), Talent programs granted from The First Affiliated Hospital of Fujian Medical University (YJRC3600).
The authors declare no competing interests.
Ethics approval and consent to participate
This study was conducted in accordance with the Declaration of Helsinki and was approved by the Ethics Committee of Fujian Medical University Union Hospital. All patients had written informed consent.
Consent for publication
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Guan, B., Xu, M., Zheng, R. et al. Novel biomarkers to predict treatment response and prognosis in locally advanced rectal cancer undergoing neoadjuvant chemoradiotherapy. BMC Cancer 23, 1099 (2023). https://doi.org/10.1186/s12885-023-11354-8