- Open Access
The response to neoadjuvant chemoradiotherapy with 5-fluorouracil in locally advanced rectal cancer patients: a predictive proteomic signature
Clinical Proteomicsvolume 15, Article number: 16 (2018)
Colorectal cancer is the third most common and the fourth most lethal cancer in the world. In the majority of cases, patients are diagnosed at an advanced stage or even metastatic, thus explaining the high mortality. The standard treatment for patients with locally advanced non-metastatic rectal cancer is neoadjuvant radio-chemotherapy (NRCT) with 5-fluorouracil (5-FU) followed by surgery, but the resistance rate to this treatment remains high with approximately 30% of non-responders. The lack of evidence available in clinical practice to predict NRCT resistance to 5-FU and to guide clinical practice therefore encourages the search for biomarkers of this resistance.
From twenty-three formalin-fixed paraffin-embedded (FFPE) biopsies performed before NRCT with 5-FU of locally advanced non-metastatic rectal cancer patients, we extracted and analysed the tumor proteome of these patients. From clinical data, we were able to classify the twenty-three patients in our cohort into three treatment response groups: non-responders (NR), partial responders (PR) and total responders (TR), and to compare the proteomes of these different groups.
We have highlighted 384 differentially abundant proteins between NR and PR, 248 between NR and TR and 417 between PR and TR. Among these proteins, we have identified many differentially abundant proteins identified as having a role in cancer (IFIT1, FASTKD2, PIP4K2B, ARID1B, SLC25A33: overexpressed in TR; CALD1, CPA3, B3GALT5, CD177, RIPK1: overexpressed in NR). We have also identified that DPYD, the main degradation enzyme of 5-FU, was overexpressed in NR, as well as several ribosomal and mitochondrial proteins also overexpressed in NR. Data are available via ProteomeXchange with identifier PXD008440.
From these retrospective study, we implemented a protein extraction protocol from FFPE biopsy to highlight protein differences between different response groups to RCTN with 5-FU in patients with locally advanced non-metastatic rectal cancer. These results will pave the way for a larger cohort for better sensitivity and specificity of the signature to guide decisions in the choice of treatment.
Colorectal cancer (CRC) is the third most common cancer (9.7%) in the world behind lung (13%) and breast (12%) cancers and the fourth most lethal (8.5%) behind lung (19%), liver (9.1%) and stomach (8.8%) cancers . This high mortality can be explained by a late detection of the cancer, so patients are often diagnosed with locally advanced tumor completely invading the wall of the rectum (T3) or peripheral tissues or organs (T4), and often regional lymph node metastases are found (N1 or N2). In fact, the 5-year survival is estimated at 93, 78, 66 and 8% for stages I, II, III and IV respectively . In North America, the standard treatment of locally advanced rectum cancer is neoadjuvant radiotherapy with concomitant chemotherapy with the intravenous 5-fluorouracil (5-FU) or its oral analogue Capecitabine (Xeloda®) followed by total mesorectal excision (TME) surgery . This neoadjuvant treatment has been shown to reduce tumor infiltration and decrease the tumor charge (down-staging effect), which is important to increase the complete resection rate (R0, healthy circumferential margin) of surgery and thus improve loco-regional tumor control, disease-free survival (DFS) as well as quality of life by preserving sphincter, urinary and sexual functions [4,5,6,7,8,9,10].
While 70% of patients show an objective response (OR), including up to 15% with a complete remission [10,11,12], a subset of patients with the same tumor stage and having been treated with the same technique do not respond favourably, suggesting a resistant phenomenon. Moreover, considering the severe secondary effects observed such as radiodermatitis and proctitis, there is a strong need in clinic to have reliable biomarkers in order to predict the outcome of the treatment and to identify the sub-group of patients with resistant phenotype, permitting to avoid long, costly and ineffective procedures but with an adequate personalized treatment.
Many tests exist to detect CRC such as faecal occult blood test (FOBT), double contrast barium enema (DCBE) or colonoscopy, but biopsy remains one of the best ways to get a definitive diagnosis. A new blood test, the ColonSentry® test, was recently developed and is based on the expression profile of seven genes ANXA3, CLEC4D, LMNB1, PRRG4, TNFAIP6, VNN1 and IL2RB . This rapid and non-invasive test represents the world’s first commercially available test for stratifying current risk for CRC. Of all these tests, it appears that the majority of them consist of screening tests but the clinical markers currently used for rectal cancer are not able to predict the individual response of the patient that would allow personalization of the treatment [14, 15]. Relevant biological factors to direct the patient to personalized treatment are lacking to guide the clinician’s decisions. The vast majority of predictive and prognostic biomarkers identified in rectal cancer generally involve genetic mutations such as the activating mutation in KRAS (G12D) which is a predictive biomarker of resistance to EGFR-directed therapy with the monoclonal antibodies cetuximab or panitumumab . Furthermore, response to radiotherapy implicates a lot of factors in cell cycle arrest, DNA damage repair and apoptosis. So far, no reliable biomarker in clinic to predict the response pattern to neoadjuvant radio-chemotherapy, which is well associated with disease-free survival.
In this study, we focused on profiling the tumor protein signature of different groups of responders to NRCT with 5-FU: non-responders, partial and total responders (NR, PR and TR respectively) in patients with locally advanced rectal cancer. The purposes of this study are (1) to extract the protein content from rectal formalin-fixed paraffin-embedded (FFPE) biopsies and to study it by mass spectrometry, (2) to identify the differentially abundant proteins (3) to identify, from these proteins, the signaling pathways predominantly represented in each group and finally, (4) to discuss potential predictive biomarkers of the response to the NRCT with 5-FU.
Study design and clinical protocol
This is a retrospective study. The eligibility includes all patients diagnosed between September 2013 and September 2015 in the Centre Hospitalier Universitaire de Sherbrooke (Québec, Canada) with a pathologically confirmed locally advanced rectal adenocarcinoma (T3–T4, any N, M0/any T, N1–N2, M0), treated by standard neoadjuvant radiotherapy and chemotherapy followed by the surgery. All participants gave informed consent before enrolment in this study. The radiotherapy is composed by 45 Gy in 25 fractions, 5 days per week for 5 weeks, ± 5.4 Gy boost by 3D technique. The concomitant chemotherapy was performed by 1600 mg/m2 intravenous 5-FU or oral capecitabine (Fig. 1a). All secondary effects were reviewed and noted according to the Common Terminology Criteria for Adverse Events (CTCAE) version 4.0. The TNM stage before the treatment (iTNM) was used according to the routine clinical examination including magnetic resonance imaging (MRI) for all cases (Table 1) and the post-treatment TNM stage (ypTNM) was defined according to the report of pathology (Table 3). The cancer stage was finally defined according to the American Joint Committee on Cancer (AJCC). The tumor tissues of initial biopsy were used for the study. Patients were then classified into three distinct groups depending on their pathologic response to treatment evaluated with the AJCC tumor regression grading (TRG) system . TRG0 corresponds to a complete regression (non-viable cancer cells), TRG1 corresponds to a near complete regression (single or small groups of tumor cells), TRG2 corresponds to a moderate regression (residual cancer outgrown by fibrosis) and TRG3 corresponds to a minimal/absent regression (minimal or no tumor cells killed). For the sake of clarity, we have associated TRG0 to the total responders (TR), TRG1 and TRG2 to the partial responders (PR) and TRG3 to the non-responders (NR). The study was approved by ethic committee of Cancer Research Center, CHUS (CRC-CHUS). Patients monitoring was performed on a regular basis to index recurrences.
Biopsy sample collection
Biopsies were collected, rapidly fixed in formalin and sent to the pathology laboratory of the Centre Hospitalier Universitaire de Sherbrooke (Québec, Canada), according to the clinical pathological process. They have undergone washings (alcohols and solvents) in order to be degreased and dehydrated before their inclusion in paraffin. Cuts of approximately 3 μm in thickness were made on the paraffin blocks and analysed after staining with haematoxylin and eosin (H&E stain). The areas of interest were selected by the pathologist (Fig. 2a) and punches were made on the paraffin blocks (Fig. 2b) using a Tissue Microarray (TMA) Master II (3DHISTECH) (Fig. 1b).
Dewaxing and homogenisation of the FFPE biopsies
100% xylene was added (100 µl for 1 mg of biopsy) for 1 h under agitation. Xylene is replaced by 50% ethanol to dehydrate the biopsy for 1 h under agitation. After dehydration, FFPE biopsies were homogenised in MS grade water (300 µl for 1 mg of remaining biopsy) using a TissueRuptor (Qiagen). Biopsies were then incubated in 2.5% SDS, 50 mM Tris pH 7.5, 10 mM dithiothréitol at 98 °C for 1 h. After centrifugation at 13,000 rpm for 1 min, the supernatant was collected (Fig. 1c).
Concentration, reduction and alkylation
Solubilised biopsies were concentrated (250–50 µl) using Nanosep 3 K Omega (Pall Life Sciences) and treated with 1X Laemmli and 2 mM dithiothréitol at 98 °C for 5 min (reduction). 50 mM iodoacetamide was added for 30 min and incubated in the dark for alkylation (Fig. 1c).
In-gel separation and tryptic digestion
50 µg of concentrated biopsies were separated on NuPAGE gel 4–12% bis–Tris 1.5 mm × 10 well (Invitrogen). After migration, the gel was stained with Coomassie Blue (SimplyBlue SafeStain) for 1 h under agitation and then discoloured overnight under agitation. Each lane was divided into four gel slices. Each gel slice was cut into cubes (ca. 1 × 1 mm) and transferred in microcentrifuge tube. Then, they were washed with MS grade water for 15 min under agitation then with 50% acetonitrile for 15 min under agitation until gel slice shrunk and became opaque. The gel slices were then subjected to a wash cycle (20 mM ammonium bicarbonate for 15 min under agitation and 20 mM ammonium bicarbonate/100% acetonitrile (1:1, vol/vol) for 15 min under agitation) until the gel became clear. The samples were then dehydrated by washing with 100% acetonitrile for 5 min under agitation until the precipitate became white and solid. Acetonitrile was removed and dehydrated samples were digested in 12.5 ng/ml Trypsin Gold (Promega) modified in 50 mM acetic acid and 20 mM ammonium bicarbonate, at 37 °C overnight under agitation (Fig. 1c).
After the tryptic digestion, the supernatant was collected in a new microcentrifuge tube. Digested samples were next incubated with 100% acetonitrile (volume equal to the digestion volume in the previous step) for 30 min at 37 °C under agitation. The supernatant was collected and digested samples underwent two incubations in 1% formic acid for 20 min at 37 °C under agitation for extracting peptides from gel. The samples were then dehydrated by washing with 100% acetonitrile for 10 min at 37 °C under agitation until the precipitate became white and solid.
Peptides were resuspended in 30 µl of 0.1% trifluoroacetic acid. With a C18 column Zip-Tip (Thermo Scientific) previously washed with 100% acetonitrile and 0.1% trifluoroacetic acid, peptides were loaded into the column by doing 10 up and down, washed with 0.1% trifluoroacetic acid and eluted with 30 µl of 50% acetonitrile/1% formic acid by doing up and down in a new microcentrifuge tube. Desalted peptides were then dried in a speedvac and resuspended in 1% formic acid (Fig. 1c).
High performance liquid chromatography separation, coupled to mass spectrometry (HPLC–MS/MS)
Trypsin-digested purified peptides were separated using a Dionex Ultimate 3000 nanoHPLC system. 2 µg of the sample in 1% formic acid (v:v) were loaded with a constant flow of 4 µL/min on an Acclaim PepMap100 C18 column (0.3 mm id × 5 mm). After trap enrichment, peptides were eluted in a PepMap C18 nanocolumn (75 µm × 50 cm) with a linear gradient of 5–35% solvent B over 240 min with a constant flow of 200 nL/min.
Mass spectrometry analysis
The HPLC system was coupled to an OrbiTrap QExactive mass spectrometer (Thermo Fisher Scientific Inc) via an EasySpray source. The spray voltage was set to 2.0 kV and the temperature of the column was set to 40 °C. Full scan MS survey spectra (m/z 350–1600) in profile mode were acquired in the Orbitrap with a resolution of 70,000 after accumulation of 1,000,000 ions. The ten most intense peptide ions from the preview scan in the Orbitrap were fragmented by collision induced dissociation (normalized collision energy 35% and resolution of 17,500) after the accumulation of 50,000 ions. Maximal filling times were 250 ms for the full scans and 60 ms for the MS/MS scans. Precursor ion charge state screening was enabled and all unassigned charge states as well as singly, 7 and 8 charged species were rejected. The dynamic exclusion list was restricted to a maximum of 500 entries with a maximum retention period of 40 s and a relative mass window of 10 ppm. The lock mass option was enabled for survey scans to improve mass accuracy.
Quantification and bioinformatics analysis
Data were processed, searched and quantified using the MaxQuant software package version 18.104.22.168 with the protein database from UniProtKB (Homo sapiens, 16/07/2013, 88,354 entries) . The following settings were used for the MaxQuant analysis: 2 miscleavages allowed; fixed modification was carbamidomethylation on cysteine; enzyme was trypsin (K/R not before P); variable modifications included in the analysis were methionine oxidation and protein N-terminal acetylation; mass tolerance of 7 ppm for precursor ions and 20 ppm for fragment ions. We performed quantification based on protein intensities. First, we calculated the sum of the protein intensities obtained for each biopsy. We then divided these protein intensities sums by the lowest protein intensity sum of the cohort to obtain a ratio. Each protein intensity for each biopsy was then normalized to this ratio so that the sums of the protein intensities for each biopsy were equal. Finally, the analysis of mass spectrometry data were next carried out with the Perseus software package  and subjected to a normalization (z-score) and an ANOVA statistical test (comparison of the three groups) or Student’s T test (comparison of two groups), both corrected with a p value of 0.05. The mass spectrometry proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE [20, 21] partner repository with the dataset identifier PXD008440.
The different groups of patients were first compared according to clinicopathological parameters. Due to a small number of patients per group, Fisher’s exact test was done with the R software . When the p value was less than 0.05, the result was considered to have statistical significance.
Baseline population characteristics
As the first step of study for the purpose of technique, 23 participants were enrolled and studied. Of the 23 patients, 19 were men (83%) and 4 were women (17%). The median age was 58.6 ± 11.1 years (31–75). Ten patients (43%) had low rectal cancer (within 5 cm from the anal verge), 11 patients (48%) had middle rectal cancer (between 5 and 10 cm from the anal verge) and 2 patients (9%) had proximal cancer (over of 10 cm from the anal verge) (Table 1). Tumor differentiation degrees defining by imaging were 5 well differentiated (22%), 9 moderately differentiated (39%) and 9 not specified (39%). TNM stages before NRCT (iTNM) were 3 stage IIA (13%), 1 stage IIIA (4%), 2 stage IIIA-IIIB (9%), 4 stage IIIB (17%), 12 stage IIIB-IIIC (52%) and 1 stage IIIC (4%) (Table 1).
Surgical treatment and complications
Neoadjuvant radiotherapy was delivered to the pelvis with a dose of 45 Gy followed or not by a boost of 5.4 Gy for 21 patients (91%) and 1 patient (4%) respectively or with a dose of 50.5 (4%). Neoadjuvant chemotherapy was concurrently performed on all patients with Xeloda® (oral) or 5-fluorouracil (intravenous) treatment (96 and 4% respectively). After NRCT, 16 patients (70%) had suffered from radiodermatitis (severity degrees are shown in Table 2). 12% of patients had suffered from proctitis and no cases of colitis or mucositis had been listed. On average, surgery delay after the beginning of treatment was 8.6 ± 2.8 weeks (1.4–13.1). Among the 23 patients, 3 patients refused the surgery. In these cases, the response is measured by MRI imaging after the treatment. The others 20 patients had a low anterior surgery (85%), or abdomino-perineal surgery (8%) even an endo-anal surgery (4%) (Table 2).
Neoadjuvant chemoradiotherapy post-treatment results
After NRCT and surgery, the response to NRCT has been reported by pathological analysis for all patients operated and noted as standard pathological TNM classification (ypTNM) (Table 3). According to the definition of response and down-staging effects, all patients operated have been attributed in three distinct groups: NR, PR and TR. Among the 23 patients, there were 7 NR (30%), 10 PR (43%) and 6 TR (26%). Perineural and lymphovascular invasion (PNI and LVI) are searched on the surgical piece: 2 and 8 cases of PNI and LVI are listed respectively.
Identification of differentially abundant proteins between groups
We were able to identify 5329 proteins among the 23 patients of our study with on average more than 4000 proteins per sample. Using Perseus software package , we had first checked distribution of protein intensities for each duplicate biopsy and Pearson correlation coefficient r for each duplicate biopsy to ensure low variability within the same biopsy (n = 1 and n = 2) (Fig. 3).
From that, we were able to identify 402 differentially abundant proteins between the three groups of responders (Additional file 1: Table S1) (Fig. 4a), 384 between NR and PR (180 over-represented in NR, 204 in PR) (Additional file 2: Table S2) (Fig. 4b), 248 between NR and TR (139 over-represented in NR, 109 in TR) (Additional file 3: Table S3) (Fig. 4c) and 417 between PR and TR (294 over-represented in PR and 123 in TR) (Additional file 4: Table S4) (Fig. 4d). Thanks to the principal component analysis (PCA) realized with the Perseus software package  (Fig. 4e), we can observe the groupings of the three responder groups NR, PR and TR. We also performed our analyses by grouping our patients into only two groups: responders and NR, (1) NR versus PR + TR, or (2) NR + PR versus TR. In the first case (1), we found 282 differentially abundant proteins, 139 of which were over-represented in NR and 143 in responders (PR + TR) (Additional file 5: Table S5). In the second case (2), 368 proteins were differentially abundant, of which 93 were over-represented in TR and 275 in non-responders (NR + PR) (Additional file 6: Table S6) (Table 4).
In terms of protein identification (Additional file 7: Table S7), among the proteins expressed in more than 75% of the TR (≥ 5 patients) and in less than 15% of the NR (≤ 1 patient), we identified the following proteins: PVR, IFIT1, F12, FASTKD2, BPGM, PIP4K2B, MCMBP, CHD1L, ENSA, PVRL1, TRMT5, CNOT10, SLC25A33, FTO, AMDHD2, PHPT1, SLC5A6, LSM12, TSPAN6, CBX1, NOSIP, TSSC4, ARID1B, ALDH3A1, AAMDC, GTF2E1, SNX7 and STX10. In contrast, among the proteins present in more than 75% of the NR (≥ 5 patients) and in less than 15% of the TR (≤ 1 patient), we identified the following proteins (numbers in parentheses representing the relative abundance increase factor of NR over TR for each protein): CALD1 (4.4), CPA3 (9.4), B3GALT5 (55.9), CNNM4 (51.3), MTIF3 (4.1) and CD177 (2). By grouping our patients in two groups, for the first case (1), among the proteins expressed in more than 75% of the responders (PR + TR; ≥ 14 patients) and in less than 15% of the NR (≤ 1 patient), we identified FASTKD2, SLC25A33, FTO and AMDHD2 (proteins already mentioned above). In contrast, among the proteins present in less than 15% of the responders (PR + TR; ≤ 2 patients) and in more than 75% of the NR (≥ 5 patients), we only identified RIPK1 (3.7). For the second case (2), among the proteins expressed in more than 75% of the non-responders (NR + PR; ≥ 14 patients) and in less than 15% of the TR (≤ 1 patient), we only identified CALD1 (already mentioned above). In contrast, among the proteins present in less than 15% of the non-responders (NR + PR; ≤ 2 patients) and in more than 75% of the TR (≥ 5 patients), we identified IFIT1, ABLIM1, CNOT10 (proteins already mentioned above) and ERBB2.
From Human Protein Atlas available from http://www.proteinatlas.org, we analysed the expression of the seven most abundant proteins in non-responders patients: CALD1, CPA3, B3GALT5, CNNM4, MTIF3, CD177 and RIPK1. These seven candidates are expressed in the normal rectal tissue at the protein level (Fig. 5a) and at the gene level (Fig. 5b). For CPA3 and CALD1, protein expression was not detected in the rectum but is expressed in the small bowel and colon respectively. In addition, there is a strong expression of their mRNAs in the rectal tissue. From Colorectal Cancer Atlas available from http://colonatlas.org/index.html , we were able to verify the expression of these seven proteins in 22 colorectal cancer cell lines as well as the presence of gene sequence variants (Fig. 5c). Many mutations have been identified for these genes but none are listed for involvement in deregulated signaling pathways in colorectal cancer.
We also observed the proteins involved in the metabolism of 5-FU. In the differentially abundant proteins, we find in particular the protein DPYD (dihydropyrimidine dehydrogenase), the enzyme involved in the hepatic catabolism of 5-FU. This protein is found to be more abundant in NR compared to TR (P = 3.18 × 10−2), in PR compared to TR (P = 2.63 × 10−2) and in NR + PR compared to TR (P = 1.87 × 10−2). We also find the protein TYMP (thymidine phosphorylase), the enzyme involved in the transformation of 5-FU into FdUMP (fluorodeoxyuridine monophosphate), which is the TS inhibitor (thymidylate synthase), the main target of 5-FU. This enzyme is found to be more abundant in PR compared to NR (P = 3.98 × 10−5), in PR compared to TR (P = 6.81 × 10−3) and in PR + TR compared to NR (P = 2.99 × 10−3).
From these differentially abundant proteins identified above, we looked for the biological processes predominantly represented in each group of responders for each comparison (Additional files 1–6: Tables S1–S6). These analyses were carried out using the DAVID (Database for Annotation, Visualization and Integrated Discovery) v6.8 software (https://david.ncifcrf.gov/) [24, 25]. Among these biological processes, there is a strong representation of ribosomal proteins (RP; RPS and RPL) and mitochondrial proteins (NDUF proteins, included in complexe I of the respiratory chain) in NR compared to the other groups as we can see in Fig. 6 obtained with the STRING software (https://string-db.org/) .
Neoadjuvant radio-chemotherapy with 5-FU is the standard treatment protocol for patients with locally advanced non-metastatic rectal cancer. However, there is still a very high resistance rate to this treatment (30% on average) and the mechanism for that individual resistance is not clearly understood. Studies looking for biomarkers predictive of this resistance have not yet led to the use of these biomarkers in the clinical decisions. FFPE biopsies have long been the subject of proteomic studies since they allow inexpensive and indefinitely storage at room temperature. Indeed, it has been shown that their protein composition after 3 or 15 years of storage was unchanged  and that a 92% overlap in the protein composition between a FFPE and fresh-frozen biopsies was obtained . This model has even been described as a “treasure” for retrospective studies . From the literature, we have developed a protein extraction protocol from FFPE biopsies of patients with locally advanced non-metastatic rectal cancer (Fig. 1) which allowed us to extract a satisfactory quantity of proteins identified by mass spectrometry, with an average of 4000 proteins identified by samples (Additional file 7: Table S7). In addition, the results presented in Fig. 3 with the distributions of protein intensities and multi-scatter plots for each biopsy allowed us to ensure low variability within duplicate samples. Protein compositions were finally compared between different treatment response groups: NR, PR and TR.
Firstly, we were interested in the proteins mainly abundant in each of these groups. For TR, we identified thirty proteins expressed in more than 75% of individuals in this response group and in less than 15% of NR: PVR, IFIT1, F12, FASTKD2, BPGM, PIP4K2B, MCMBP, CHD1L, ENSA, PVRL1, TRMT5, CNOT10, SLC25A33, FTO, AMDHD2, PHPT1, SLC5A6, LSM12, TSPAN6, CBX1, NOSIP, TSSC4, ARID1B, ALDH3A1, AAMDC, GTF2E1, SNX7, STX10, ABLIM1 and ERBB2.
Many of them have been identified in the past for cancer. In particular, the case of IFIT1 (interferon-induced protein with tetratricopeptide repeats 1) can be discussed. A study reported that high levels of IFIT1 mRNA were associated with radiation resistance in breast cancer . In parallel, another team studied the protein expression of IFIT1 and associated a high level of this protein with better local survival without relapse . They also concluded that protein expression may be a better measure for IFIT1 as it has been shown that this protein inhibits translation and cell growth in a dose-dependent manner and can regulate its own translation [31, 32].
Concerning FASTKD2 (FAST kinase domain-containing protein 2, mitochondrial), it was shown to be regulated by NRIF3 (nuclear receptor-interacting factor 3) and the DIF-1 (interferon regulatory factor 2-binding protein 2) complex and this regulation modulated apoptosis in mammary and pancreatic cancer cells [33, 34]. Indeed, the absence of FASTKD2 repression seems to lead to apoptosis in these cells . A lower proportion of FASTKD2 in non-responder individuals in our study may be responsible for resistance to treatment-induced apoptosis in these individuals.
We also identified PIP4K2B (phosphatidylinositol 5-phosphate 4-kinase type-2 beta), of which a low expression has been associated with low survival, high grade, and increased tumor size in breast cancer . Knockdown experiments of this protein have indeed promoted the epithelium-mesenchyme transition (EMT) via activation of the TGF-β pathway and resulted in a decrease in the tumor suppressor protein E-cadherin (CDH1) .
ARID1B (AT-rich interactive domain-containing protein 1B) is part of an SWI/SNF complex with among others ARID1A, its mutually exclusive homologue. This complex allows the remodeling of chromatin in ATP-dependent manner and plays an important role in cell proliferation, differentiation, development and tumor suppression. ARID1A is the third most significantly mutated gene in human rectal cancer, with a frequency of 39% in MSI type cancers [36,37,38]. It has also been shown that a decreased ARID1A expression results in inhibition of 5-FU-induced apoptosis in colorectal cancer cell lines . Moreover, when ARID1A is mutated, there is an increase in the probability of losing the expression of ARID1B . In our cohort, ARID1B is predominantly abundant in patients responding to 5-FU neoadjuvant chemoradiotherapy, which would imply the potential role of ARID1B in 5-FU resistance.
Finally, SLC25A33 (solute carrier family 25 member 33) is a mitochondrial transporter for pyrimidines and is essential for the metabolism of mitochondrial DNA and RNAs . Since 5-FU is an anti-pyrimidine analogue , it would be possible that SLC25A33 participates in the import of 5-FU into the mitochondria to carry out its cytotoxicity by incorporating itself into the mitochondrial DNA and RNAs . The absence of SLC25A33 in non-responder patients in our cohort raises the potential role of this protein in participating in 5-FU resistance mechanisms.
At the opposite, for NR, we identified seven proteins expressed in more than 75% of individuals in this response group and in less than 15% of TR: CALD1 (×4.4), CPA3 (×9.4), B3GALT5 (×55.9), CNNM4 (×51.3), MTIF3 (×4.1), CD177 (×2) and RIPK1 (×3.7).
Concerning CALD1 (caldesmon), this protein seems to display contradictory roles in cancer. On one side, CALD1 has been identified as a potential repressor of invasion into colorectal cancer cells . In contrast, CALD1 has been described as a key component of TGF-β-directed EMT by its overexpression . Recently, it was concluded that elevated expression of CALD1 in stromal cells of CRC patients was predictive of robustly shorter disease-free intervals .
For CPA3 (carboxypeptidase A3), it has been identified in a molecular signature of eleven genes for extracapsular spread in oral squamous cell carcinoma. This signature establishes a poor prognosis of outcome in patients without lymph node metastases . B3GALT5 (beta-1.3-galactosyltransferase 5) has been identified as a predictor of postoperative recurrence and survival in patients with hepatic carcinoma and its high expression has been associated with advanced stages and poor outcome . CD177 antigen has also been identified as a poor prognostic factor in patients with non-small cell lung cancer . For RIPK1 (receptor-interacting serine/threonine-protein kinase 1), this protein transduces inflammatory signals and cell death after binding of death receptors, activation of pathogen recognition receptors and DNA damage. The balance between its pro-death and pro-survival functions is very complex and still need to be discussed [50,51,52].
Another study substantially similar to ours has also focused on the identification of proteins predictive of the response to 5-FU neoadjuvant chemoradiotherapy in patients with rectal cancer . Using a 2D-DIGE quantitative proteomic approach, they analyzed the proteome of frozen biopsies and showed several proteins that could be predictive of patient nonresponse: fibrinogen β chain, three isoforms of actin, Serpin B9, peroxiredoxin 4 and two isoforms of cathepsin D. Although our patients have followed the same therapeutic regimen, we do not find the same proteins. This can be explained in several ways. (1) We used FFPE biopsies while they used frozen biopsies and their 2D-DIGE separation technique is also different from ours. (2) For patient grouping, we used the AJCC system (TRG0–TRG3) while Repetto and collaborators used the Mandard system (TRG1–TRG5) . Based on their grouping, they consider that TRG1 and TRG2 form the group of “good” response whereas if we rely on our grouping, we would have classified their TRG2 patients in the group of partial response and not in the total response. This shows the importance of specifying the starting material, the extraction technique and the patient classification system used.
Secondly, as indicated in the results, we also identified proteins involved in 5-FU metabolism in the differentially abundant proteins: DPYD and TYMP, found to be overexpressed in NR and responders, respectively. DPYD is the enzyme responsible for the hepatic elimination of 5-FU at approximately 80% . It allows the transformation of natural uracil and thymine pyrimidine bases into the inactive metabolites 5-FUH2 and 5-FTH2 (dihydrouracil and dihydrothymine), respectively . Patients with partial or total deficiency for this enzyme will be at high risk of toxicity to 5-FU treatment that may result in death . It has already been shown that high levels of DPYD mRNA correlate with 5-FU resistance  which is consistent with our results from the protein point of view. Regarding TYMP, also called PD-ECGF (platelet-derived endothelial cell growth factor), its role in 5-FU resistance is more controversial as reported by Elamin et al. . This protein stimulates metastasis, invasion, angiogenesis and cell death evasion on the one hand, and on the other hand catalyses the transformation reaction of 5-FU into FdUMP, a compound that will competitively inhibit TS, the main target of 5-FU . It has already been shown that high levels of TYMP mRNA correlate with 5-FU resistance . On the contrary, in our study cohort, it is the treatment-responsive patients group who have higher protein levels in TYMP compared to NR.
Thirdly, we were interested in the biological processes mainly represented in each group, thanks to the DAVID v6.8 software. What stands out most of our analyses are the ribosomal (RP) and mitochondrial (MP) proteins mostly represented in NRs. Among the RPs, it finds in particular the proteins composing the small (40S) and the large subunit (60S) of the ribosomes, RPS and RPL respectively. Among the MPs, the group most represented is that of NDUF proteins. These proteins encoded by nuclear DNA are part of the complex I (NADH dehydrogenase [ubiquinone]) of the mitochondrial respiratory chain. Interestingly, a study by Marin-Vicente et al.  on the effects of 5-FU on the proteome of colorectal cancer cell lines revealed that the predominantly down-regulated proteins after this treatment were MPs and RPs. They also showed that the effect on these proteins was not due to a general stress-induced cellular response since treatment with raltitrexed, an inhibitor of TS, did not produce the same effect . They concluded that the early effects on RPs and MPs were probably due to the effects of 5-FU on pre-rRNA maturation and mitochondrial ribosome biogenesis and/or tRNAs, respectively . In the presence of 5-FU, the pre-rRNAs do not transform into mature rRNAs. There is then an accumulation of pre-rRNAs that will be polyadenylated and degraded by the exosome. RPs will also accumulate in the nucleoplasm, interact with MDM2 and engage the 5-FU-mediated apoptotic response, and will be degraded. A higher representation of RPs in the NR patients’ tumor proteome in our cohort may be responsible for a lower response to 5-FU treatment. Regarding NDUF proteins, they are identified in a Zhang et al.  study on the proteogenomic characterization of colorectal cancer. This team has established signatures for the different proteomic subtypes established from many factors and NDUF proteins are found down-regulated within several of these signatures but not discussed . A higher representation of NDUF proteins within the tumor proteome of NR individuals in our cohort may be implicated in a lower sensitivity to 5-FU treatment.
In conclusion, our retrospective study aims to identify a predictive protein signature of the NRCT response to radiotherapy with 5-FU in patients with locally advanced non-metastatic rectal cancer from FFPE biopsies obtained in pre-treatment. Our results allowed us to validate our protein extraction protocol for the final realization of MS experiments. From the 23 patients in our cohort, we were able to establish three response groups, 6 NR, 10 PR and 7 TR, and compared the tumor proteomes. We have identified many differentially abundant proteins already identified as having a role in cancer, such as the earlier discussed proteins IFIT1, FASTKD2, PIP4K2B, ARID1B and SLC25A33 (more abundant in TR) or CALD1, CPA3, B3GALT5, CD177 and RIPK1 (more abundant in NR). We have also demonstrated proteins of the metabolism of 5-FU, DPYD and TYMP, and obtained results similar to those obtained in the literature in particular for DPYD, whose strong expression is associated with a poorer response to 5-FU. Finally, the detection of RPs and MPs seems to open the way for potential new mechanisms of resistance to 5-FU. By increasing the size of our cohort and moving to targeted proteomics experiments to detect and quantify the proteins of interest, we expect to define a protein signature to develop a test that could determine whether patients should be treated by RCTN with 5-FU.
AT-rich interactive domain-containing protein 1B
common terminology criteria for adverse events
database for annotation, visualization and integrated discovery
interferon regulatory factor 2-binding protein 2
double contrast barium enema
epidermal growth factor receptor
FAST kinase domain-containing protein 2, mitochondrial
faecal occult blood test
high performance liquid chromatography
interferon-induced protein with tetratricopeptide repeats 1
initial tumor, nodes, metastasis
magnetic resonance imaging
nuclear receptor-interacting factor 3
principal component analysis
platelet-derived endothelial cell growth factor (= TYMP)
phosphatidylinositol 5-phosphate 4-kinase type-2 beta
receptor-interacting serine/threonine-protein kinase 1
60S ribosomal protein, large subunit
40S ribosomal protein, small subunit
solute carrier family 25 member 33
total mesorectal excision
tumor, nodes, metastasis
transcripts per million
tumor regression grading system
thymidine phosphorylase (= PD-ECGF)
post-treatment tumor, nodes, metastasis
International Agency for Research on Cancer—World Health Organization. GLOBOCAN 2012: estimated cancer incidence, mortality and prevalence worldwide in 2012. http://globocan.iarc.fr.
Stintzing S. Management of colorectal cancer. F1000Prime Rep. 2014;6:108. https://doi.org/10.12703/P6-108.
Cellini F, Valentini V. Current perspectives on preoperative integrated treatments for locally advanced rectal cancer: a review of agreement and controversies. Oncology. 2012;26:730–5.
Gérard J, Conroy T, Bonnetain F, Bouché O, Chapet O, Closon-Dejardin M, et al. Preoperative radiotherapy with or without concurrent fluorouracil and leucovorin in T3–4 rectal cancers: results of FFCD 9203. J Clin Oncol. 2006;24:4620–5.
Sauer R, Liersch T, Merkel S, Fietkau R, Hohenberger W, Hess C, 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–33.
Sauer R, Becker H, Hohenberger W, Rödel C, Wittekind C, Fietkau R, et al. Preoperative versus postoperative chemoradiotherapy for rectal cancer. N Engl J Med. 2004;351:1731–40.
Bosset J, Collette L, Calais G, Mineur L, Maingon P, Radosevic-Jelic L, et al. Chemotherapy with preoperative radiotherapy in rectal cancer. N Engl J Med. 2006;355:1114–23.
Ortholan C, Romestaing P, Chapet O, Gerard J. Correlation in rectal cancer between clinical tumor response after neoadjuvant radiotherapy and sphincter or organ preservation: 10-year results of the Lyon R 96-02 randomized trial. Int J Radiat Oncol Biol Phys. 2012;83:e165–71.
Gerard J, Chapet O, Nemoz C, Hartweig J, Romestaing P, Coquard R, et al. Improved sphincter preservation in low rectal cancer with high-dose preoperative radiotherapy: the Lyon R96-02 randomized trial. J Clin Oncol. 2004;22:2404–9.
Roh MS, Colangelo LH, O’Connell MJ, Yothers G, Deutsch M, Allegra CJ, et al. Preoperative multimodality therapy improves disease-free survival in patients with carcinoma of the rectum: NSABP R-03. J Clin Oncol. 2009;27:5124–30.
Brændengen M, Tveit K, Berglund Å, Birkemeyer E, Frykholm G, Påhlman L, et al. Randomized phase III study comparing preoperative radiotherapy with chemoradiotherapy in nonresectable rectal cancer. J Clin Oncol. 2008;26:3687–94.
García-Flórez LJ, Gómez-Álvarez G, Frunza AM, Barneo-Serra L, Martínez-Alonso C, Fresno-Forcelledo MF. Predictive markers of response to neoadjuvant therapy in rectal cancer. J Surg Res. 2015;194:120–6.
Marshall KW, Mohr S, El Khettabi F, Nossova N, Chao S, Bao W, et al. A blood-based biomarker panel for stratifying current risk for colorectal cancer. Int J Cancer. 2010;126(5):1177–86.
Agostini M, Crotti S, Bedin C, Cecchin E, Maretto I, D’Angelo E, et al. Predictive response biomarkers in rectal cancer neoadjuvant treatment. Front Biosci. 2014;6:1–22.
Ryan JE, Warrier SK, Lynch AC, Ramsay RG, Phillips WA, Heriot AG. Predicting pathological complete response to neoadjuvant chemoradiotherapy in locally advanced rectal cancer: a systematic review. Color Dis. 2016;18:234–46.
Lièvre A, Blons H, Laurent-Puig P. Oncogenic mutations as predictive factors in colorectal cancer. Oncogene. 2010;29:3033–43.
Kim SH, Chang HJ, Kim DY, Park JW, Baek JY, Kim SY, et al. What is the ideal tumor regression grading system in rectal cancer patients after preoperative chemoradiotherapy? Cancer Res Treat. 2016;48:1–12.
Cox J, Mann M. MaxQuant enables high peptide identification rates, individualized p.p.b.-range mass accuracies and proteome-wide protein quantification. Nat Biotechnol. 2008;26:1367–72.
Tyanova S, Temu T, Sinitcyn P, Carlson A, Hein MY, Geiger T, et al. The Perseus computational platform for comprehensive analysis of (prote)omics data. Nat Methods. 2016;13:731–40.
Vizcaíno JA, Csordas A, Del-Toro N, Dianes JA, Griss J, Lavidas I, et al. 2016 update of the PRIDE database and its related tools. Nucleic Acids Res. 2016;44:D447–56.
Deutsch EW, Csordas A, Sun Z, Jarnuczak A, Perez-Riverol Y, Ternent T, et al. The ProteomeXchange consortium in 2017: Supporting the cultural change in proteomics public data deposition. Nucleic Acids Res. 2017;45:D1100–6.
R Core Team. R: the R project for statistical computing. R foundation for statistical computing, Vienna, Autria. [Internet]. 2017 [cited 2017 Nov 27]. https://www.r-project.org/.
Chisanga D, Keerthikumar S, Pathan M, Ariyaratne D, Kalra H, Boukouris S, et al. Colorectal cancer atlas: an integrative resource for genomic and proteomic annotations from colorectal cancer cell lines and tissues. Nucleic Acids Res. 2016;44:D969–74.
Huang DW, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009;37:1–13.
Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2008;4:44–57.
Jensen LJ, Kuhn M, Stark M, Chaffron S, Creevey C, Muller J, et al. STRING 8—a global view on proteins and their functional interactions in 630 organisms. Nucleic Acids Res. 2009;37:D412–6.
Hood BL, Conrads TP, Veenstra TD. Mass spectrometric analysis of formalin-fixed paraffin-embedded tissue: unlocking the proteome within. Proteomics. 2006;6:4106–14.
Sprung RW, Brock JWC, Tanksley JP, Li M, Washington MK, Slebos RJC, et al. Equivalence of protein inventories obtained from formalin-fixed paraffin-embedded and frozen tissue in multidimensional liquid chromatography-tandem mass spectrometry shotgun proteomic analysis. Mol Cell Proteom. 2009;8:1988–98.
Maes E, Broeckx V, Mertens I, Sagaert X, Prenen H, Landuyt B, et al. Analysis of the formalin-fixed paraffin-embedded tissue proteome: pitfalls, challenges, and future prospectives. Amino Acids. 2013;45:205–18.
Weichselbaum RR, Ishwaran H, Yoon T, Nuyten DSA, Baker SW, Khodarev N, et al. An interferon-related gene signature for DNA damage resistance is a predictive marker for chemotherapy and radiation for breast cancer. Proc Natl Acad Sci USA. 2008;105(47):18490–5.
Danish HH, Goyal S, Taunk NK, Wu H, Moran MS, Haffty BG. Interferon-induced protein with tetratricopeptide repeats 1 (IFIT1) as a prognostic marker for local control in T1-2 N0 breast cancer treated with breast-conserving surgery and radiation therapy (BCS + RT). Breast J. 2017;19:231–9.
Guo J, Hui DJ, Merrick WC, Sen GC. A new pathway of translational regulation mediated by eukaryotic initiation factor 3. EMBO J. 2000;19:6891–9.
Das S, Yeung KT, Mahajan MA, Samuels HH. Fas Activated Serine-Threonine Kinase Domains 2 (FASTKD2) mediates apoptosis of breast and prostate cancer cells through its novel FAST2 domain. BMC Cancer. 2014;14:852.
Yeung KT, Das S, Zhang J, Lomniczi A, Ojeda SR, Xu C, et al. A novel transcription complex that selectively modulates apoptosis of breast cancer cells through regulation of FASTKD2. Mol Cell Biol. 2011;31:2287–98.
Keune W, Sims AH, Jones DR, Bultsma Y, Lynch JT, Jirström K, et al. Low PIP4K2B expression in human breast tumors correlates with reduced patient survival : a role for PIP4K2B in the regulation of E-cadherin expression. Cancer Res. 2013;73:6913–25.
Mathur R, Alver BH, Roman AKS, Wilson BG, Wang X, Agoston AT, et al. ARID1A loss impairs enhancer-mediated gene regulation and drives colon cancer in mice. Nat Genet. 2017;49:296–302.
Cajuso T, Hänninen UA, Kondelin J, Gylfe AE, Tanskanen T, Katainen R, et al. Exome sequencing reveals frequent inactivating mutations in ARID1A, ARID1B, ARID2 and ARID4A in microsatellite unstable colorectal cancer. Int J Cancer. 2014;135:611–23.
Zhang B, Wang J, Wang X, Zhu J, Liu Q, Shi Z, et al. Proteogenomic characterization of human colon and rectal cancer. Nature. 2014;513:382–7.
Xie C, Fu L, Han Y, Li Q, Wang E. Decreased ARID1A expression facilitates cell proliferation and inhibits 5-fluorouracil-induced apoptosis in colorectal carcinoma. Tumor Biol. 2014;35:7921–7.
Helming K, Wang X, Wilson BG, Vazquez F, Haswell JR, Manchester HE, et al. ARID1B is a specific vulnerability in ARID1A-mutant cancers. Nat Med. 2014;10:730–6.
Di Noia MA, Todisco S, Cirigliano A, Rinaldi T, Agrimi G, Iacobazzi V, et al. The human SLC25A33 and SLC25A36 genes of solute carrier family 25 encode two mitochondrial pyrimidine nucleotide transporters. J Biol Chem. 2014;289:33137–48.
Heidelberger C, Chaudhuri NK, Danneberg P, Mooren D, Griesbach L, Duschinsky R, et al. Fluorinated pyrimidines, a new class of tumour-inhibitory compounds. Nature. 1957;179:663–6.
Parker W, Cheng Y. Metabolism and mechanism of action of 5-fluorouracil. Pharmacol Ther. 1990;48:381–95.
Yoshio T, Morita T, Kimura Y, Tsujii M, Hayashi N, Sobue K. Caldesmon suppresses cancer cell invasion by regulating podosome/invadopodium formation. FEBS Lett. 2007;581:3777–82.
Morita T, Mayanagi T, Sobue K. Dual roles of myocardin-related transcription factors in epithelial-mesenchymal transition via slug induction and actin remodeling. J Cell Biol. 2007;179:1027–42.
Calon A, Lonardo E, Berenguer-Llergo A, Espinet E, Hernando-Momblona X, Iglesias M, et al. Stromal gene expression defines poor-prognosis subtypes in colorectal cancer. Nat Genet. 2015;47:320–9.
Wang W, Lim WK, Leong HS, Chong FT, Lim TKH, Tan DSW, et al. An eleven gene molecular signature for extra-capsular spread in oral squamous cell carcinoma serves as a prognosticator of outcome in patients without nodal metastases. Oral Oncol. 2015;51:355–62.
Kuo HH, Lin RJ, Hung JT, Hsieh CB, Hung TH, Lo FY, et al. High expression FUT1 and B3GALT5 is an independent predictor of postoperative recurrence and survival in hepatocellular carcinoma. Sci Rep. 2017;7:1–13.
Sakabe T, Azumi J, Haruki T, Umekita Y, Nakamura H, Shiota G. CD117 expression is a predictive marker for poor prognosis in patients with non small cell lung cancer. Oncol Lett. 2017;13(5):3703–8.
Dondelinger Y, Jouan-Lanhouet S, Divert T, Theatre E, Bertin J, Gough PJ, et al. NF-κB-independent role of IKKα/IKKβ in preventing RIPK1 kinase-dependent apoptotic and necroptotic cell death during TNF signaling. Mol Cell. 2015;60:63–76.
Gentle IE, Wong WWL, Evans JM, Bankovacki A, Cook WD, Khan NR, et al. In TNF-stimulated cells, RIPK1 promotes cell survival by stabilizing TRAF2 and cIAP1, which limits induction of non-canonical NF-kB and activation of caspase-8. J Biol Chem. 2012;286:13282–91.
Green DR. Another face of RIPK 1. EMBO Rep. 2015;16:674–5.
Repetto O, De Re V, De Paoli A, Belluco C, Alessandrini L, Canzonieri V, et al. Identification of protein clusters predictive of tumor response in rectal cancer patients receiving neoadjuvant chemoradiotherapy. Oncotarget. 2017;8:28328–41.
Heggie GD, Sommadossi J-P, Cross DS, Huster WJ, Diasio RB. Clinical pharmacokinetics of 5-fluorouracil and its metabolites in plasma, urine, and bile. Cancer Res. 1987;47(8):2203–6.
Aparicio T, Ducreux M, Chaussade S. 5-fluorouracile : données sur le métabolisme et place actuelle dans le traitement des cancers digestifs. Gastroenterol Clin Biol. 2002;26:38–47.
Caudle KE, Thorn CF, Klein TE, Swen JJ, McLeod HL, Diasio RB, et al. Clinical pharmacogenetics implementation consortium guidelines for dihydropyrimidine dehydrogenase genotype and fluoropyrimidine dosing. Clin Pharmacol Ther. 2013;94:640–5.
Salonga D, Danenberg KD, Johnson M, Metzger R, Groshen S, Tsao-Wei DD, et al. Colorectal tumors responding to 5-fluorouracil have low gene expression levels of dihydropyrimidine dehydrogenase, thymidylate synthase, and thymidine phosphorylase. Clin Cancer Res. 2000;6:1322–7.
Elamin YY, Rafee S, Osman N, Byrne KJ, Gately K. Thymidine phosphorylase in cancer; enemy or friend? Cancer Microenviron. 2016;9:33–43.
Usuki K, Saras J, Waltenberger J, Miyazono K, Pierce G, Thomason A, et al. Platelet-derived endothelial cell growth factor has thymidine phosphorylase activity. Biochem Biophys Res Commun. 1992;184:1311–6.
Metzger R, Danenberg K, Leichman CG, Salonga D, Schwartz EL, Wadler S, et al. High basal level gene expression of thymidine phosphorylase (platelet-derived endothelial cell growth factor) in colorectal tumors is associated with nonresponse to 5-fluorouracil. Clin Cancer Res. 1998;4:2371–6.
Marin-Vicente C, Lyutvinskiy Y, Romans Fuertes P, Zubarev R, Visa N. The effects of 5-fluorouracil on the proteome of colon cancer cells. J Proteome Res. 2013;12:1969–79.
AC contributed to data acquisition and analysis as well as to drafting the manuscript. CSW contributed to selection of patients and obtention of their consent forms. SG and PGG contributed to selection of sample collection. AAM and VL participated to the development of the extraction protocol. FMB contributed to correction and final approval of the manuscript, and supervised the study. All authors read and approved the final manuscript.
This work was funded to F. M. B. by the Natural Sciences and Engineering Research Council of Canada (Grant Number 418404), and by a Université de Sherbrooke—Merck Sharp & Dohme grant. F. M. B. is also a FRQS Junior II scholar. Only the authors listed on the manuscript contributed towards the article.
The authors declare that they have no competing interests.
Availability of data and materials
All the datasets included in the present manuscript are presented in the main manuscript or additional supporting files. The complete data, including the raw files and all extracted peptides information can be found at the ProteomeXchange Consortium (http://proteomecentral.proteomexchange.org) via the PRIDE partner repository with the dataset identifier PXD008440.
Consent for publication
No individual personal data are included.
Ethics approval and consent to participate
We performed a post hoc analysis of a prospective cohort study including 23 rectal cancer patients who were hospitalized at the Centre Hospitalier Universitaire of Sherbrooke (CHUS—Québec, Canada) between September 2013 and September 2015. The study was conducted according to the principles expressed in the Declaration of Helsinki, and the study protocol was approved by the local ethics committee (Project #14-142). Before enrolment, informed consent was obtained from patients, their relatives, or their legal guardians.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.