Next Article in Journal
Isolation of High-Purity Betanin from Red Beet and Elucidation of Its Antioxidant Activity against Peroxynitrite: An in vitro Study
Previous Article in Journal
Physiological Analysis and Genetic Mapping of Short Hypocotyl Trait in Brassica napus L.
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Data-Independent Acquisition Mass Spectrometry Analysis of FFPE Rectal Cancer Samples Offers In-Depth Proteomics Characterization of the Response to Neoadjuvant Chemoradiotherapy

by
Aleksandra Stanojevic
1,
Martina Samiotaki
2,
Vasiliki Lygirou
3,
Mladen Marinkovic
4,5,
Vladimir Nikolic
6,
Suzana Stojanovic-Rundic
4,5,
Radmila Jankovic
1,
Antonia Vlahou
3,
George Panayotou
2,
Remond J. A. Fijneman
7,
Sergi Castellví-Bel
8,9,10,
Jerome Zoidakis
3,11 and
Milena Cavic
1,*
1
Department of Experimental Oncology, Institute for Oncology and Radiology of Serbia, Pasterova 14, 11000 Belgrade, Serbia
2
Institute for Bioinnovation, Biomedical Sciences Research Center “Alexander Fleming”, Fleming 34, 166 72 Vari, Greece
3
Department of Biotechnology, Biomedical Research Foundation, Academy of Athens, 4 Soranou Ephessiou Street, 115 27 Athens, Greece
4
Clinic for Radiation Oncology and Diagnostics, Department of Radiation Oncology, Institute for Oncology and Radiology of Serbia, Pasterova 14, 11000 Belgrade, Serbia
5
Faculty of Medicine, University of Belgrade, Dr Subotica 8, 11000 Belgrade, Serbia
6
Clinic for Medical Oncology, Institute for Oncology and Radiology of Serbia, Pasterova 14, 11000 Belgrade, Serbia
7
Department of Pathology, The Netherlands Cancer Institute, Plesmanlaan 121, 1066 CX Amsterdam, The Netherlands
8
Gastroenterology Department, Fundació Clínic per la Recerca Biomèdica-Institut d’Investigacions Biomèdiques August Pi i Sunyer (FRCB-IDIBAPS), C/del Rosselló, 149, 08036 Barcelona, Spain
9
Centro de Investigación Biomédica en Red de Enfermedades Hepáticas y Digestivas (CIBERehd) Almagro, 3, 28029 Madrid, Spain
10
Hospital Clínic, University of Barcelona, C/del Villarroel, 170, 08036 Barcelona, Spain
11
Department of Biology, National and Kapodistrian University of Athens, Panepistimíou 30, 106 79 Athens, Greece
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2023, 24(20), 15412; https://doi.org/10.3390/ijms242015412
Submission received: 25 August 2023 / Revised: 30 September 2023 / Accepted: 3 October 2023 / Published: 21 October 2023
(This article belongs to the Special Issue Looking Closer to See Bigger: Challenges in Single-Cell Proteomics)

Abstract

:
Locally advanced rectal cancer (LARC) presents a challenge in identifying molecular markers linked to the response to neoadjuvant chemoradiotherapy (nCRT). This study aimed to utilize a sensitive proteomic method, data-independent mass spectrometry (DIA-MS), to extensively analyze the LARC proteome, seeking individuals with favorable initial responses suitable for a watch-and-wait approach. This research addresses the unmet need to understand the response to treatment, potentially guiding personalized strategies for LARC patients. Post-treatment assessment included MRI scans and proctoscopy. This research involved 97 LARC patients treated with intense chemoradiotherapy, comprising radiation and chemotherapy. Out of 97 LARC included in this study, we selected 20 samples with the most different responses to nCRT for proteome profiling (responders vs. non-responders). This proteomic approach shows extensive proteome coverage in LARC samples. The analysis identified a significant number of proteins compared to a prior study. A total of 915 proteins exhibited differential expression between the two groups, with certain signaling pathways associated with response mechanisms, while top candidates had good predictive potential. Proteins encoded by genes SMPDL3A, PCTP, LGMN, SYNJ2, NHLRC3, GLB1, and RAB43 showed high predictive potential of unfavorable treatment outcome, while RPA2, SARNP, PCBP2, SF3B2, HNRNPF, RBBP4, MAGOHB, DUT, ERG28, and BUB3 were good predictive biomarkers of favorable treatment outcome. The identified proteins and related biological processes provide promising insights that could enhance the management and care of LARC patients.

1. Introduction

Colorectal cancer (CRC) is the third most common type of cancer worldwide, with almost two million newly diagnosed cases in 2020 [1]. High mortality rates place CRC in second place after lung cancer [1]. In Serbia, the situation is similar. According to data from the Institute of Public Health “Milan Jovanovic Batut” of the Republic of Serbia, CRC holds second place by incidence and mortality rates, with around 5000 new cases and 3000 deaths annually [2]. In the majority of cases, it is diagnosed in advanced stages where limited treatment options are available and survival is poor. Our group and others have invested efforts into profiling the diagnostic, prognostic, and predictive factors for CRC and anal cancer in an effort to provide better research strategies for treatment and overall management [3,4,5,6,7]. However, current early detection and screening programs, as well as treatment options, need further improvement on a global level. Rectum is the most distal part of the digestive tract, located between the sigmoid colon and anal canal. Colon and rectal cancers have been traditionally considered as a single disease entity, and rectal cancer (RC) represents around 35% of diagnosed CRC cases. Rectal cancer has distinct environmental and genetic risk factors that differentiate it from colon cancer [8]. Its incidence has been reported to increase in the 18- to 50-year-old age group, especially in younger adults [9].
Locally advanced rectal cancer (LARC) is the most diagnosed type of RC, which includes stage II (T3/4N0M0) and III (T1-4 N+ M0), according to the Union for International Cancer Control (UICC) [10]. The standard treatment for LARC is neoadjuvant chemoradiotherapy (nCRT) followed by radical surgery (total mesorectal excision). nCRT was established as the gold standard of LARC treatment after 2004 as a result of two studies, CAO/ARO/AIO-94 and EORTC 22921, which compared it to previously used adjuvant radiotherapy with or without chemotherapy. According to the CAO/ARO/AIO-94, neoadjuvant RT dramatically reduced the rates of local failure, while the 11-year follow-up update showed that the long-term overall survival rate was about 60% [11]. The EORTC 22921 trial showed that the use of chemotherapy with neoadjuvant radiation reduced local recurrence rates but had no effect on distant progression-free survival. Additionally, nCRT contributed to sphincter preservation and improved the patient’s quality of life [12]. However, only 20–30% of patients experience a complete clinical or pathological response to nCRT, while some patients will experience poor response or will have distant progression during nCRT [13,14]. Characterization of mechanisms of response to therapy and the search for predictive biomarkers to nCRT is an unmet need in LARC.
The watch-and-wait approach was introduced because of the need for close follow-up of LARC patients with complete clinical response to nCRT and allowed the extension of periods between neoadjuvant therapy and surgery, thus lowering morbidity related to surgery [15,16,17,18]. No biomarker has yet been validated in this setting. Patients with a favorable response to nCRT would be candidates for a less invasive surgical approach or would be enrolled in a watch-and-wait approach in the case of a complete clinical response (cCR). That would increase their quality of life and contribute to the overall reduction of treatment costs [18].
In this study, we aimed to perform in-depth proteomics characterization of preoperative LARC biopsy samples by employing data-independent acquisition mass spectrometry (DIA-MS) to unravel new tissue molecular features that might lead to different responses to nCRT.

2. Results

2.1. Study Design and Rationale

Neoadjuvant chemoradiotherapy is a standard treatment for locally advanced rectal cancer, but understanding the mechanism of response to therapy is still an unmet need. For that purpose, we examined our cohort of 97 LARC patients whose clinical profile was well characterized. Ninety-seven patients with LARC treated at the Institute for Oncology and Radiology of Serbia from 2018 to 2019 were included. Inclusion criteria comprised histopathologically confirmed LARC (T3-T4N0 or any T stage N+), a distant margin up to 12 cm from the anal verge, and ECOG performance status ≤ 2. Patients received long-course chemoradiotherapy with concurrent chemotherapy. Tumor response was assessed after surgery according to the classification by Mandard, and patients were divided into responders (TRG 1 and 2—complete and near complete response, respectively) and non-responders (TRG 3, 4, and 5—moderate, poor, and no response, respectively) based on pathohistological criteria. Our study involved the analysis of 20 patient samples exhibiting a range of responses to achieve a comprehensive understanding of the diverse molecular features potentially associated with response. The observed results will help us better understand the mechanisms behind the response, enabling the selection of protein candidates that can serve as predictive biomarkers for favorable or unfavorable responses. Patients expected to have a favorable response may be considered for a less invasive surgical approach or enrollment in a watch-and-wait strategy in case of complete clinical response (cCR). This approach can lead to good local control without the morbidity associated with radical surgery after neoadjuvant treatment, subsequently improving the quality of life in this group of patients. Conversely, non-responders may be candidates for intensified neoadjuvant treatment and earlier surgical intervention without delay following initial treatment, the introduction of targeted therapy, or other treatment adoptions.

2.2. Proteomic Comparison of Responders and Non-Responders

The use of DIA-MS allowed the identification and quantification of more than 3000 proteins per sample Figure 1, a significant increase when compared to the 1000 proteins identified by DDA-MS in rectal cancer FFPE samples [19]. In total, 4269 proteins were identified in 20 rectal cancer FFPE samples (Table S1), while 1923 proteins were identified in all samples. After 50% filtering and log2 transformation, the Welch t-test was applied. The number of identified proteins before and after 50% is shown in Figure 1. Statistical analysis indicated 915 DEPs with significant differences (p < 0.05; S0 = 0.1) between responders and non-responders. When two groups were compared, DEPs included 700 proteins overexpressed in non-responders and 215 overexpressed in responders (Table S2 and Figure 2b). Using a more stringent statistical setting (p < 0.01; S0 = 0.1), 384 DEPs were found between responders and non-responders, 81 of which were upregulated in responders compared to non-responders, and 303 proteins with the opposite trend (Table S3).
Principal Component Analysis (PCA) (Figure 2a) shows the separation of the patients depending on their proteomic profile, indicating that responders had significantly different proteomics profiles than non-responders. On the PCA plot, we observe two samples that are outliers from the group of non-respondents (67 and 70) as well as sample 65 from the group of respondents. A potential reason for this behavior is the fact that these samples have a lower number of identified proteins compared to others, which may be associated with the lower quality of the sample and divergent proteomic profile. For comparison, a volcano plot (−log10(p-value) vs. Welch t-test difference) was created to graphically show the proteome changes between the two groups of samples (Figure 2b).
Hierarchical clustering analysis was performed on Z-score normalized data with stricter Welch’s t-test statistics (p < 0.01), while distance for analysis was performed using Euclidean distance (Figure 3). Hierarchical clustering analysis (Figure 3) revealed details on protein abundance among two groups of responders and non-responders; clustering also provided clear grouping results depending on the response, with samples 67 and 70 as an outlier. This sample had the lowest amount of detected proteins, which can be a possible reason for aberrant classification.
Furthermore, a hierarchical clustering analysis was performed for the top 10 proteins from the group of responders and non-responders based on the strength of statistical significance (Figure 4). Based on the obtained results, we conclude that the tissue proteomic profile differs significantly between patients and thus enables their clear classification. Proteins encoded by genes SMPDL3A, PCTP, LGMN, SYNJ2, NHLRC3, GLB1, and RAB43 showed high predictive potential of unfavorable treatment outcome, while RPA2, SARNP, PCBP2, SF3B2, HNRNPF, RBB4, MAGOHB, DUT, ERG28, and BUB3 were good predictive biomarker of favorable treatment outcome. Hemoglobin-related proteins (HBB, HBA1, and HBD) can be possible contaminants originating from FFPE tissue samples. The mentioned proteins are good candidates with great potential for validation on a larger group of patients, which would confirm their discriminatory potential.

2.3. Pathway Enrichment Analysis

Enrichment pathway analysis was performed on proteins that were significantly differentially expressed (p < 0.05; S0 = 0.1). Initially, all proteins were chosen regardless of grouping with the goal of better understanding the signaling pathways involved in treatment response (Figure 5a). Further enrichment analysis was carried out on two groups (responders and non-responders) independently in an attempt to explain discrepancies in treatment response. The findings were analyzed and represented based on their biological significance to RC biology. To keep the analysis output concise, only the leading terms of each pathway are shown. Results indicated that some of the leading signaling pathways that correlate with response to nCRT in patients with LARC include the metabolism of RNA, MYC targets, neutrophil degranulation, cellular transport, and response to stimuli. The responder group was characterized by signaling pathways related to cell cycle signaling (metabolism of RNA, synthesis of DNA, DNA strand elongation, mitochondrial translation initiation, chromosome maintenance), as well as MYC targets scores, regulation of expression of SLITs and ROBOs, mTOR1 signaling pathway, and unfolded protein response (Figure 5c). The non-responder group was characterized by signaling pathways related to vesicle-mediated transport, neutrophil degranulation, hemostasis, coagulation, heme metabolism, post-translational modifications, as well as the metabolism of vitamins, cofactors, and lipids. Signaling pathways related to epithelial–mesenchymal transition and hypoxia, which have been associated with an increased risk of metastasis, were also found to be important in non-responders (Figure 5b).

2.4. STRING In Silico Analysis

Data obtained reveal several protein-rich groups with several members having high levels of interactions in the responder and non-responder group (PPI enrichment p-value: <1.0 × 10−16), indicating that proteins interact with one another more frequently than would be predicted by a randomly selected group of proteins from the genome with the same size and degree distribution [20]. After Cytoscape analysis using MCODE extension, six clusters were detected. We can conclude that there is a strong interaction between proteins involved in information RNA processing and genes whose protein products participate in transcription, especially when it comes to pre-mRNA processing (Figure S1a) and factors involved in the process of alternative splicing (Figure S1b). A high degree of interaction is also associated with proteins that participate in the formation of snRNA molecules (Figure S1b). Another group of proteins that are clearly distinguished based on STRING analysis are the ribosomal proteins of the RPL family (Figure S1e) and MRPL family (Figure S1f), as well as PA2G4, which represents the proliferation factor (Figure S1e). All the mentioned groups of proteins are characterized by a high mutual connection. Another group of proteins includes factors involved in the DNA replication process as well as factors for the organization of the proteasomal system (Figure S1c,d). On the other hand, the group of DEPs overrepresented in the non-responder group is characterized by a much larger number of proteins that are less closely related to each other. One of the detected clusters included proteins that control cell death, proliferation, and signal transduction (Figure S2a), as well as proteins involved in lipid metabolism (Figure S2b). Cluster 3 included proteins that are related to the mitochondrial electron transport chain (Figure S2c), while cluster 4 included proteins related to retrograde electron transport (Figure S2d), followed by purine and pyrimidine catabolism (Figure S2d) and the HLA class of proteins included in antigen processing and presentation (Figure S2e).

2.5. Shortlisting of Potential Biomarkers Based on Transcriptomics Data

The proteomics results obtained were further examined to discover promising predictive biomarkers of response to neoadjuvant chemoradiotherapy in patients with LARC. The differential expression of proteins identified in our study was confirmed in transcriptomics datasets. The ROC curve was considered to assess the performance of predictive biomarkers for response to chemoradiotherapy. For this purpose, DEPs obtained after DIA-MS/MS were analyzed using ROCplotter software (https://www.rocplot.org, accessed on 1 May 2022), and genes with AUC > 0.7, ROC p-value < 0.05, and Mann–Whitney p-value < 0.05 were categorized as promising biomarkers. Out of a total of 915 DEPs, 23 genes met all three criteria. The responder group had the following proteins whose expression was also confirmed at the mRNA level: CRKL, LAP3, THTPA, PES1, PPP2R5E, IFI30, C17orf75, QDPR, RRM2B, USO1, GLRX ARAF, CTBS, and SNRPD3. Moreover, the non-responder group had the following proteins, COPB1, MGLL, HAS1, TALDO1, DNAH9, KDELR3, HLA-DPB1, RBP3 and STAP2, as presented in Table 1 and Table 2. ROC curves and their discriminatory potential for the mentioned genes are shown in Figure 6 and Figure 7. Linking proteomics with transcriptomics data can lead to the discovery of promising rectal cancer biomarkers that are easy, cost-effective, and fast to detect.

2.6. Search for Drug Targets

In order to search for drug targets among proteins that are differentially overexpressed in the group of non-responders, GeneCards and DrugBank databases were searched. Some of the proteins that are differentially overexpressed in the group of non-responders versus responders are drug targets used in the treatment of some other pathological conditions including QPRT, CLCA4, ATG4B, and PTGS2 (Table 3) [22,23]. The question of the discriminatory effect of gene expression is raised, as well as whether the use of these drugs can be used as a part of initial treatment and would lead to a better response to therapy by treating patients with locally advanced rectal cancer together with standard neoadjuvant chemoradiotherapy.

3. Discussion

Understanding the molecular features associated with response to neoadjuvant chemoradiotherapy is an unmet clinical need in LARC. The DIA-MS offers unprecedented proteome coverage for FFPE samples. DIA-MS enabled the in-depth study of the proteome from FFPE tissue samples, which represented a major challenge because of damage caused by the fixation protocol [24]. By exploring the dynamic phenotypic characteristics of tumor cells before therapy and tumor response to therapy, DIA-MS allows us to characterize response mechanisms and thus enable patient monitoring and more effective treatment. FFPE samples are routinely used for DDA-MS analysis, and the use of FFPE samples for the DIA-MS method is increasing [25]. In contrast to fresh frozen (FF) tissue, FFPE samples undergo protein cross-linking during standard preservation protocol, and due to that, it is challenging to analyze native proteins. Comparing the results obtained when it comes to FF tissue versus FFPE samples indicates a high correlation between the results, which makes FFPE samples a good alternative to FF samples [26,27]. On the other side, FFPE samples are more suitable for retrospective studies and are easily accessed during everyday clinical routines [25]. Due to technical limitations, results obtained in the previous studies showed a restricted number of identified proteins, while DIA-MS/MS offered an in-depth characterization of rectal cancer tissue, enabling molecular characterization and profiling of the response [19,28,29]. Our study included analysis of discrete and well-characterized clinical samples of rectal cancer in order to identify the maximum number of different molecular features potentially associated with response.
In total, 4269 proteins were identified in 20 rectal cancer FFPE samples. Principal Component Analysis (PCA) indicated that responders had a significantly different proteomics profile than non-responders. Statistical analyses comparing the two groups resulted in the identification of 915 differentially expressed proteins (215 in responders and 700 in non-responders) (p < 0.05).
The therapy approach used for LARC includes radiotherapy in combination with 5-fluorouracil-based chemotherapy followed by surgery. Radiotherapy primarily exerts its effects by damaging DNA through the generation of molecular fragments like free radicals and excited molecules. Cells are most sensitive to radiation during mitosis and the early G1 phase, while they are most resistant during the S phase of the cell cycle [30]. The abovementioned effects are first reflected in the cell cycle, and the overcoming of toxic effects can be carried out through the reparation of the resulting damage. 5-fluorouracil is a standard chemotherapeutic drug (antimetabolite) that is metabolized by the liver. Within cells, 5-FU undergoes metabolic conversion, resulting in three active metabolites: 5-FdUMP and 5-FdUTP, which damage DNA, and 5-FUTP, which integrates into RNA, exerting an antiproliferative effect. Our finding confirmed the significance of DNA metabolism and related signaling pathways in treatment response. Patients with a poor response exhibit deregulated pathways related to DNA strand elongation and synthesis.
Based on hierarchical clustering analysis, we conclude that the tissue proteomic profile differs significantly between patients, which enables their clear classification (Figure 4). Additionally, some of those proteins have been shown to be included in the development and progress of many other types of cancers, but none of them were investigated in terms of rectal cancer therapy resistance. All of the top 10 proteins from both groups showed high predictive potential. It has been shown that protein legumain (LGMN) is overexpressed in breast, prostate, and liver cancer and that its role is significant in cancer development, progression, and invasion [31,32]. This protein was also shown to be significantly associated with the development of peritoneal metastases. In our study, this protein was shown to be a promising predictive biomarker of poor response to therapy. This protein is a member of several pathways shown to be significant in response to nCRT, including mTOR, coagulation, adaptive immune system, and lipid metabolism. Pathway mTORC1 was altered in the group of good responders. Apart from LGMN, there is BUB3, which is shown to be overexpressed in the responder group and included in the mTOR signaling pathway as well. Analysis of DEPs provided a potential scenario that included the downregulation of genes related to the mTORC1 pathway in responders or overexpression of mTORC1-related genes in non-responders. Considering available data, both scenarios will lead to poor response to treatment. It has been demonstrated that treatment resistance in a variety of cancer types is correlated with the stimulation of mTOR signaling pathways [33,34]. Glycolysis, glycoprotein and lipid synthesis, mitochondrial and lysosome biogenesis, and metabolic balance all depend on mTORC1. Translation is directly impacted by the regulation of numerous transcriptional factors. It regulates the production of nucleotides and the metabolism of glucose through proteins in the metabolic pathway. Additionally, mTORC1 regulates the assembly of the proteasome and autophagy [35]. Glucose metabolism, apoptosis, cell migration, cytoskeletal reorganization, and cell proliferation are regulated by mTORC2 [36,37]. According to our results, signaling pathways related to protein and lipid synthesis, glycolysis, mitochondrial biogenesis, and lysosome biogenesis correlate with poor response to nCRT, which is in compliance with previous findings. An in-depth characterization of mTORC1 signaling pathways is needed to shed light on the exact biochemical mechanism that leads to good/poor responses to treatment.
Our results indicated that the immune response might play an important role in predicting response to therapy. Signaling pathways associated with IL-12 were found to correlate with a good response to therapy, while signaling pathways related to the adaptive immune system were related to a poor response. Heeran et al. highlighted the association of inflammation with obesity status in rectal cancer patients in terms of lowering the level of inflammatory factors released from TME [38]. Our data correlated adipogenesis with poor response treatment. As increased adipose tissue synthesis is directly associated with a rise in body weight, it could also result in reduced therapeutic response. A study by Lee et al. has shown that obesity represents an independent predictor of cCR, which contradicts our results [39]. In some studies, no clear correlation was found between obesity and rectal cancer treatment outcomes [40]. Next to it, Synaptojanin 2 (SYNJ2) was shown to be a good marker of poor prognosis in lung cancer. This protein is included in signaling pathways altered in the non-responder group, including vesicle-mediated transport and lipids. When it comes to predictive markers of a good response, none of the presented proteins were investigated in terms of predictive potential. It was shown that RB-binding protein 4 (RBBP4) is associated with a poor prognosis of colon cancer, while our results suggest a favorable effect. This protein is part of two signaling pathways that are alerted in good response. Slit–Robo signaling plays an important role in angiogenesis. The vertebrate Robo4 gene, which has been associated with regulating angiogenesis and blood artery permeability, has highly specific endothelial cell expression [41,42]. In our cohort, it was shown that the regulation of SLITs and ROBOs is highly correlated with a good response to nCRT. Slit/Robo signaling has both pro- and anti-angiogenic functions; therefore, its effect on angiogenesis depends on the environment. It has been demonstrated that Slit2 stimulates angiogenesis through the Robo1 receptor but that it also inhibits endothelial migration through the Robo4 receptor, which aids in a positive therapeutic response. It was also reported that Split3 promotes angiogenesis [43]. Apart from SYNJ2, MAGOHB and DNA-directed RNA polymerase I subunit RPA2 were shown to be overexpressed in responder group. In terms of good response to treatment and favorable outcomes, inhibition of endothelial migration and angiogenesis is preferred. Based on the results of research previously conducted by our group, we found that hematological parameters, including neutrophil-to-monocyte ratio, initial basophil, eosinophil, and monocyte counts, are significantly different between the responder and non-responder groups. According to MRI findings, we realized that non-responders are more often presented with extramural vascular invasion [44]. In support of all of those observations, our results indicated the importance of epithelial–mesenchymal transition in patients with a poor response to therapy. Hypoxia has also been linked with neutrophil degranulation. In hypoxic conditions, degranulation occurs, and released factors affect tumor progression [45]. According to our results, the neutrophil degranulation pathway was affected in both the responder and non-responder groups, which highlights its importance in this process and warrants further functional studies. And some proteins of this pathway that showed good predictive potential are coded by gene HNRNPF in the responder group and NHL repeat-containing protein 3 (NHLRC3) and Beta-galactosidase (GLB1) in the non-responder group.
In rectal cancer patients receiving neoadjuvant CRT, it was shown that high expression of FGFR2 was associated with an advanced tumor stage, a poor treatment response, and lower survival [46]. The DIA-MS approach indicated the importance of FGFR2 alternative splicing in good response to treatment, and the protein encoded by gene HNRNPF was shown to have good predictive potential, which implies that exploring its variants might be useful for the prediction of a good outcome.
The highest correlation with poor response to treatment was shown for signaling pathways related to vesicle-mediated transport and endocytosis, while proteins encoded by genes LGMN, GLB1, and SYNJ2 and Ras-related protein Rab-43 (RAB43) are part of listed signaling pathways and good predictors of poor response. Proteins and other cargo must be carried through the cell via a cellular transport mechanism in which the transported materials are conveyed in membrane-bound vesicles. The vesicle lumen or the vesicle membrane is where transported compounds are contained [47]. By carrying biomolecules (proteins, lipids, deoxyribonucleic acid, and ribonucleic acid) throughout the tumor microenvironment, exosomes released from cell membrane play a crucial role in tumor proliferation, differentiation, metastasis, and resistance to chemotherapy and radiation [48]. The results obtained in this study indicated a great potential for exploring intercellular communication in the tumor microenvironment as well as in the tumor when profiling response to therapy. The synergistic effect of these inter-relations has not yet been clarified and would also be validated by our group in future functional studies.
Deoxyuridine triphosphatase (DUT) is associated with a shorter DFS in patients with CRC [49]. This protein is important for DNA strand elongation and unfolded protein response and is an MYC target. In addition to DUT, SF3B2 and BUB3 were shown to be overexpressed in responder group. Signaling pathways related to MYC targets have also been found to be significant for the good response to nCRT in our LARC cohort. MYC is an oncogene and transcription factor that regulates cell-cycle-related signaling pathways, supporting their crucial effect in response to treatment [50]. Signaling pathways, which include chromosome maintenance and telomere extension by telomerase, were found to be affected in good response. This observation suggested that changes at the chromosome level can lead to a good response to treatment and increase the efficiency of RT. Short telomeres can lead to chromosomal instability and the formation of cancer, while on the other hand, long telomeres, due to their length, can undergo a higher number of divisions and thus increase the probability of transformation of a normal cell into a malignant cell, this phenomenon called telomere length paradox [51]. When it comes to the response to therapy, the correlation between telomere length and cancer treatment is not fully clarified [52]. Additionally, mTOR, a key regulator of cell growth and division in healthy conditions, can be inappropriately activated in tumor cells and thus promote tumor cell growth, metastasis, and invasion of fresh, healthy tissues [53]. Cytokines and growth factors are released during immunological response against harm that CRT may cause, and they play a big part in the generation of ROS, including superoxide, hydrogen peroxide, and nitrogen (II)-oxide [54]. Interleukins (IL-2, IL-12) have been found in preclinical studies that might help to modulate the antitumor response and radiosensitize cells [55]. IL-12 achieves its antitumor activity by promoting the immune response via the activation of natural killer cells and cytotoxic T cells and exerting an anti-angiogenic effect [56,57]. In a study performed by Heeran et al., increased levels of IL-12 were detected in the blood of LARC patients compared to healthy individuals. This suggests the potential for promoting an immune response that may lead to improved treatment outcomes [38]. Our study showed that the protein encoded by the gene HNRNPF has predictive potential for good treatment outcomes.
The data in the literature from the past several decades identified altered glycosylation as a sign of malignancy. It was shown that glycosylation acts as a mediator of the inflammatory response [58]. In gliomas and laryngeal carcinoma, a correlation between glycosylation and radioresistance has already been demonstrated [59]. According to our results, glycan degradation and asparagine N-linked glycosylation might play an important role in the poor response of LARC patients to nCRT, proteins encoded by genes PCTP and GLB1 shown to have high potential in the prediction of poor response to treatment.
Linking the transcriptomic and proteomic profiles of the cell is an important parameter for understanding the molecular basis of the response to therapies. The software ROC plotter (www.rocplot.org, accessed on 3 May 2022) showed a low correlation between gene expression and the proteome profile of the tested samples. Out of 915 differentially expressed proteins, only 23 showed promising discriminatory potential when it comes to gene expression. It should be kept in mind that the transcriptomics data were obtained from samples of different demographics and ethnic origins; therefore, further validation on patient samples within our cohort should be performed. Genes that have shown a transcriptional discriminatory potential for predicting a good response have also not been investigated in rectal carcinoma so far. Published results indicate that PPP2R5E, PES1, RRM2B, GLRX, and CRKL are important in the prognosis and development of colorectal cancer [60,61,62,63,64,65]. When it comes to the prediction of poor response, some of the aforementioned genes showed a poor prognostic potential risk for the development of CRC. The HLA-DPB1 gene has been shown to be a good predictive marker of response to nCRT in patients with rectal cancer [66]. These results indicate that there is great potential in examining the level of expression of the mentioned genes. Exploiting this potential would lead to rapid, inexpensive, and easy methods for predicting response to therapy in patients with LARC. The translation of research from protein to gene expression and confirmation of the obtained candidates would enable a more cost-effective approach and, thus, a more efficient selection of patients when it comes to predicting the response to neoadjuvant chemoradiotherapy. Gene expression analysis shows high sensitivity when it comes to FFPE tissue analysis and would be more likely to be performed during everyday practice. At the same time, further analysis of the gene expression profile and correlation with a proteomic profile of the tissue would enable a more detailed investigation of the mechanism behind response treatment, which is still unclear.
Based on all obtained results, we conclude that there is a statistically significant difference between the proteomic profiles of LARC patients who respond well and those who respond poorly to nCRT. As all patients require surgery after nCRT per current guidelines, profiling of adequate biomarkers of response is a pressing matter. Further validation of target signaling pathways detected in this study that might have an effect on the response to nCRT is planned on a prospective cohort of LARC patients at the Institute for Oncology and Radiology of Serbia to ensure more efficient and cost-effective treatment of patients while maximizing their quality of life.

4. Materials and Methods

4.1. Patient Cohort Characteristics and Treatment

A total of 97 LARC patients treated at the Institute for Oncology and Radiology of Serbia in the period of 2018–2019 were included in this study. The inclusion criteria were histopathologically verified adenocarcinoma of the rectum, with a distant margin up to 12 cm from the anal verge by rigid proctoscopy, ECOG performance status ≤ 2. LARC was defined as T3–T4N0 or any T stage N+, according to clinical and histological criteria of the 8th edition of the TNM classification of malignant tumors [67]. Pretreatment evaluation included an abdominal and pelvic MRI scan and a computed tomography (CT) scan or X-ray of the chest. All patients were treated with long-course chemoradiotherapy (CRT). Radiotherapy (RT) was delivered with a total dose of 50.4 Gy in 28 fractions (conventionally fractioned 1.8 Gy/fr), using the technique with 3 or 4 radiation areas (all areas as recommended by the International Committee of Radiation Units and Measurements (ICRU, 50/62 per day) [68,69]. Concomitant chemotherapy was initiated on the first day of RT and administered during the first and fifth weeks of RT. The chemotherapy regimen included 5-FU (350 mg/m2 daily) and Leucovorin (25 mg/m2 daily). A complete patient medical database has been prepared from official records.
Patients were assessed for tumor response between the 6th and 8th weeks after CRT completion with pelvic MRI scan and rigid proctoscopy. The pathohistological response after surgery was assessed according to tumor regression grading (TRG) categories by Mandard [70]. According to the TRG status, the patients were divided into two groups: responders (TRG 1–2) and non-responders (TRG 3–5). Our study included analysis of extreme candidates in order to achieve the maximum range of different molecular features potentially associated with response. Twenty-four formalin-fixed paraffin-embedded (FFPE) biopsy samples were taken at the moment of disease diagnosis and were collected and used for proteomic analysis. After the quality control check, four samples were excluded from further analysis, and finally, twenty samples were processed (9 responders and 11 non-responders). Characteristics of the study cohort are shown in Table 4.
The project was approved by the Ethics Committee of the Institute for Oncology and Radiology of Serbia (approval No. 2211-01 from 11 June 2020), and all patients signed informed consent. All experiments have been performed in accordance with the Helsinki Declaration of 1975, as revised in 2013.

Protein Extraction from FFPE Tissue Samples

From each of the 20 FFPE LARC samples (9 responders and 11 non-responders), 10 sections 15–20 μm thick were cut using a microtome and transferred to 2 mL Eppendorf tubes. From all sections, 3 were selected that contained the largest amount of tissue. Samples were deparaffinized using Xylene in two steps of 5 min and 1 min successively. The samples were centrifuged, and in the next step, the pellet tissue was rehydrated using several dilutions of ethanol and finally washed with double distilled water. The tissue pellet was air-dried and dissolved in FASP protein extraction buffer (100 mM Tris-HCl pH 7.6, SDS 4%, 100 mM 1,4-Dithioerythtriol (DTE)). The tissue was homogenized, while the disintegration of the cell membrane was achieved by sonication (three cycles of 5 s with 36% power). The samples were heated at 90 °C for 1 h to extract the cellular proteins into the solution. The supernatant containing extracted proteins was transferred to a new Eppendorf tube, and an appropriate amount of ammonium bicarbonate (ABC) buffer was added. The sample was concentrated using a 3 kDa cut-off Amicon filter.

4.2. Protein Digestion and Preparation for LC-MS/MS Analysis

The total volume of concentrated proteins was added on SDS PAGE (5% stacking gel, 12% separating gel). Preparative SDS PAGE was performed, and the gel was fixed, washed, and stained with Coomassie colloidal dye (File S1). Protein bands were cut from the gel for each sample separately, chopped, and transferred to Eppendorf tubes. The strips were decolorized with a solution of 40% Acetonitrile and 50 mM NH4HCO3 until parts of the gel became completely transparent. The samples were reduced with 10 mM DTE in 100 mM NH4HCO3 and alkylated 10 mg/mL Iodoacetamide in 100 mM NH4HCO3 and then washed with 100 mM NH4HCO3, destain solution, and water, respectively. The samples were dried in a speed vac until transparent crystals formed. Each sample was treated with trypsin solution, which enables the cutting of the polypeptides after the amino acids lysine and arginine. The peptides formed were extracted with NH4HCO3 solution followed by incubation in a 1:1 solution of 10% formic acid and Acetonitrile. The peptide solution was purified using PVDF filters (Merck Millipore, Darmstadt, Germany). The samples were dried in a Speedvac and prepared for further processing.

4.3. LC-MS/MS Analysis

Samples were run in two technical replicates on a liquid chromatography–tandem mass spectrometry (LC-MS/MS) setup consisting of an Ultimate 3000 RSLC online with a Thermo Q Exactive HF-X Orbitrap mass spectrometer. Peptide solutions were directly injected and separated on a 25 cm long analytical C18 column (PepSep, 1.9 μm beads, 75 µm ID) using a gradient of 7% to 36% Buffer B (0.1% formic acid in 80% Acetonitrile) for 70 min, followed by an increase to 95% in 5 min, a second increase to 99% in 0.5 min, and then kept constant for equilibration for 14.5 min. A full MS was acquired in profile mode using a Q Exactive HF-X Hybrid Quadropole-Orbitrap mass spectrometer, operating in the scan range of 375–1400 m/z using 120 K resolving power with an AGC of 3 × 106 and max IT of 60 ms. Data-independent analysis followed, using 8 Th windows (39 loop counts) with 15 K resolving power with an AGC of 3 × 105, max IT of 22 ms, and a normalized collision energy (NCE) of 26. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE [71,72] partner repository with the dataset identifier PXD040451.

4.4. MS Data Analysis

Orbitrap raw data were analyzed in DIA-NN 1.8 (Data-Independent Acquisition by Neural Networks) by searching against the reviewed Human Uniprot database (retrieved 4/21) in the library free mode of the software, allowing up to two tryptic missed cleavages. Human Uniprot Database includes 27,246 proteins and 21,442 genes with 10,241,864 precursors generated. A spectral library was created from the DIA runs and used to reanalyze them. Parameters regarding peak generation and analysis are defined in the DIA-NN algorithm. DIA-NN default settings have been used with oxidation of methionine residues and acetylation of the protein N-termini set as variable modifications and carbamidomethylation of cysteine residues as fixed modification. N-terminal methionine excision was also enabled. A maximum number of variable modifications is set to 3. Both ends are fully tryptic, allowing up to two tryptic missed cleavages. The match between runs feature was used for all analyses, and the output (precursor) was filtered at 0.01 FDR. Retention time alignment is performed in DIA-NN. Correction for the mass accuracy is performed for each sample in DIA-NN automatically. Filtering of the quality is based on the false discovery rate of 0.01 at peptide and protein levels. Finally, the protein inference was performed on the level of genes using only proteotypic peptides. The analysis was set for at least one unique peptide per protein. The generated results were processed statistically and visualized in the Perseus software v1.6.15.0 (Max Planck Institute of Biochemistry, Munich, Germany) and GraphPad Prism 8.0.1. With the help of DIA-MS and DIA-NN processing, we identified 5756 groups of proteins, or 4875 unique proteins, that were used for subsequent analysis. Raw data were filtered based on a minimum of 50% valid values in at least one group of responder/non-responder and log2 transformed. After initial processing, non-human genes were excluded from further consideration, while missing values were replaced with imputed values that correspond to the limit of detection (LOD). The modified Student’s t-test, known as Welch’s Test for Unequal Variances, is used. In general, for samples with unequal variance, the adjusted degrees of freedom tend to increase the test power. Differentially expressed proteins (DEPs) were classified as proteins with p < 0.05 (Unequal Welch t-test with S0 cut off 0.1). Proteins with a Welch t-difference above 0 were classified as overexpressed in responders compared to non-responders, while proteins with a Welch t-test difference lower than 0 were classified as overexpressed in non-responders compared to responders. Visualization of the obtained results was performed using PCA, volcano plot, and hierarchical clustering. PCA was plotted by https://www.bioinformatics.com.cn/en accessed on 15 September 2023, a free online platform for data analysis and visualization, while volcano plots were plotted using https://huygens.science.uva.nl/VolcaNoseR/ accessed on 15 September 2023.

4.5. Pathway Enrichment Analysis

To understand the mechanism of response to treatment, pathway enrichment analysis was performed on DEPs between responders and non-responders using Metascape software (MSBio v3.5.20220422). Enrichment analysis parameters were set on a minimum of three genes overlapping between pathways and the input lists. Kyoto Encyclopedia of Genes and Genomes (KEGG) and Reactome and Hallmark (MSigDB) ontologies were used for correlation. Only statistically significant pathways (p-value ≤ 0.05 and minimum enrichment score above 1.5) were taken into account. The obtained results were considered and represented based on biological relevance with respect to RC biology. As a result, the leading term from each group is provided for simplification.

4.6. STRING In Silico Analysis

The STRING analysis network of DEPs overrepresented in responder/non-responder group was built based on the highest confidence (0.9) evidence from experimental interaction data, co-expression data, gene fusions, gene co-occurrence, gene neighborhood, and predictive and knowledge text mining. For easier data processing of 700 DEPs in non-responder group, disconnected nodes in the network were not presented. The analysis was performed using STRING v.11.5 and corresponding images and data downloaded in the original form with statistical significance set at p < 0.05 [20].

4.7. Shortlisting of Potential Biomarkers

Shortlisting of potential biomarkers was performed using the ROCplotter (www.rocplot.org, accessed on 3 May 2022), an online tool that uses the transcriptome data of a large set of rectal cancer patients (N = 284) to find gene expression-based predictive biomarkers. A single database was created by combining published gene expression data from accessible datasets with treatment information. Receiver operating characteristic (ROC) curve analysis was performed to assess the predictive accuracy of each gene [73,74]. The observed cohort included 42 patients treated with 5-fluorouracil and radiotherapy and was categorized into responders (N = 23) and non-responders (N = 19) according to the Response Evaluation Criteria in Solid Tumors (RECIST) criteria. Using a score method devised to assess each probe set for specificity, coverage, and degradation resistance, the optimal microarray probe set to represent a gene was chosen using the JetSet tool [75]. ROC curve with p-value < 0.05 was considered to evaluate the prediction ability of genes that showed a significant difference between the two groups.

5. Conclusions

Data-independent acquisition mass spectrometry (DIA-MS) analysis of FFPE LARC samples offered unprecedented in-depth proteomics characterization. Proteins encoded by genes SMPDL3A, PCTP, LGMN, SYNJ2, NHLRC3, GLB1, and RAB43 showed high predictive potential of unfavorable treatment outcome, while RPA2, SARNP, PCBP2, SF3B2, HNRNPF, RBBP4, MAGOHB, DUT, ERG28, and BUB3 were good predictive biomarkers of favorable treatment outcome. These newly identified molecular features associated with response to nCRT might prove useful for the construction of predictive panels to improve the management and care of LARC patients.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/ijms242015412/s1.

Author Contributions

Conceptualization, M.C., A.S., V.N., M.M. and S.S.-R.; methodology, A.S., J.Z., V.L. and M.S.; software, A.S., M.S. and V.L.; formal analysis, A.S., V.L., M.S., J.Z. and M.C.; investigation, M.C., A.S., V.N. and M.M.; resources, M.S., G.P., J.Z. and A.V.; data curation, M.S., A.S. and V.L.; writing—original draft preparation, A.S., M.C. and M.S.; writing—review and editing, M.C., M.M., R.J.A.F., S.C.-B., J.Z. and M.S.; visualization, A.S. and V.L.; supervision, M.C., R.J., S.S.-R., A.V., G.P. and J.Z.; project administration, M.C. and J.Z.; funding acquisition, M.C., R.J.A.F., S.C.-B. and J.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This study was funded by the Horizon Europe STEPUPIORS Project (HORIZON-WIDERA-2021-ACCESS-03, European Commission, Agreement No. 101079217) and the Ministry of Education and Science of the Republic of Serbia (Agreement No. 451-03-47/2023-01/200043). This article is based upon work from COST Action CA17118, supported by COST (European Cooperation in Science and Technology).

Institutional Review Board Statement

The project was approved by the Ethics Committee of the Institute for Oncology and Radiology of Serbia (approval No. 2211-01 from 11 June 2020), and all patients signed informed consent. All experiments have been performed in accordance with the Helsinki Declaration of 1975, as revised in 2013.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository (https://www.ebi.ac.uk/pride/) with the dataset identifier PXD040451.

Acknowledgments

We acknowledge the support of this work by the project “The Greek Research Infrastructure for Personalised Medicine (pMED-GR)” (MIS 5002802), which is implemented under the Action “Reinforcement of the Research and Innovation Infrastructure”, funded by the Operational Programme “Competitiveness, Entrepreneurship, and Innovation” (NSRF 2014-2020) and co-financed by Greece and the European Union (European Regional Development Fund).

Conflicts of Interest

The authors declare that they have no conflict of interest.

References

  1. Sung, H.; Ferlay, J.; Siegel, R.L.; 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, 209–249. [Google Scholar] [PubMed]
  2. Institute of Public Health of Serbia “Dr Milan Jovanović Batut”. Malignant Tumours in Republic of Serbia; Institute of Public Health of Serbia “Dr Milan Jovanović Batut”: Belgrade, Serbia, 2022. [Google Scholar]
  3. Jakovljevic, K.; Malisic, E.; Cavic, M.; Krivokuca, A.; Dobricic, J.; Jankovic, R. KRAS and BRAF mutations in Serbian patients with colorectal cancer. J. BUON 2012, 17, 575–580. [Google Scholar] [PubMed]
  4. Brotto, K.; Malisic, E.; Cavic, M.; Krivokuca, A.; Jankovic, R. The Usability of Allele-Specific PCR and Reverse-Hybridization Assays for KRAS Genotyping in Serbian Colorectal Cancer Patients. Dig. Dis. Sci. 2013, 58, 998–1003. [Google Scholar] [PubMed]
  5. Čavić, M.; Krivokuća, A.; Boljević, I.; Brotto, K.; Jovanović, K.; Tanić, M.; Filipović, L.; Zec, M.; Mališić, E.; Janković, R.; et al. Pharmacogenetics in cancer therapy–8 years of experience at the Institute for Oncology and Radiology of Serbia. J. BUON 2016, 21, 1287–1295. [Google Scholar] [PubMed]
  6. Stojanovic-Rundic, S.; Marinkovic, M.; Cavic, M.; Karapandzic, V.P.; Gavrilovic, D.; Jankovic, R.; de Voer, R.M.; Castellvi-Bel, S.; Krivokapić, Z. The role of haematological parameters in predicting the response to radical chemoradiotherapy in patients with anal squamous cell cancer. Radiol. Oncol. 2021, 55, 449–458. [Google Scholar] [PubMed]
  7. Nikolic, N.; Radosavljevic, D.; Gavrilovic, D.; Nikolic, V.; Stanic, N.; Spasic, J.; Cacev, T.; Castellvi-Bel, S.; Cavic, M.; Jankovic, G. Prognostic Factors for Post-Recurrence Survival in Stage II and III Colorectal Carcinoma Patients. Medicina 2021, 57, 1108. [Google Scholar]
  8. Glynne-Jones, R.; Wyrwicz, L.; Tiret, E.; Brown, G.; Rödel, C.; Cervantes, A.; Arnold, D.; ESMO Guidelines Committee. Rectal cancer: ESMO Clinical Practice Guidelines for diagnosis, treatment and follow-up. Ann. Oncol. 2017, 28, iv22–iv40. [Google Scholar]
  9. Stoffel, E.M.; Murphy, C.C. Epidemiology and Mechanisms of the Increasing Incidence of Colon and Rectal Cancers in Young Adults. Gastroenterology 2020, 158, 341–353. [Google Scholar]
  10. Brierley, D.J.; Gospodarowicz, K.M.; Wittekind, C. TNM Classification of Malignant Tumours, 8th ed.; Wiley: New York, NY, USA, 2017. [Google Scholar]
  11. Sauer, R.; Liersch, T.; Merkel, S.; Fietkau, R.; Hohenberger, W.; Hess, C.; Becker, H.; Raab, H.R.; Villanueva, M.T.; Witzigmann, H.; et al. Preoperative versus postoperative chemoradiotherapy for locally advanced rectal cancer: Results of the German CAO/ARO/AIO-94 randomized phase III trial after a median follow-up of 11 years. J. Clin. Oncol. 2012, 30, 1926–1933. [Google Scholar]
  12. Oronsky, B.; Reid, T.; Larson, C.; Knox, S.J. Locally advanced rectal cancer: The past, present, and future. Semin. Oncol. 2020, 47, 85–92. [Google Scholar]
  13. Petrelli, F.; Trevisan, F.; Cabiddu, M.; Sgroi, G.; Bruschieri, L.; Rausa, E.; Ghidini, M.; Turati, L. Total Neoadjuvant Therapy in Rectal Cancer: A Systematic Review and Meta-analysis of Treatment Outcomes. Ann. Surg. 2020, 271, 440–448. [Google Scholar] [PubMed]
  14. Loos, M.; Quentmeier, P.; Schuster, T.; Nitsche, U.; Gertler, R.; Keerl, A.; Kocher, T.; Friess, H.; Rosenberg, R. Effect of preoperative radio(chemo)therapy on long-term functional outcome in rectal cancer patients: A systematic review and meta-analysis. Ann. Surg. Oncol. 2013, 20, 1816–1828. [Google Scholar] [PubMed]
  15. Bernier, L.; Balyasnikova, S.; Tait, D.; Brown, G. Watch-and-Wait as a Therapeutic Strategy in Rectal Cancer. Curr. Colorectal Cancer Rep. 2018, 14, 37–55. [Google Scholar] [PubMed]
  16. Dossa, F.; Chesney, T.R.; Acuna, S.A.; Baxter, N.N. A watch-and-wait approach for locally advanced rectal cancer after a clinical complete response following neoadjuvant chemoradiation: A systematic review and meta-analysis. Lancet Gastroenterol. Hepatol. 2017, 2, 501–513. [Google Scholar]
  17. Smith, J.J.; Strombom, P.; Chow, O.S.; Roxburgh, C.S.; 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, e185896. [Google Scholar]
  18. López-Campos, F.; Martín-Martín, M.; Fornell-Pérez, R.; García-Pérez, J.C.; Die-Trill, J.; Fuentes-Mateos, R.; López-Durán, S.; Domínguez-Rullán, J.; Ferreiro, R.; Riquelme-Oliveira, A.; et al. Watch and wait approach in rectal cancer: Current controversies and future directions. World J. Gastroenterol. 2020, 26, 4218–4239. [Google Scholar]
  19. Babic, T.; Lygirou, V.; Rosic, J.; Miladinov, M.; Djikic Rom, A.; Baira, E.; Stroggilos, R.; Pappa, E.; Zoidakis, J.; Krivokapic, Z.; et al. Pilot proteomic study of locally advanced rectal cancer before and after neoadjuvant chemoradiotherapy indicates high metabolic activity in non-responders’ tumor tissue. Proteom. Clin Appl. 2023, 17, 2100116. [Google Scholar]
  20. Szklarczyk, D.; Gable, A.L.; Nastou, K.C.; Lyon, D.; Kirsch, R.; Pyysalo, S.; Doncheva, N.T.; Legeay, M.; Fang, T.; Bork, P.; et al. The STRING database in 2021: Customizable protein-protein networks, and functional characterization of user-uploaded gene/measurement sets. Nucleic Acids Res. 2021, 49, D605–D612. [Google Scholar]
  21. Uhlen, M.; Zhang, C.; Lee, S.; Sjöstedt, E.; Fagerberg, L.; Bidkhori, G.; Benfeitas, R.; Arif, M.; Liu, Z.; Edfors, F.; et al. A pathology atlas of the human cancer transcriptome. Science 2017, 357, eaan2507. [Google Scholar] [CrossRef]
  22. Stelzer, G.; Rosen, N.; Plaschkes, I.; Zimmerman, S.; Twik, M.; Fishilevich, S.; Stein, T.I.; Nudel, R.; Lieder, I.; Mazor, Y.; et al. The GeneCards suite: From gene data mining to disease genome sequence analyses. Curr. Protoc. Bioinform. 2016, 2016, 1.30.1–1.30.33. [Google Scholar]
  23. Safran, M.; Rosen, N.; Twik, M.; BarShir, R.; Stein, T.I.; Dahary, D.; Fishilevich, S.; Lancet, D. The GeneCards Suite. Practical Guide to Life Science Databases; Springer Nature: Singapore, 2021; pp. 27–56. [Google Scholar]
  24. Kim, Y.J.; Sweet, S.M.M.; Egertson, J.D.; Sedgewick, A.J.; Woo, S.; Liao, W.L.; Merrihew, G.E.; Searle, B.C.; Vaske, C.; Heaton, R.; et al. Data-Independent Acquisition Mass Spectrometry to Quantify Protein Levels in FFPE Tumor Biopsies for Molecular Diagnostics. J. Proteome Res. 2019, 18, 426–435. [Google Scholar] [PubMed]
  25. Weke, K.; Kote, S.; Faktor, J.; Al Shboul, S.; Uwugiaren, N.; Brennan, P.M.; Goodlett, D.R.; Hupp, T.R.; Dapic, I. DIA-MS proteome analysis of formalin-fixed paraffin-embedded glioblastoma tissues. Anal. Chim. Acta 2022, 1204, 339695. [Google Scholar] [PubMed]
  26. Amarnani, A.; Capri, J.R.; Souda, P.; Elashoff, D.A.; Lopez, I.A.; Whitelegge, J.P.; Singh, R.R. Quantitative Proteomics Using Formalin-fixed, Paraffin-embedded Biopsy Tissues in Inflammatory Disease. J. Proteom. Bioinform. 2019, 12, 104. [Google Scholar] [CrossRef] [PubMed]
  27. Gao, X.H.; Li, J.; Gong, H.F.; Yu, G.Y.; Liu, P.; Hao, L.Q.; Liu, L.J.; Bai, C.G.; Zhang, W. Comparison of Fresh Frozen Tissue With Formalin-Fixed Paraffin-Embedded Tissue for Mutation Analysis Using a Multi-Gene Panel in Patients with Colorectal Cancer. Front. Oncol. 2020, 10, 310. [Google Scholar]
  28. Repetto, O.; De Re, V.; De Paoli, A.; Belluco, C.; Alessandrini, L.; Canzonieri, V. Identification of Protein Clusters Predictive of Tumor Response in Rectal Cancer Patients Receiving Neoadjuvant Chemo-Radiotherapy. Oncotarget 2017, 8, 28328. Available online: www.impactjournals.com/oncotarget/ (accessed on 15 April 2022).
  29. Chauvin, A.; Wang, C.S.; Geha, S.; Garde-Granger, P.; Mathieu, A.A.; Lacasse, V.; Boisvert, F.M. The response to neoadjuvant chemoradiotherapy with 5-fluorouracil in locally advanced rectal cancer patients: A predictive proteomic signature. Clin. Proteom. 2018, 15, 1–16. [Google Scholar]
  30. Pawlik, T.M.; Keyomarsi, K. Role of cell cycle in mediating sensitivity to radiotherapy. Int. J. Radiat. Oncol. Biol. Phys. 2004, 59, 928–942. [Google Scholar]
  31. Khan, S.U.; Khan, I.M.; Khan, M.U.; Ud Din, M.A.; Khan, M.Z.; Khan, N.M.; Liu, Y. Role of LGMN in tumor development and its progression and connection with the tumor microenvironment. Front. Mol. Biosci. 2023, 10, 1121964. [Google Scholar] [CrossRef]
  32. Wang, Y.; Zhang, S.; Wang, H.; Cui, Y.; Wang, Z.; Cheng, X.; Li, W.; Hou, J.; Ji, Y.; Liu, T. High Level of Legumain Was Correlated With Worse Prognosis and Peritoneal Metastasis in Gastric Cancer Patients. Front. Oncol. 2020, 10, 966. [Google Scholar] [CrossRef]
  33. Yao, Y.; Zhou, D.; Shi, D.; Zhang, H.; Zhan, S.; Shao, X.; Sun, K.; Sun, L.; Wu, G.; Tian, K.; et al. GLI1 overexpression promotes gastric cancer cell proliferation and migration and induces drug resistance by combining with the AKT-mTOR pathway. Biomed. Pharmacother. 2019, 111, 993–1004. [Google Scholar]
  34. Murugan, A.K. mTOR: Role in cancer, metastasis and drug resistance. Semin. Cancer Biol. 2019, 59, 92–111. [Google Scholar] [PubMed]
  35. Janku, F.; McConkey, D.J.; Hong, D.S.; Kurzrock, R. Autophagy as a target for anticancer therapy. Nat. Rev. Clin. Oncol. 2011, 8, 528–539. [Google Scholar] [PubMed]
  36. Saxton, R.A.; Sabatini, D.M. mTOR Signaling in Growth, Metabolism, and Disease. Cell 2017, 168, 960–976. [Google Scholar] [PubMed]
  37. Mossmann, D.; Park, S.; Hall, M.N. mTOR signalling and cellular metabolism are mutual determinants in cancer. Nat. Rev. Cancer 2018, 18, 744–757. [Google Scholar]
  38. Heeran, A.B.; Dunne, M.R.; Morrissey, M.E.; Buckley, C.E.; Clarke, N.; Cannon, A.; Donlon, N.E.; Nugent, T.S.; Durand, M.; Dunne, C.; et al. The protein secretome is altered in rectal cancer tissue compared to normal rectal tissue, and alterations in the secretome induce enhanced innate immune responses. Cancers 2021, 13, 1–18. [Google Scholar]
  39. Lee, S.Y.; Kim, C.H.; Kim, Y.J.; Kwak, H.D.; Ju, J.K.; Kim, H.R. Obesity as an independent predictive factor for pathologic complete response after neoadjuvant chemoradiation in rectal cancer. Ann. Surg. Treat Res. 2019, 96, 116–122. [Google Scholar] [CrossRef]
  40. Clark, W.; Siegel, E.M.; Chen, Y.A.; Zhao, X.; Parsons, C.M.; Hernandez, J.M.; Weber, J.; Thareja, S.; Choi, J.; Shibata, D. Quantitative measures of visceral adiposity and body mass index in predicting rectal cancer outcomes after neoadjuvant chemoradiation. J. Am. Coll. Surg. 2013, 216, 1070–1081. [Google Scholar] [CrossRef]
  41. Huminiecki, L. Magic roundabout is an endothelial-specific ohnolog of ROBO1 which neo-functionalized to an essential new role in angiogenesis. PLoS ONE 2019, 14, e0208952. [Google Scholar]
  42. Jones, C.A.; Nishiya, N.; London, N.R.; Zhu, W.; Sorensen, L.K.; Chan, A.C.; Lim, C.J.; Chen, H.; Zhang, Q.; Schultz, P.G.; et al. Slit2-Robo4 signalling promotes vascular stability by blocking Arf6 activity. Nat. Cell. Biol. 2009, 11, 1325–1331. [Google Scholar]
  43. Wu, M.F.; Liao, C.Y.; Wang, L.Y.; Chang, J.T. The role of Slit-Robo signaling in the regulation of tissue barriers. Tissue Barriers 2017, 5, e1331155. [Google Scholar]
  44. Marinkovic, M.; Stojanovic-Rundic, S.; Stanojevic, A.; Ostojic, M.; Gavrilovic, D.; Jankovic, R.; Maksimovic, N.; Stroggilos, R.; Zoidakis, J.; Castellví-Bel, S.; et al. Exploring novel genetic and hematological predictors of response to neoadjuvant chemoradiotherapy in locally advanced rectal cancer. Front. Genet. 2023, 14, 1245594. [Google Scholar] [PubMed]
  45. Lodge, K.M.; Cowburn, A.S.; Li, W.; Condliffe, A.M. The impact of hypoxia on neutrophil degranulation and consequences for the host. Int. J. Mol. Sci. 2020, 21, 1183. [Google Scholar] [PubMed]
  46. Li, C.F.; He, H.L.; Wang, J.Y.; Huang, H.Y.; Wu, T.F.; Hsing, C.H.; Lee, S.W.; Lee, H.H.; Fang, J.L.; Huang, W.T.; et al. Fibroblast growth factor receptor 2 overexpression is predictive of poor prognosis in rectal cancer patients receiving neoadjuvant chemoradiotherapy. J. Clin. Pathol. 2014, 67, 1056–1061. [Google Scholar] [PubMed]
  47. Gillespie, M.; Jassal, B.; Stephan, R.; Milacic, M.; Rothfels, K.; Senff-Ribeiro, A.; Griss, J.; Sevilla, C.; Matthews, L.; Gong, C.; et al. The reactome pathway knowledgebase 2022. Nucleic Acids Res. 2022, 50, D687–D692. [Google Scholar] [PubMed]
  48. Shi, J.; Zhang, Y.; Yao, B.; Sun, P.; Hao, Y.; Piao, H.; Zhao, X. Role of exosomes in the progression, diagnosis, and treatment of gliomas. Med. Sci. Monit. Int. Med. J. Exp. Clin. Res. 2020, 26, e924023-1–e924023-16. [Google Scholar] [CrossRef]
  49. Nobili, S.; Napoli, C.; Landini, I.; Morganti, M.; Cianchi, F.; Valanzano, R.; Tonelli, F.; Cortesini, C.; Mazzei, T.; Mini, E. Identification of potential pharmacogenomic markers of clinical efficacy of 5-fluorouracil in colorectal cancer. Int. J. Cancer 2011, 128, 1935–1945. [Google Scholar] [CrossRef]
  50. Kim, J.; Lee, J.H.; Iyer, V.R. Global identification of Myc target genes reveals its direct role in mitochondrial biogenesis and its E-box usage in vivo. PLoS ONE 2008, 3, e1798. [Google Scholar]
  51. Aviv, A.; Anderson, J.J.; Shay, J.W. Mutations, Cancer and the Telomere Length Paradox. Trends Cancer 2017, 3, 253–258. [Google Scholar]
  52. Gallicchio, L.; Gadalla, S.M.; Murphy, J.D.; Simonds, N.I. The Effect of Cancer Treatments on Telomere Length: A Systematic Review of the Literature. JNCI J. Natl. Cancer Inst. 2018, 110, 1048–1058. [Google Scholar]
  53. Zou, Z.; Tao, T.; Li, H.; Zhu, X. MTOR signaling pathway and mTOR inhibitors in cancer: Progress and challenges. Cell Biosci. 2020, 10, 1–11. [Google Scholar]
  54. Schieber, M.; Chandel, N.S. ROS function in redox signaling and oxidative stress. Curr. Biol. 2014, 24, R453–R462. [Google Scholar]
  55. Vuletić, A.; Mirjačić Martinović, K.; Tišma Miletić, N.; Zoidakis, J.; Castellvi-Bel, S.; Čavić, M. Cross-Talk Between Tumor Cells Undergoing Epithelial to Mesenchymal Transition and Natural Killer Cells in Tumor Microenvironment in Colorectal Cancer. Front. Cell. Dev. Biol. 2021, 9, 750022. [Google Scholar]
  56. Anufrieva, K.S.; Shender, V.O.; Arapidi, G.P.; Pavlyukov, M.S.; Shakhparonov, M.I.; Shnaider, P.V.; Butenko, I.O.; Lagarkova, M.A.; Govorun, V.M. Therapy-induced stress response is associated with downregulation of pre-mRNA splicing in cancer cells. Genome Med. 2018, 10, 49. [Google Scholar] [PubMed]
  57. Herskind, C.; Fleckenstein, K.; Lohr, J.; Li, C.Y.; Wenz, F.; Lohr, F. Antitumorale wirkung von interferonen und interleukinen in kombination mit strahlentherapie. Teil II: Strahlenbiologische und immunologische strategien. Strahlenther. Onkol. 2004, 180, 331. [Google Scholar]
  58. Gray, M.A.; Stanczak, M.A.; Mantuano, N.R.; Xiao, H.; Pijnenborg, J.F.A.; Malaker, S.A.; Miller, C.L.; Weidenbacher, P.A.; Tanzo, J.T.; Ahn, G.; et al. Targeted glycan degradation potentiates the anticancer immune response in vivo. Nat. Chem. Biol. 2020, 16, 1376–1384. [Google Scholar] [CrossRef]
  59. Dong, X.; Luo, Z.; Wang, Y.; Meng, L.; Duan, Q.; Qiu, L.; Peng, F.; Shen, L. Altered O-glycosylation is associated with inherent radioresistance and malignancy of human laryngeal carcinoma. Exp. Cell. Res. 2018, 362, 302–310. [Google Scholar] [CrossRef]
  60. Park, T. Crk and CrkL as Therapeutic Targets for Cancer Treatment. Cells 2021, 10, 739. [Google Scholar] [CrossRef]
  61. Xie, W.; Qu, L.; Meng, L.; Liu, C.; Wu, J.; Shou, C. PES1 regulates sensitivity of colorectal cancer cells to anticancer drugs. Biochem. Biophys. Res. Commun. 2013, 431, 460–465. [Google Scholar] [CrossRef]
  62. Cristóbal, I.; Manso, R.; Rincón, R.; Caramés, C.; Senin, C.; Borrero, A.; Martínez-Useros, J.; Rodriguez, M.; Zazo, S.; Aguilera, O.; et al. PP2A inhibition is a common event in colorectal cancer and its restoration using FTY720 shows promising therapeutic potential. Mol. Cancer Ther. 2014, 13, 938–947. [Google Scholar] [CrossRef]
  63. Masuda, H.; Baggerly, K.A.; Wang, Y.; Zhang, Y.; Gonzalez-Angulo, A.M.; Meric-Bernstam, F.; Valero, V.; Lehmann, B.D.; Pietenpol, J.A.; Hortobagyi, G.N.; et al. Differential Response to Neoadjuvant Chemotherapy among 7 Triple-Negative Breast Cancer Molecular Subtypes. Clin. Cancer Res. 2013, 19, 5533–5540. [Google Scholar] [CrossRef]
  64. Waitz, R.; Solomon, S.B.; Petre, E.N.; Trumble, A.E.; Fassò, M.; Norton, L.; Allison, J.P. Potent induction of tumor immunity by combining tumor cryoablation with anti-CTLA-4 therapy. Cancer Res. 2012, 72, 430–439. [Google Scholar] [CrossRef] [PubMed]
  65. Cha, M.-K.; Kim, I.-H. Preferential overexpression of glutaredoxin3 in human colon and lung carcinoma. Cancer Epidemiol. 2009, 33, 281–287. [Google Scholar] [CrossRef] [PubMed]
  66. Qian, L.; Lai, X.; Gu, B.; Sun, X. An Immune-Related Gene Signature for Predicting Neoadjuvant Chemoradiotherapy Efficacy in Rectal Carcinoma. Front. Immunol. 2022, 13, 784479. [Google Scholar] [CrossRef]
  67. O’Sullivan, B.; Brierley, J.; Byrd, D.; Bosman, F.; Kehoe, S.; Kossary, C.; Piñeros, M.; Van Eycken, E.; Weir, H.K.; Gospodarowicz, M. The TNM classification of malignant tumours—Towards common understanding and reasonable expectations. Lancet Oncol. 2017, 18, 849–851. [Google Scholar] [CrossRef] [PubMed]
  68. Landberg, T.; Chavaudra, J.; Dobbs, J.; Hanks, G.; Johansson, K.A.; Möller, T. Report 50. J. Int. Comm. Radiat. Units Meas. 2016, 26. [Google Scholar] [CrossRef]
  69. Landberg, T.; Chavaudra, J.; Dobbs, J.; Gerard, J.P.; Hanks, G.; Horiot, J.C. Report 62. J. Int. Comm. Radiat. Units Meas. 2016, 32. [Google Scholar] [CrossRef]
  70. Thies, S.; Langer, R. Tumor regression grading of gastrointestinal carcinomas after neoadjuvant treatment. Front. Oncol. 2013, 3, 262. [Google Scholar] [CrossRef]
  71. Perez-Riverol, Y.; Bai, J.; Bandla, C.; García-Seisdedos, D.; Hewapathirana, S.; Kamatchinathan, S.; Kundu, D.J.; Prakash, A.; Frericks-Zipper, A.; Eisenacher, M.; et al. The PRIDE database resources in 2022: A hub for mass spectrometry-based proteomics evidences. Nucleic Acids Res. 2022, 50, D543–D552. [Google Scholar] [CrossRef]
  72. Deutsch, E.W.; Bandeira, N.; Perez-Riverol, Y.; Sharma, V.; Carver, J.J.; Mendoza, L.; Kundu, D.J.; Wang, S.; Bandla, C.; Kamatchinathan, S.; et al. The ProteomeXchange consortium at 10 years: 2023 update. Nucleic Acids Res. 2023, 51, D1539–D1548. [Google Scholar] [CrossRef]
  73. Fekete, J.T.; Győrffy, B. ROCplot.org: Validating predictive biomarkers of chemotherapy/hormonal therapy/anti-HER2 therapy using transcriptomic data of 3,104 breast cancer patients. Int. J. Cancer 2019, 145, 3140–3151. [Google Scholar] [CrossRef]
  74. El Khouli, R.H.; Macura, K.J.; Barker, P.B.; Habba, M.R.; Jacobs, M.A.; Bluemke, D.A. Relationship of temporal resolution to diagnostic performance for dynamic contrast enhanced MRI of the breast. J. Magn. Reson. Imaging 2009, 30, 999–1004. [Google Scholar] [CrossRef] [PubMed]
  75. Li, Q.; Birkbak, N.J.; Gyorffy, B.; Szallasi, Z.; Eklund, A.C. Jetset: Selecting the Optimal Microarray Probe Set to Represent a Gene. 2011. Available online: http://www.biomedcentral.com/1471-2105/12/474 (accessed on 2 May 2022).
Figure 1. Number of unique proteins identified before and after 50% filtering.
Figure 1. Number of unique proteins identified before and after 50% filtering.
Ijms 24 15412 g001
Figure 2. Principal Component Analysis (a) and a volcano plot (b) of differentially abundant proteins in terms of response to treatment. PCA (a) indicated that responders (blue) compared to non-responders (red) had significantly different proteomics profiles. Proteomics profiles of responders vs. non-responders were compared, and results were represented by volcano plot as −log10 (Welch t-test p-value) plotted against the difference (responders–non-responders). Genes upregulated in responders were colored blue, and those upregulated in non-responders were colored red. Top 10 with highest statistical difference are colored in both groups, and top 5 with highest Welch t-test difference.
Figure 2. Principal Component Analysis (a) and a volcano plot (b) of differentially abundant proteins in terms of response to treatment. PCA (a) indicated that responders (blue) compared to non-responders (red) had significantly different proteomics profiles. Proteomics profiles of responders vs. non-responders were compared, and results were represented by volcano plot as −log10 (Welch t-test p-value) plotted against the difference (responders–non-responders). Genes upregulated in responders were colored blue, and those upregulated in non-responders were colored red. Top 10 with highest statistical difference are colored in both groups, and top 5 with highest Welch t-test difference.
Ijms 24 15412 g002
Figure 3. Hierarchical cluster analysis was performed on the proteomic data. The results indicate the existence of two clusters of patients based on proteomic analysis, which are highly correlated with the clinical response to therapy with two outliers (responders are written in blue, and non-responders are written in pink).
Figure 3. Hierarchical cluster analysis was performed on the proteomic data. The results indicate the existence of two clusters of patients based on proteomic analysis, which are highly correlated with the clinical response to therapy with two outliers (responders are written in blue, and non-responders are written in pink).
Ijms 24 15412 g003
Figure 4. Hierarchical clustering analysis revealed details on the protein abundance of the top 20 candidates among two groups of responders and non-responders. Y axis represents the top 10 overexpressed DEPs in both the responder (pink) and non-responder (blue).
Figure 4. Hierarchical clustering analysis revealed details on the protein abundance of the top 20 candidates among two groups of responders and non-responders. Y axis represents the top 10 overexpressed DEPs in both the responder (pink) and non-responder (blue).
Ijms 24 15412 g004
Figure 5. Enrichment pathway analysis of the differentially expressed proteins between responders and non-responders. Pathway enrichment analysis was conducted using the Metascape program. All statistically enriched terms were first identified by software, and then accumulative hypergeometric p-values and enrichment factors were computed and used for filtering. The remaining relevant terms were subsequently hierarchically clustered into a tree based on the similarity of their gene memberships as measured by Kappa statistics. The obtained results were considered and represented based on biological relevance with respect to rectal cancer biology. KEGG, Reactome, Corum, and Hallmark databases were used for pathways enrichment analysis. (a) Bar graph of enriched signaling pathways of the differentially expressed proteins (915 DEPSs) (colored by p-values). (b) Bar graph of enriched signaling pathways of the differentially expressed proteins overexpressed in non-responder group (700 DEPs) (colored by p-values). (c) Bar graph of enriched signaling pathways of the differentially expressed proteins overexpressed in responder group (215 DEPs) (colored by p-values).
Figure 5. Enrichment pathway analysis of the differentially expressed proteins between responders and non-responders. Pathway enrichment analysis was conducted using the Metascape program. All statistically enriched terms were first identified by software, and then accumulative hypergeometric p-values and enrichment factors were computed and used for filtering. The remaining relevant terms were subsequently hierarchically clustered into a tree based on the similarity of their gene memberships as measured by Kappa statistics. The obtained results were considered and represented based on biological relevance with respect to rectal cancer biology. KEGG, Reactome, Corum, and Hallmark databases were used for pathways enrichment analysis. (a) Bar graph of enriched signaling pathways of the differentially expressed proteins (915 DEPSs) (colored by p-values). (b) Bar graph of enriched signaling pathways of the differentially expressed proteins overexpressed in non-responder group (700 DEPs) (colored by p-values). (c) Bar graph of enriched signaling pathways of the differentially expressed proteins overexpressed in responder group (215 DEPs) (colored by p-values).
Ijms 24 15412 g005
Figure 6. ROC plot of genes overexpressed in non-responder group.
Figure 6. ROC plot of genes overexpressed in non-responder group.
Ijms 24 15412 g006
Figure 7. ROC plots of genes overexpressed in responder group.
Figure 7. ROC plots of genes overexpressed in responder group.
Ijms 24 15412 g007
Table 1. Shortlisted DEPs with characteristics of promising biomarkers enriched in the responder group compared to the non-responder group. Genes with an ROC p-value less than 0.05, an AUC greater than 0.7, and a Mann–Whitney p-value less than 0.05 were characterized as promising biomarkers.
Table 1. Shortlisted DEPs with characteristics of promising biomarkers enriched in the responder group compared to the non-responder group. Genes with an ROC p-value less than 0.05, an AUC greater than 0.7, and a Mann–Whitney p-value less than 0.05 were characterized as promising biomarkers.
GenesGene NameBiological FunctionFOLD ChangeMann–Whitney Test p-ValueAUCROC p-Value
CRKLCRK-like Proto-OncogeneActivate the RAS and JUN kinase signaling pathways and transform fibroblasts in a RAS-dependent fashion1.40.00630.7488.00 × 10−4
LAP3Leucine Aminopeptidase 3Predicted to enable peptidase activity; involved in proteolysis.1.50.00590.7469.30 × 10−4
THTPAThiamine TriphosphataseEncodes an enzyme that catalyzes the biosynthesis of thiamine disphophate (vitamin B1) by hydrolysis of thiamine triphosphate1.40.00890.7381.20 × 10−3
PES1Pescadillo Ribosomal Biogenesis Factor 1Encodes a nuclear protein that contains a breast-cancer-associated gene 1 (BRCA1) C-terminal interaction domain1.20.00960.7361.40 × 10−3
PPP2R5EProtein Phosphatase 2 Regulatory Subunit B’EpsilonProtein phosphatase 2A is implicated in the negative control of cell growth and division1.20.010.731.60 × 10−3
IFI30IFI30 Lysosomal Thiol ReductaseThis enzyme has an important role in MHC class II-restricted antigen processing1.50.0110.7282.30 × 10−3
C17orf75Chromosome 17 Open Reading Frame 75Involved in intracellular protein transport and vesicle tethering to Golgi1.30.0150.7223.00 × 10−3
QDPRQuinoid Dihydropteridine ReductaseThis gene encodes the enzyme dihydropteridine reductase, which catalyzes the NADH-mediated reduction of quinonoid dihydrobiopterin1.30.0130.7194.10 × 10−3
RRM2BRibonucleotide Reductase Regulatory TP53 Inducible Subunit M2BThis heterotetrameric enzyme catalyzes the conversion of ribonucleoside diphosphates to deoxyribonucleoside diphosphates1.30.0160.7193.60 × 10−3
GLRXGlutaredoxinThis enzyme highly contributes to the antioxidant defense system1.70.020.7094.70 × 10−3
USO1USO1 Vesicle Transport FactorPeripheral membrane protein that recycles between the cytosol and the Golgi apparatus during interphase1.20.0220.7076.60 × 10−3
ARAFA-Raf Proto-Oncogene, Serine/Threonine KinaseInvolved in negative regulation of apoptotic process, regulation of TOR signaling, and regulation of cellular protein metabolic process1.20.0240.7066.80 × 10−3
CTBSChitobiaseLysosomal glycosidase is involved in degradation of asparagine-linked oligosaccharides on glycoproteins1.30.0240.7056.70 × 10−3
SNRPD3Small Nuclear Ribonucleoprotein D3 PolypeptideThis gene encodes a core component of the spliceosome, which is a nuclear ribonucleoprotein complex that functions in pre-mRNA splicing1.20.0270.78.30 × 10−3
Table 2. Shortlisted DEPs with characteristics of promising biomarkers enriched in the non-responder group compared to the responder group. Genes with a ROC p-value less than 0.05, an AUC greater than 0.7, and a Mann–Whitney p-value less than 0.05 were characterized as promising biomarkers.
Table 2. Shortlisted DEPs with characteristics of promising biomarkers enriched in the non-responder group compared to the responder group. Genes with a ROC p-value less than 0.05, an AUC greater than 0.7, and a Mann–Whitney p-value less than 0.05 were characterized as promising biomarkers.
GenesGene NameBiological FunctionFOLD ChangeMann–Whitney Test p-ValueAUCROC p-Value
COPB1COPI Coat Complex Subunit Beta 1This gene encodes a protein subunit of the coatomer complex to mediate biosynthetic protein transport from the ER via the Golgi up to the trans-Golgi network1.10.00610.7494.00 × 10−4
MGLLMonoglyceride LipaseCatalyzes the conversion of monoacylglycerides to free fatty acids and glycerol, and gene expression may play a role in cancer tumorigenesis and metastasis1.30.00830.741.10 × 10−3
HAS1Hyaluronan Synthase 1Essential to hyaluronan synthesis, it has a structural role in tissue architectures and regulates cell adhesion, migration, and differentiation.1.90.0120.7291.80 × 10−3
TALDO1Transaldolase 1Transaldolase 1 is a key enzyme of the nonoxidative pentose phosphate pathway, providing ribose-5-phosphate for nucleic acid synthesis and NADPH for lipid biosynthesis.1.20.0120.7292.30 × 10−3
DNAH9Dynein Axonemal Heavy Chain 9This gene encodes the heavy chain subunit of axonemal dynein, a large multi-subunit molecular motor.1.40.0130.7242.90 × 10−3
KDELR3KDEL Endoplasmic Reticulum Protein Retention Receptor 3This gene encodes a member of the KDEL endoplasmic reticulum protein retention receptor family1.60.0150.7224.10 × 10−3
HLA-DPB1Major Histocompatibility Complex, Class II, DP Beta 1It plays a central role in the immune system by presenting peptides derived from extracellular proteins1.80.020.7095.20 × 10−3
RBP3Retinol-binding Protein 3Interphotoreceptor retinol-binding protein is found primarily in the interphotoreceptor matrix of the retina between the retinal pigment epithelium and the photoreceptor cells *1.80.0220.7085.30 × 10−3
STAP2Signal-transducing Adaptor Family Member 2This gene encodes the substrate of breast tumor kinase, an Src-type non-receptor tyrosine kinase1.20.0250.7036.90 × 10−3
* According to the Human Protein Atlas, several cases of colorectal cancers displayed moderate cytoplasmic and membranous staining [21].
Table 3. Drug targets among top proteins that are differentially overexpressed in the group of non-responders.
Table 3. Drug targets among top proteins that are differentially overexpressed in the group of non-responders.
GeneProtein Name Encoded by GeneDrugDRUGBANK IDDrug TypeUsageDrug Approval StatusWelch
t-Test
p Value
QPRTQuinolinate PhosphoribosyltransferaseNiacinDB00627B vitaminHyperlipidemia, dyslipidemia, hypertriglyceridemiaApproved, investigational, nutraceutical4.43 × 10−5
CLCA4Chloride Channel Accessory 4TalniflumateDB09295Small molecule, CaCC blocker, and Cl−-/HCO3exchange inhibitorCystic fibrosis, chronic obstructive pulmonary disease (COPD), and asthmaExperimental1.11 × 10−2
ATG4BAutophagy-related 4B Cysteine PeptidaseEsomeprazoleDB00736Proton pump inhibitorGastroesophageal reflux disease (GERD) for gastric protection to prevent recurrence of stomach ulcers or gastric damageApproved, investigational1.89 × 10−4
ATG4B NimodipineDB00393Calcium channel blockerActs primarily on vascular smooth muscle cells; improves the neurologic outcome following subarachnoid hemorrhage from ruptured intracranial aneurysm.Approved, investigational
PTGS2Prostaglandin Endoperoxide Synthase 2Diclofenac, Acetylsalicylic acid, Ibuprofen, Rofecoxib, AcetaminophenDB00586, DB00945, DB01050, DB00533, DB00316COX inhibitors, anti-inflammatory agentsTherapy for acute and chronic pain and inflammation from a variety of causesApproved2.08 × 10−4
Table 4. Clinical data of the study cohort.
Table 4. Clinical data of the study cohort.
CharacteristicsResponders
N (%)
Non-Responders
N (%)
p-Value *
GenderMale3 (33.3)7 (63.6)
4 (36.4)
0.4
Female6 (66.7)
Age (years)Mean (SD)64.0 (6.7)62.4 (10.3)1.0
Median (Range)66.0 (50–72)64.0 (48–83)
UICC stagingII0.0 (0)2.0 (18.2)0.5
III9.0 (100)9.0 (81.8)
Tumor grade17.0 (77.8)8.0 (72.7)0.8
22.0 (22.2)3.0 (27.3)
30.0 (0.0)0.0 (0.0)
Tumor localizationInferior rectum (<5 cm)5.0 (55.6)9.0 (81.8)0.4
Mid rectum (5–10 cm)4.0 (44.4)2.0 (18.2)
Superior rectum (>10 cm)0.0 (0)0.0 (0)
Acute toxicitywithout1.0 (12.5)2.0 (18.2)0.7
with8.0 (87.5)9.0 (81.8)
Tumor Regression Grade (TRG)18.0 (88.89)/
21.0 (11.11)/
3/2.0 (18.18)
4/7.0 (63.64)
5/2.0 (18.18)
Total 9 (100)11 (100)
* Chi-squared test with Yates correction (two-tailed).
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Stanojevic, A.; Samiotaki, M.; Lygirou, V.; Marinkovic, M.; Nikolic, V.; Stojanovic-Rundic, S.; Jankovic, R.; Vlahou, A.; Panayotou, G.; Fijneman, R.J.A.; et al. Data-Independent Acquisition Mass Spectrometry Analysis of FFPE Rectal Cancer Samples Offers In-Depth Proteomics Characterization of the Response to Neoadjuvant Chemoradiotherapy. Int. J. Mol. Sci. 2023, 24, 15412. https://doi.org/10.3390/ijms242015412

AMA Style

Stanojevic A, Samiotaki M, Lygirou V, Marinkovic M, Nikolic V, Stojanovic-Rundic S, Jankovic R, Vlahou A, Panayotou G, Fijneman RJA, et al. Data-Independent Acquisition Mass Spectrometry Analysis of FFPE Rectal Cancer Samples Offers In-Depth Proteomics Characterization of the Response to Neoadjuvant Chemoradiotherapy. International Journal of Molecular Sciences. 2023; 24(20):15412. https://doi.org/10.3390/ijms242015412

Chicago/Turabian Style

Stanojevic, Aleksandra, Martina Samiotaki, Vasiliki Lygirou, Mladen Marinkovic, Vladimir Nikolic, Suzana Stojanovic-Rundic, Radmila Jankovic, Antonia Vlahou, George Panayotou, Remond J. A. Fijneman, and et al. 2023. "Data-Independent Acquisition Mass Spectrometry Analysis of FFPE Rectal Cancer Samples Offers In-Depth Proteomics Characterization of the Response to Neoadjuvant Chemoradiotherapy" International Journal of Molecular Sciences 24, no. 20: 15412. https://doi.org/10.3390/ijms242015412

APA Style

Stanojevic, A., Samiotaki, M., Lygirou, V., Marinkovic, M., Nikolic, V., Stojanovic-Rundic, S., Jankovic, R., Vlahou, A., Panayotou, G., Fijneman, R. J. A., Castellví-Bel, S., Zoidakis, J., & Cavic, M. (2023). Data-Independent Acquisition Mass Spectrometry Analysis of FFPE Rectal Cancer Samples Offers In-Depth Proteomics Characterization of the Response to Neoadjuvant Chemoradiotherapy. International Journal of Molecular Sciences, 24(20), 15412. https://doi.org/10.3390/ijms242015412

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop