- Research article
- Open Access
Molecular response to induction chemotherapy and its correlation with treatment outcome in head and neck cancer patients by means of NMR-based metabolomics
BMC Cancer volume 21, Article number: 410 (2021)
The aim of this prospective study is to identify the biomarkers associated with the effects of induction chemotherapy (iCHT) in terms of the favorable/weaker response to the treatment in locally advanced head and neck squamous cells carcinomas (LA-HNSCC).
The studied group consisted of 53 LA-HNSCC patients treated with iCHT. The treatment tolerance was measured by the Common Terminology Criteria for Adverse Events (CTCAE). The response to the treatment was evaluated by the clinical, fiberoptic and radiological examinations made before and after iCHT (the TNM Classification of Malignant Tumors was used for classifying the extent of cancer spread). Proton nuclear magnetic resonance (1H NMR) serum spectra of the samples collected before and after iCHT were acquired with a 400 MHz spectrometer and analyzed using the multivariate and univariate statistical methods.
The molecular response to iCHT involves an increase of the serum lipids which is accompanied by the simultaneous decrease of alanine, glucose and N-acetyl-glycoprotein (NAG). Furthermore, in males, the iCHT induced changes in the lipid signals and NAG significantly correlate with the regression of the primary tumor. The OPLS-DA multivariate model identified two subgroups of the patients with a weaker metabolic and clinical response. The first one consisted of the patients with a significantly lower initial nodal stage, the second one showed no differences in the initial clinical and metabolic statuses.
The NMR-based metabolomic study of the serum spectra revealed that iCHT induces the marked changes in the LA-HNSCC patients’ metabolic profiles and makes it possible to stratify the patients according to their response to iCHT. These effects are sex dependent. Further studies on a larger scale accounting for sex and the clinical and metabolic factors are warranted.
iCHT in the treatment of LA-HNSCC [1,2,3] is based on two standard regimens, which are the TPF schedule (docetaxel, cisplatin and 5-fluorouracil) and the PF schedule (cisplatin and 5-fluorouracil). TPF demonstrates a less favorable toxicity profile and a better response rate [4, 5]. Severe hematological adverse events (e.g. leukopenia, neutropenia, anemia) are more often reported in TPF [6,7,8], while non-hematological adverse events (e.g. mucositis, nausea, vomiting) seem to be higher in PF [6,7,8,9]. It is suggested that sex is among other factors, like age, performance status, disease status, lesion size, medication exposure and comorbidities, that may contribute, via the genetic, hormonal and molecular pathways, to the chemotherapy efficacy and its toxicity [10, 11].
Despite the improvements in treatment and new therapies, the survival rates in LA-HNSCC remain stagnant [12,13,14]. Thus, there is the urgent need for novel approaches that can be used to monitor both the tolerability and the treatment efficacy. Especially body fluid and liquid biopsy biomarkers appear to hold a strong potential for development the tools for such applications that could direct anticancer therapies in HNSCC . In recent years, there is an increasing number of the reports based on the metabolomic approaches attempting to monitor chemotherapy side effects and its efficacy in various cancer types [16,17,18,19,20,21]. In head and neck cancers Ye at al. revealed the important metabolic differences in the iCHT efficacy between the responsive and unresponsive patients using the chromatographic methods . They compared the pre- and post-chemotherapy serum samples and found that induction chemotherapy induces various metabolic changes, e.g. in carbohydrate metabolism, amino acid metabolism, fatty acid metabolism and steroid metabolism. Lactic acid, glutamic acid, and aspartic acid were even defined as the potential biomarkers of the iCHT efficacy . On the other hand, Jelonek et al.  analyzed the chromatographic profiles of the sera of the HNSCC patients treated with concurrent chemo-radiotherapy or radiotherapy alone with or without prior induction chemotherapy. They detected the differences in the molecular responses to the various treatment modalities, hypothetically associated with a different risk of toxicity, however no metabolic abundance changes between the pre- and post-treatment samples in case of iCHT were observed, presumably due to restricting their study only to the first iCHT cycle.
In the current work we decided to perform the NMR-based study of the influence of iCHT on the serum metabolome of the LA-HNSCC patients to identify the NMR biomarkers of the iCHT-induced metabolic changes and to assess their relationships with the response to treatment.
It may be hypothesized that such NMR approach may help to differentiate the well responding patients from the nonresponding ones, to explain the causes of the metabolic differences and to determine the influence of sex on the therapy response.
Characteristics of patient groups
The study was approved by Ethical Committee of the Maria Sklodowska-Curie National Research Institute of Oncology, and written informed consent was obtained from each patient. The retrospective study was conducted on the HNSCC patients treated oncologically in the 1st Radiation and Clinical Oncology Department of Maria Sklodowska-Curie National Research Institute of Oncology, Gliwice Branch. The study group consisted of 35 men and 18 women, all Caucasians, at the median age of 57 (22–74) years. Induction chemotherapy is mainly used in the advanced cancer stages (LA-HNSCC), therefore the patients in the current study where staged to TNM IV (81.1%) and III (18.9%). There were no patients with metastases (M0 = 100%), all patients were HPV negative.
The protocols (regimens) of iCHT were as follows:
3 cycles of TPF administered every 21 days (docetaxel 75 mg/m2, followed by cisplatin 100 mg/m2 on day 1, and 5-fluorouracil 1000 mg/m2 per day administered as a continuous 24-h infusion for 4 days) followed by chemoradiotherapy delivered as a sequential therapy - 21 patients 
4 cycles of TPF administered every 21 days (docetaxel 75 mg/m2 of body-surface area, followed by cisplatin 75 mg/m2 on day 1, 5-fluorouracil 750 mg/m2 per day administered as a continuous 24-h infusion for 4 days) followed by radiotherapy as a sequential therapy - 2 patients 
3 cycles of PF administered every 21 days (cisplatin 100 mg/m2 on day 1, 5-fluorouracil 1000 mg/m2 as a continuous 24-h infusion for 4 days) followed by chemoradiotherapy with cisplatin delivered as a sequential therapy – 28 patients [4, 23]
3 cycles of paclitaxel and carboplatin (PC) administered every 21 days (175 mg/m2) and carboplatin at a dose calculated using the Calvert formula area under the curve of 5 followed by chemoradiotherapy with carboplatin delivered as a sequential therapy – 2 patients .
Induction chemotherapy was administered every 3 weeks (defined as one cycle) for up to four cycles, unless progressive disease, unacceptable toxic effects, or withdrawal from the study occurred earlier. The patients received one (9 patients), two (6 patients), three (36 patients) or four iCHT cycles (2 patients). The patients with an incomplete iCHT were treated with the TPF (7 patients) and PF (8 patients) regimens. The administration of the antiemetics, antibiotics, steroids, hydration, diuretics was mandatory during each cycle.
In aim of quantifying the patients’ general well-being and the activities of daily life their performance status was measured to determine whether they can receive chemotherapy. All of them had 2 or less performance status by WHO/ECOG/Zubrod scale, the adequate hematologic, renal, and hepatic functions, without serious comorbidities and was able to receive iCHT. During chemotherapy, the patients were monitored clinically and with the laboratory tests on day 1 of each cycle, before starting the treatment. The iCHT timing depended on the treatment tolerance measured by the Common Terminology Criteria for Adverse Events (CTCAE) - the standardized classification of adverse effects of cancer therapy . CTCAE uses a range of grades from 1 to 5 (1 – Mild, 2 – Moderate. 3 – Severe, 4 - Life-threatening, 5 – Death). The most common causes of the delays in cycles’ administration were neutropenia or/and leukopenia grade 3 and more.
Assessment of treatment response to induction chemotherapy
The treatment response to induction chemotherapy was evaluated by the clinical, fiberoptic and radiological examinations. The TNM Classification of Malignant Tumors v. 7 was used before and after iCHT as a globally recognized standard for classifying the extent of cancer spread (T - tumor, N – regional involved lymph nodes, M - distant metastasis). The prefix modifiers were used as follows: “c” – a stage determined from the evidence acquired before the treatment (including the clinical examination, imaging, endoscopy, biopsy, surgical exploration), “y” – a stage assessed after iCHT in the same way as “c”. In our study most of the patients were the significant responders (yTNM lower than cTNM and an explicit percentage volumetric radiological remission) and iCHT was followed by radiotherapy or concurrent chemoradiotherapy in compliance with a schedule.
The primary lesions were measured within the head and neck in the images taken before (preCHT) and after induction chemotherapy (postCHT) in the radiology department of the Maria Sklodowska-Curie National Research Institute of Oncology. If an incomplete cycle of induction chemotherapy and immediate subsequent radical radiochemotherapy was performed (7 patients), the CT scan for the radiotherapy planning was used after induction chemotherapy for a sake of comparison. All types of imaging capable of obtaining a multiplanar reconstruction (MPR) were considered (CT, MR and PET-CT scans).
All volumetric measurements were performed by a qualified radiologist using the ‘MR Segmentation’ tool which is available with Siemens Healthcare syngo.via software, version VB20A. The following procedure was applied for each lesion. After the image was viewed with the ‘MR Basic’ system, the ‘MPR view’ was enabled. Then the ‘MR Segmentation tool’ was selected from the ‘VOI selection’ list to evaluate an arbitrary VOI in a semiautomatic way. The segmentation was performed using the ‘brush based segmentation’ tool. The middle part of the lesion was marked by tracing it with a green line, its boundary was enclosed by tracing it with a red line. The lesion area was defined in this way on every third layer. The starting and ending layers of the volume segmentation were marked with a solid red area only. Finally, the lesion volume was calculated, measured in cubic centimeters, and limited by a green outline.
The percentage of tumor regression was calculated using the following formula: 100-(postCHT tumor volume)*100/(preCHT tumor volume). The clinical response to iCHT was expressed as a difference of cT-yT and cN-yN.
Serum samples collection
The overnight fasting blood samples from the peripheral vein were collected before (preCHT) the first cycle of iCHT and within a week preceding the start of radiation therapy (postCHT), at the median of 45 days after the end of iCHT. The blood samples were incubated for 30 min at room temperature and then centrifuged (1000×g, 10 min) to remove a clot, and stored frozen at − 80 °C until the NMR measurements are performed. The total number of 130 samples were collected.
Sample preparation for NMR spectroscopy
The serum samples were thawed in two steps (at 4 °C and at room temperature) and mixed with the phosphate buffer (pH 7.4) containing D2O and TSP. The aliquots of 600 μl of the solution were transferred into 5 mm Wilmad WG-1235-7 NMR tubes (Wilmad Labglass, USA) and kept at 4 °C until the NMR analysis.
The same measurement protocol as in our previous metabolomic studies [26, 27] was applied. The 1H NMR spectra were acquired on a Bruker 400 MHz Avance III spectrometer (Bruker Biospin, Rheinstetten, Germany) equipped with a 5 mm PABBI probe. The quality control tests were performed at every measurement day. The NMR probe tuning and matching, shimming, determination of the transmitter offset value for the water pulse presaturation and 90° pulse adjustments were always made for each sample. The receiver gain was set to 90.5 and the temperature to 310 K for all experiments. Four different 1H NMR spectra – NOESY (1D nuclear Overhauser enhancement spectroscopy), CPMG (Carr-Purcell-Meiboom-Gill), diffusion edited (DIFF) and J-resolved (JRES) – were acquired for each serum sample. The characteristics of the acquired spectra as well as the pulse sequence parameters are given in Supplementary Material (Table S1).
All 1D spectra were processed with a line broadening of 0.3 Hz and automatically phase corrected (in Topspin software from Bruker Biospin), referenced to the methyl doublet of alanine at 1.5 ppm and bucketed over the region 9.0–0.5 ppm with the bucket width set to 0.002 ppm using AMIX software (Bruker Biospin). The spectrum region of water (5.15–4.38 ppm, d = 0.77 ppm) was removed from the analysis in order to prevent variation in each sample. No normalization was applied. This is the standard processing protocol used in our metabolomic lab [26, 27].
The metabolites identification was done based on the comparisons with the reference compounds library (in Chenomx NMR Suite Professional (Chenomx Inc., Edmonton, Canada)), as well as on the multiplicity and scalar couplings information extracted from the 2D JRES spectra, and using the information from Human Metabolome Database (http://www.hmdb.ca/) and the available literature.
The low molecular weight metabolites were quantified based on the 1D positive projections of the JRES spectra. The diffusion edited spectra were used for quantification of the lipid signals. The integrals were measured in the spectral regions defined individually for each metabolite using the “sum all points in region” method in AMIX (Bruker Biospin) software.
Data analysis and validation of multivariate model
The data analysis was carried out using SIMCA-P+ (Umetrics, v. 14) software. The post-processed NMR spectra were Pareto scaled and the distinction between the preCHT and postCHT groups was done using the orthogonal partial least squares discriminant analysis (OPLS-DA). The OPLS-DA results are presented graphically in two types of plots, which detailed description and interpretation is available in . The validation of the OPLS-DA model was carried out using the internal cross-validation.
The statistical significance of the estimated predictive power of the OPLS-DA model was tested using permutation testing and ANOVA of the cross-validated residuals (cv-ANOVA) as well as the external test set consisted of 24 serum samples acquired from 24 (LA-HNSCC) patients. The blood samples were collected after iCHT (23 patients) and in case of one patient only the pre-treatment sample was included (the post-treatment one was excluded due to an insufficient blood sample volume). The detailed description of the external test group is given in Supplementary Material (Table S2).
The metabolites’ integrals were evaluated for their statistical significance with the Wilcoxon signed-rank (WSR), Mann-Whitney U (MWU), Kruskal-Wallis ANOVA (KWA) and Spearman’s rank correlations (SRC) tests. The significance threshold was set at 0.05. Statistica software (Statsoft, v. 12) was applied for the univariate statistics.
The detailed baseline characteristics of the patients under the investigation is presented in Table 1.
The 400 MHZ 1H NMR CPMG median spectra from the preCHT and postCHT groups are shown in Fig.S1.
Distinction between the preCHT and postCHT samples using OPLS-DA
The results from the OPLS-DA analysis are presented in Fig. 1. The OPLS-DA scores plot (Fig. 1a) shows a fair separation between the pre- and postCHT classes, with a small number of the mixing samples. The metabolites important for a class separation (lipids, alanine and glucose) are identified based on the s-line plot (Fig. 1c). The summary of the OPLS-DA analysis is presented in Table 2. The results of the validation procedure of the OPLS-DA model are presented in Fig. 1b (the predicted (tPS) scores plot for the external test set), Fig. S2 (the permutation plot) and Table 2 (the misclassification table for the external test set and cv-ANOVA).
Analysis of the mixing samples – investigation of demographic and clinical factors
The OPLS-DA model is able to distinguish the preCHT and postCHT samples, although the separation is not very distinct due to the cloud of the mixing samples seen in the OPLS-DA scores plot center (Fig. 1a). These mixing samples are shown in Fig. 2, a modified version of Fig. 1a, as two groups marked in red and green (if for a particular patient a postCHT sample is identified as a mixing one, the corresponding preCHT sample from this patient is also colored), whereas the non-mixing ones are greyed out (the Grey group). We identified 10 postCHT cases shifted to the left side of the plot, towards the negative values (the area occupied by the preCHT samples), and 9 preCHT cases shifted to the right side, towards the positive values (the area occupied by the postCHT samples). The mixing postCHT samples and their corresponding preCHT samples shifted to the left side of the plot are colored in red (the Red group), while the mixing preCHT samples and their corresponding postCHT samples shifted to the right side of the plot are colored in green (the Green group). Such visualization is helpful to trace the trajectories of the iCHT induced changes in the serum metabolic profiles. As seen, for both Green and Red groups the weaker metabolic response, measured as a distance between the adequate preCHT-postCHT t  coordinates, is observed than for the Grey, non-mixing group. This distance is the largest for the Grey group (KWA with the p-values of 0.027 and 0.034 for the comparison with the Red and Green ones) indicating the largest amplitude of the metabolic alterations in the grey patients. Moreover, except for one green patient (marked in Fig. 2 by the arrow), the trajectories are correct, i.e. directed from left (preCHT) to right (postCHT).
The box plot representation of the relative changes in the metabolites identified by OPLS-DA and the WSR test as important to differentiate the preCHT and postCHT samples in the mixing and non-mixing patients is presented in Fig. 3. The lipid signals are significantly increased after iCHT only in the non-mixing patients (p < 0.000, the Grey group, WSR). The Green patients start and end the treatment with the higher lipid levels as compared to the values detected for the other groups: the preCHT signal at 5.3 ppm is significantly higher (p = 0.04, KWA) than for the grey, non-mixing patients, whereas the postCHT lipid signals at 0.9 and 5.3 ppm are significantly higher than in the Red group (p = 0.011 and 0.011 respectively, KWA), however the treatment-related increase in this group is negligible. Similarly, glucose is significantly decreased after iCHT only in the non-mixing patients (p = 0.03, WSR). The decrease of alanine is significant in the Grey (p < 0.000, WSR) and Red (p = 0.03, WSR) groups, while the NAG is significantly decreased in the Grey (p = 0.005, WSR) and Green (p = 0.04, WSR) ones.
The baseline demographic and clinical status of the mixing patients was carefully investigated. The Green patients had significantly lower initial nodal (cN) stage when compared to the Red and the non-mixing (Grey) ones (KWA with the p-values of 0.029 and 0.047, respectively). There were no significant differences in the cT and cTNM stages, tumor sites, sex, age. The treatment related factors were also evaluated, showing no significant differences in the number of the iCHT cycles, overall treatment time (OTT) as well as utilization of the particular iCHT scheme (TPF/PF) between the Grey, Red and Green subgroups. Assessment of the iCHT induced toxicity revealed mostly the grade 1 or 2 (low-grade) toxicity according to the CTCAE. The most common low-grade hematological adverse events were: anemia, leukopenia, thrombocytopenia and neutropenia affecting 55, 34, 23 and 15% of the patients, respectively. The high grade (grade 3 or 4) adverse events were observed in 26% of the patients (neutropenia), 15% of the patients (leukopenia) and 4% of the patients (thrombocytopenia). The intensification of the side effects were at the same levels (KWA the p-values close to 1) within the subgroups.
Influence of baseline demographic and clinical factors on the metabolic profile
The pretreatment lipids at 0.9, 1.3 and 5.3 ppm showed a weak (SCR R = -0.39, − 0.35 and − 0.36, respectively) inverse correlation with the initial TNM stage, while the pretreatment alanine correlated weakly with the initial nodal stage (R = 0.37). Furthermore, the preCHT levels of glucose and alanine were significantly (MWU p-values of 0.009 and 0.001, respectively) higher in males (M) compared to females (F), and thus in TPF compared to PF due to the different M/F ratios in these groups (there were fewer women treated with TPF than with PF, although the difference was not statistically significant (p = 0.15, MWU)).
The box plot representation of the relative changes in the metabolites identified by OPLS-DA and the WSR test as important to differentiate the preCHT and postCHT samples in M and F is presented in Fig. 4. It is clearly visible that only in M the levels of glucose and alanine are lower after iCHT. There are no differences in the glucose and alanine levels between M and F when comparing the postCHT samples.
Correlation of the tumor regression after iCHT with the metabolic profile
In our study, most of the patients (87%) were the significant responders (as defined in Assessment of treatment response to induction chemotherapy), however, the degree of tumor shrinkage varied among the cases. The values of the radiologically evaluated tumor volumes, the percentages of the primary tumor regression and the initial and post-treatment clinical staging are provided in Table S3. Although the percentage of the tumor regression was higher in the Grey group compared to Red and Green, the difference was below the statistical significance. The percentage of the patients without the postCHT primary tumor and the nodal downstaging were markedly higher in Red (20%) and Green (22%) compared to Grey (8.8%), however, due to the much larger size of the Grey group than the Red and Green groups, this cannot be the basis for a statistical inference.
In order to evaluate whether the magnitude of the postCHT changes in the analyzed metabolites is correlated with a clinically or radiologically measured response to the treatment (as defined in Assessment of treatment response to induction chemotherapy), the SRC coefficients were calculated. The magnitude of the metabolic alterations was expressed as the difference between the preCHT and postCHT levels of the metabolites. The correlations were calculated for the whole group as well as separately, for the subgroups stratified according to the sex-related and metabolic responses (Grey, Red, Green). The obtained results (only with a statistical significance, i.e. with p < 0.05 and R > 3) are listed in Table 3.
For the whole studied group, the clinically observed regression of the primary tumor (cT-yT) is weakly correlated with the postCHT increase of the lipid signals (R = -0.31, SRC) and the decrease of NAG (R = 0.33, SRC). In M this correlation is stronger (|R| ~ 0.4–0.57, SRC), while in F it is absent. The increase in lipids is also correlated with cT-yT in the Grey group (R = -0.38 and − 0.46, SRC), while in the Green group, the cT-yT correlates with the postCHT decrease of NAG (R = 0.73, SRC) and alanine (R = 0.7, SRC) only. The percentage of the primary tumor regression calculated based on the radiological scans is worse explained by the metabolic profile showing only the correlation with the postCHT increase in lipids (M) (R = -0.45 and − 0.4, SRC) and the alanine decrease (Green) (R = 0.77, SRC). The clinically observed regression of the nodal stage, cN-yN, shows no correlation with the important metabolites. Furthermore, no metabolic correlations between the response to the treatment in the Red group as well as in F is observed.
In our previous works we identified the metabolic changes in the blood serum related to the toxicity of radio- and chemoradiotherapy in the HNSCC patients [26, 27], however, the potential influence of induction chemotherapy was not assessed. In the present work we examine the serum metabolic profiles acquired before the induction chemotherapy (preCHT) and within a week before the start of a sequential radio- or chemoradiotherapy (postCHT) and analyze the metabolic changes in correlation with the conditions related to patient (age, sex), treatment (OTT, iCHT regimen, toxicity) and oncological disease (the TNM re-stage and volumetric radiological remission after treatment).
The already gathered “omic” results suggest that to expand our understanding of the iCHT effects as well as of the possible mechanisms of the toxicity induced by chemoradiotherapy and of a subsequent radiation therapy new studies should focus on the molecular mechanisms and biochemical pathways [26,27,28]. As reveals from our present study, the NMR-based metabolomics sheds some light on the metabolic impact and efficacy of induction chemotherapy. To the best of our knowledge, this is the first NMR based metabolomics approach to the identification of the biomarkers of induction chemotherapy in HNSCC.
OPLS-DA mixing serum samples – the measure of response to iCHT
The OPLS-DA model with two orthogonal components successfully separated the preCHT and postCHT serum samples (Fig. 1a) as well as was able to predict the sample class from the external test group (Fig. 1b). However, there were 19 identified mixing samples located in the origin of the OPLS-DA scores plot (Fig. 1a). These mixing – Red and Green samples (with the corresponding pre/postCHT samples, 19 pairs) – were carefully checked for any clinical, demographic or treatment related conditions that may affect their metabolic profiles. The Green patients had significantly lower initial nodal stage (cN), and presumably this is a reason of their spectra to be shifted to the right in the OPLS-DA scores plot (the postCHT direction). This statement can be supported by the fact that in the study group the pretreatment lipids are inversely correlated with the initial TNM stage, while the pretreatment alanine correlates with the initial nodal stage. The Green patients had the higher levels of the lipid compounds as well as the lower level of alanine before the treatment when compared to others (Fig. 3). The iCHT-induced metabolic changes in the Green group show the correct (except for one patient) trajectory, i.e. from left to right (Fig. 2), however the response is significantly weaker than in the Grey group. The median percentage of the primary tumor regression is also lower in this group (71.8%) compared to Grey (87%), although without a statistical significance (Table S3).
In the Red subgroup, a weaker (compared to Grey) metabolic response to iCHT was also observed with only a significant decrease in the alanine postCHT level (Fig. 3) and the median percentage of the primary tumor regression of 76.4% (Table S3). However, contrary to the Green group, these patients do not differ significantly at baseline (preCHT) from the Grey group. Their glucose, alanine, NAG and lipids levels are similar to those for Grey and they do not present any differences in sex, age and the initial tumor stage.
S-line plot – the metabolites important for discrimination between the preCHT and postCHT classes
The metabolites important for discrimination between the preCHT and postCHT classes were identified based on the OPLS-DA s-line plot (Fig. 1c). The characteristic metabolic changes for the postCHT samples are the increased lipid signals with the simultaneous decreases in alanine, glucose and NAG.
The effect of chemotherapy on the serum lipids has been studied for at least three decades [20, 29,30,31,32], mainly in breast cancer [33, 34]. To the best of our knowledge no such studies has been conducted on HNSCC patients, despite the fact that the lipid profile is altered in this malignancy  as well as the main chemotherapeutic agents for HNSCC (cisplatin, docetaxel and 5-fluorouracil) are reported to induce dyslipidemia [30, 32, 33]. There is no consensus in the literature about the character of the chemotherapy-induced changes in the lipid profile. However, an increase of triglycerides (TG) and low-density lipoprotein (LDL) cholesterol with a simultaneous reduction of high-density lipoprotein (HDL) cholesterol is often reported in breast cancer [29, 33, 34], wherein an increase in TG could be a potential biomarker of effective chemotherapy . In colorectal cancer chemotherapy induces an increase in total cholesterol, TG, HDL and lowers the levels of LDL . In the patients who favorably responded to the treatment the increased levels of the blood phospholipids were detected  and those with a significant chemotherapy induced HDL elevation exhibit a better 3-year disease free survival and 7-year overall survival than those without . In other types of cancers an increase of total cholesterol and LDL cholesterol was also observed in the patients with a significant response to the treatment .
The standard NMR acquisition protocol for the measurement of the serum lipid signals (diffusion edited) and the NOESY protocol provide information about the mobile unsaturated lipids (the lipid signal at 5.3 ppm) as well as the total cholesterol (the lipid signals at 0.9 and 1.3 ppm), however without an ability to identify the particular lipoprotein subclasses, i.e. VLDL (very low density lipoprotein), LDL, HDL, etc. . Nevertheless, these NMR experiments are suitable to study a metabolic response to iCHT.
Our results show a close relationship between the increase of serum lipids and a favorable response to iCHT – the postCHT increase in the lipid signals correlates with the clinical stage reduction of the primary tumor (Table 3). Such relation is not, however, observed in F as well as in the patients with a weaker metabolic response to iCHT (the Red and Green groups). In the Red group with the lowest increase of lipids after iCHT one patient discontinued chemotherapy due to myocardial infarction, one patient had a complete pathological response and eight patients ended iCHT with a partial regression. The median degree of the tumor regression in this group was slightly lower as compared to the patients who had significantly elevated lipids (Table S3). Also, in the Green group the increase of the lipid signals was insignificant, but the green patients showed the lowest median percentage of the primary tumor regression, however, without a statistical significance. For the Grey group the appropriate regression percentage was 87.01%. In the whole study group 7 patients (6 women, 1 man) showed no postCHT reduction of the clinical tumor and nodal staging (cT = yT and cN = yN), in addition two of them showed no radiological percentage of the tumor regression after iCHT (Table S3). However, there was no specific metabolic fingerprint differentiating them from others (even when the non-responding female patients where compared to those who favorably responded).
Another metabolic effect of induction chemotherapy is lowering of the levels of glucose and alanine, two metabolites which blood levels are significantly positively correlated with each other . However, their reduction was observed only in the male patients and is due to the increased pretreatment levels of these metabolites as compared to the F group.
Possible explanation of the sex-related metabolic changes in the post iCHT HNSCC patients
The detailed investigation of the available literature data gives few possible explanations of these sex-related differences: 1) Fasting blood glucose among males is found to be higher than in females ; 2) Impaired fasting glucose in pre-diabetes and diabetes is more prevalent in men compared to women ; 3) Glucose metabolism changes during the development and progression of oral and tongue squamous cell carcinomas – the serum glucose is elevated in oral cancer patients (compared to controls) and is remarkably increased in patients with a late stage disease than in those with an early stage . In our study group, there were 3 confirmed diabetic patients (only the male ones, one with a type-1 diabetes; no information about the pre-diabetic status in the remaining patients), however their preCHT glucose and alanine levels were not higher as compared to the remaining male patients. In addition, after excluding the diabetics, the difference in the preCHT glucose and alanine levels between the M and F groups as well as the postCHT drop of glucose and alanine in M were still significant. Considering that men are more reluctant to go to the doctor , thus their disease lasts longer (is more advanced, i.e. more men with HNSCC present the nodal stages 2 and 3 as well as the TNM stage IV than women) , the metabolic alterations due to cancer may be more pronounced. It seems to be reflected in the structure of the study group where more men than women presented higher cN stage (28.3% of the male patients compared to 5.6% of the female patients presented cN3 stage and 14.3% of men compared to 33.4% of women presented cN stage of 0 or 1), with a difference close to a statistical significance (p = 0.061, MWU). In consequence, although the influence of sex and diabetic status cannot be neglected, it seems likely that the elevated preCHT glucose and alanine levels in M could be due to the severity of the disease. On the other hand, the decrease of glucose and alanine in the male patients is not correlated with the clinical or radiological response to the treatment (Table 3). However, such correlation is observed in one of the subgroups with an insignificant postCHT increase in lipids – in the Green group the decrease in alanine is significantly correlated with a favorable response to iCHT (Table 3). This is particularly surprising, because the green patients started the treatment with the lowest level of alanine.
Furthermore, glucose and alanine are the biologically correlated metabolites, i.e. alanine is one of the major substrates for gluconeogenesis during the increased demand for energy. Both these metabolites are strongly related to the patient’s nutritional status and their decreased levels in the blood serum correlate with the weight loss during the anticancer treatment . In the study group, despite the higher nodal stage, in M the levels of BMI and prealbumin are higher than in F (p = 0.022 and 0.003 respectively, MWU). However, contrary to the NMR detected metabolites, BMI, prealbumin and albumin do not seem to be affected by iCHT, showing no differences between the preCHT and postCHT samples for both M and F (p ≅ 0.5, WSR). Thus, we can assume that the observed metabolic alterations are probably not due to the differences in the nutritional status.
Inflammation marker vs. sex
NAG is an NMR marker of an inflammatory state  and it correlates well with the level of C-reactive protein (a clinical marker of inflammation) [26, 27], however in the study group such correlation was observed only for M (p < 0.05, R = 0.57 and 0.33 for preCHT and postCHT samples respectively, SRC). Despite the more prevalent higher N stage in the male patients, no differences were observed in the preCHT NAG between the M and F groups. An inflammation reduction was observed in both sexes, however, such effect was significant only in M (p = 0.001, WSR). Similarly to glucose and alanine, it may, however, be related to the initial tumor stage. Furthermore, the decrease of NAG after iCHT correlates with a favorable response to the treatment, although not in the F group (Table 3). We are not able to find the explanation for the different response of the inflammatory marker between M and F.
NMR-based metabolomics can stratify the patients according to their response to iCHT. Further studies on a larger scale accounting for the sex and the clinical and metabolic factors are warranted.
The study limitations
The findings of this study have to be seen in light of some limitations already mentioned earlier. The main is the small sizes of the studied groups. The importance of this limitation is all the more important as the groups are heterogeneous in terms of sex, which turns out to influence the observed effects. Moreover, the percentages of the patients in the Red, Green and Grey groups differ, which makes the statistical inference difficult. The remedy is increasing the groups – it will also significantly facilitate the NMR detection of the metabolites of low intensity in the serum NMR spectra, which in turn will help to determine the metabolic paths responsible for the observed effects. Undoubtedly, further studies should be performed on a larger scale and accounting for sex and the clinical and metabolic factors.
The NMR-based metabolomic study of the serum spectra revealed that iCHT induces the marked changes in the LA-HNSCC patients’ metabolic profiles. The molecular response to iCHT involves an increase of the serum lipids which is accompanied, however only in the males, by the simultaneous decreases of alanine, glucose and N-acetyl-glycoprotein (NAG). Furthermore, in M, the iCHT induced changes in the lipid signals and NAG significantly correlate with the primary tumor regression.
The metabolomic approach makes it also possible to stratify the patients according to their response to iCHT. In the weaker responding group two sub-groups were detected: the first one with a significantly lower initial nodal stage and the second one showing no differences in the initial clinical and metabolic statuses. These effects are sex dependent. Further studies on a larger scale accounting for the sex and the clinical and metabolic factors are warranted.
Availability of data and materials
The datasets generated during and/or analyzed during the current study are available from the corresponding author on a reasonable request.
- cN, cT:
nodal (cN) and primary tumor (cT) stage before induction chemotherapy
common terminology criteria for adverse events
Eastern Cooperative Oncology Group
high density lipoprotein
human papilloma virus
locally advanced head and neck squamous cell carcinoma
low density lipoprotein
Mann-Whitney U test
nuclear magnetic resonance
nuclear overhauser effect spectroscopy
orthogonal partial least squares discriminant analysis
overall treatment time
chemotherapy schedule based on paclitaxel and carboplatin
chemotherapy schedule based on cisplatin and 5-fluorouracil
parts per million
after induction chemotherapy
before induction chemotherapy
a predictive ability of the OPLS-DA model
an amount of variation in the data that is correlated to class separation
an amount of variation in the data that is uncorrelated (orthogonal) to the class separation
a fraction of the class membership (Y) variation modeled using the data matrix (X)
Spearman’s rank correlation
classification of malignant tumors
chemotherapy schedule based on docetaxel, cisplatin and 5-fluorouracil
predicted OPLS-DA score
very low density lipoprotein
World Health Organization
Wilcoxon signed rank test.
Fayette J, et al. Neoadjuvant modified TPF (docetaxel, cisplatin, fluorouracil) for patients unfit to standard TPF in locally advanced head and neck squamous cell carcinoma: a study of 48 patients. Oncotarget. 2016. https://doi.org/10.18632/oncotarget.8934.
Haddad RI, Posner M, Hitt R, Cohen EEW, Schulten J, Lefebvre JL, et al. Induction chemotherapy in locally advanced squamous cell carcinoma of the head and neck: role, controversy, and future directions. Ann Oncol. 2018;29(5):1130–40. https://doi.org/10.1093/annonc/mdy102.
de Bree R, Wolf GT, de Keizer B, Nixon IJ, Hartl DM, Forastiere AA, et al. Response assessment after induction chemotherapy for head and neck squamous cell carcinoma: from physical examination to modern imaging techniques and beyond. Head Neck. 2017;39(11):2329–49. https://doi.org/10.1002/hed.24883.
Posner MR, Hershock DM, Blajman CR, Mickiewicz E, Winquist E, Gorbounova V, et al. Cisplatin and fluorouracil alone or with docetaxel in head and neck cancer. N Engl J Med. 2007;357(17):1705–15. https://doi.org/10.1056/NEJMoa070956.
Vermorken JB, Remenar E, van Herpen C, Gorlia T, Mesia R, Degardin M, et al. Cisplatin, fluorouracil, and docetaxel in unresectable head and neck cancer. N Engl J Med. 2007;357(17):1695–704. https://doi.org/10.1056/NEJMoa071028.
Albers AE, Grabow R, Qian X, Jumah MD, Hofmann VM, Krannich A, et al. Efficacy and toxicity of docetaxel combination chemotherapy for advanced squamous cell cancer of the head and neck. Mol Clin Oncol. 2017;7(1):151–7. https://doi.org/10.3892/mco.2017.1281.
Rapidis A, Sarlis N, Lefebvre JL, Kies M. Docetaxel in the treatment of squamous cell carcinoma of the head and neck. Ther Clin Risk Manag. 2008;4(5):865–86. https://doi.org/10.2147/TCRM.S3133.
Jin T, Qin WF, Jiang F, Jin QF, Wei QC, Jia YS, et al. Cisplatin and fluorouracil induction chemotherapy with or without Docetaxel in Locoregionally advanced nasopharyngeal carcinoma. Transl Oncol. 2019;12(4):633–9. https://doi.org/10.1016/j.tranon.2019.01.002.
Galizia D, Milani A, Geuna E, Martinello R, Cagnazzo C, Foresto M, et al. Self-evaluation of duration of adjuvant chemotherapy side effects in breast cancer patients: a prospective study. Cancer Med. 2018;7(9):4339–44. https://doi.org/10.1002/cam4.1687.
Kim HI, Lim H, Moon A. Sex differences in cancer: epidemiology, genetics and therapy Biomol. Ther. 2018;26(4):335–42.
Wang J, Huang Y. Pharmacogenomics of sex difference in chemotherapeutic toxicity. Curr Drug Discov Technol. 2007;4(1):59–68. https://doi.org/10.2174/157016307781115485.
Masuda M, Toh S, Wakasaki T, Suzui M, Joe AK. Somatic evolution of head and neck cancer—biological robustness and latent vulnerability. Mol Oncol. 2013;7(1):14–28. https://doi.org/10.1016/j.molonc.2012.10.009.
Tribius S, Reemts E, Prosch C, Raguse M, Petersen C, Kruell A, et al. Global quality of life during the acute toxicity phase of multimodality treatment for patients with head and neck cancer: can we identify patients most at risk of profound quality of life decline? Oral Oncol. 2012;48(9):898–904. https://doi.org/10.1016/j.oraloncology.2012.03.011.
Santa Cruz O, Tsoutsou P, Castella C, Khanfir K, Anchisi S, Bouayed S, et al. Locoregional control and toxicity in head and neck carcinoma patients following helical Tomotherapy-delivered intensity-modulated radiation therapy compared with 3D-CRT data. Oncology. 2018;95(2):61–8. https://doi.org/10.1159/000489217.
Hsieh JC-H, et al. Review of emerging biomarkers in head and neck squamous cell carcinoma in the era of immunotherapy and targeted therapy. Head Neck. 2019;41(S1):19–45. https://doi.org/10.1002/hed.25932.
Miolo G, et al. Pharmacometabolomics study identifies circulating spermidine and tryptophan as potential biomarkers associated with the complete pathological response to trastuzumab-paclitaxel neoadjuvant therapy in HER-2 positive breast cancer. Oncotarget. 2016. https://doi.org/10.18632/oncotarget.9489.
Wei S, Liu L, Zhang J, Bowers J, Gowda GAN, Seeger H, et al. Metabolomics approach for predicting response to neoadjuvant chemotherapy for breast cancer. Mol Oncol. 2013;7(3):297–307. https://doi.org/10.1016/j.molonc.2012.10.003.
Sun Y, Kim JH, Vangipuram K, Hayes DF, Smith EML, Yeomans L, et al. Pharmacometabolomics reveals a role for histidine, phenylalanine, and threonine in the development of paclitaxel-induced peripheral neuropathy. Breast Cancer Res Treat. 2018;171(3):657–66. https://doi.org/10.1007/s10549-018-4862-3.
Jobard E, Blanc E, Négrier S, Escudier B, Gravis G, Chevreau C, et al. A serum metabolomic fingerprint of bevacizumab and temsirolimus combination as first-line treatment of metastatic renal cell carcinoma. Br J Cancer. 2015;113(8):1148–57. https://doi.org/10.1038/bjc.2015.322.
Yang K, Zhang F, Han P, Wang ZZ, Deng K, Zhang YY, et al. Metabolomics approach for predicting response to neoadjuvant chemotherapy for colorectal cancer. Metabolomics. 2018;14(9):110. https://doi.org/10.1007/s11306-018-1406-0.
Ye G, Liu Y, Yin P, Zeng Z, Huang Q, Kong H, et al. Study of induction chemotherapy efficacy in oral squamous cell carcinoma using pseudotargeted metabolomics. J Proteome Res. 2014;13(4):1994–2004. https://doi.org/10.1021/pr4011298.
Jelonek K, Krzywon A, Jablonska P, Slominska EM, Smolenski RT, Polanska J, et al. Systemic effects of radiotherapy and concurrent chemo-radiotherapy in head and neck Cancer patients—comparison of serum Metabolome profiles. Metabolites. 2020;10(2). https://doi.org/10.3390/metabo10020060.
Taylor S, et al. Randomized comparison of neoadjuvant cisplatin and fluorouracil infusion followed by radiation versus concomitant treatment in advanced head and neck cancer. J Clin Oncol. 1994;12(2):385–95. https://doi.org/10.1200/JCO.19220.127.116.115.
Dunphy FR, Dunleavy TL, Harrison BR, Trinkaus KM, Kim HJ, Stack BC, et al. Induction paclitaxel and carboplatin for patients with head and neck carcinoma. Analysis of 62 patients treated between 1994 and 1999. Cancer. 2001;91(5):940–8. https://doi.org/10.1002/1097-0142(20010301)91:5<940::AID-CNCR1083>3.0.CO;2-A.
http://evs.nci.nih.gov/ftp1/CTCAE/About.html Accessed 27 October 2020.
Boguszewicz, Ł., Hajduk A., Mrochem-Kwarciak J., Skorupa A., Ciszek M., Heyda A., Składowski K., Sokół M. 1H NMR based metabolomic approach to monitoring of the head and neck cancer treatment toxicity. Metabolomics. 2016; https://doi.org/https://doi.org/10.1007/s11306-016-1043-4, 12, 6.
Boguszewicz Ł, Bieleń A, Mrochem-Kwarciak J, Skorupa A, Ciszek M, Heyda A, et al. NMR-based metabolomics in real-time monitoring of treatment induced toxicity and cachexia in head and neck cancer: a method for early detection of high risk patients. Metabolomics. 2019;15(8):110. https://doi.org/10.1007/s11306-019-1576-4.
Alsahafi E, Begg K, Amelio I, Raulf N, Lucarelli P, Sauter T, et al. Clinical update on head and neck cancer: molecular biology and ongoing challenges. Cell Death Dis. 2019;10(8):540. https://doi.org/10.1038/s41419-019-1769-9.
Alexopoulos CG, Pournaras S, Vaslamatzis M, Avgerinos A, Raptis S. Changes in serum lipids and lipoproteins in cancer patients during chemotherapy. Cancer Chemother Pharmacol. 1992;30(5):412–6. https://doi.org/10.1007/BF00689971.
Stathopoulos GP, Stergiou GS, Perrea-Kostarelis DN, Dontas IA, Karamanos BG, Karayiannacos PE. Influence of 5-fluorouracil on serum lipids. Acta Oncol. 1995;34(2):253–6. https://doi.org/10.3109/02841869509093964.
Wang Y, et al. Predictive value of chemotherapy-related high-density lipoprotein cholesterol (HDL) elevation in patients with colorectal cancer receiving adjuvant chemotherapy: an exploratory analysis of 851 cases. Oncotarget. 2016. https://doi.org/10.18632/oncotarget.10145.
Wang G, Su C, Yin T. Paclitaxel and platinum-based chemotherapy results in transient dyslipidemia in cancer patients. Mol Clin Oncol. 2017;6(2):261–5. https://doi.org/10.3892/mco.2016.1107.
Sharma M, Tuaine J, McLaren B, Waters DL, Black K, Jones LM, et al. Chemotherapy agents alter plasma lipids in breast cancer patients and show differential effects on lipid metabolism genes in liver cells. PLoS One. 2016;11(1):e0148049. https://doi.org/10.1371/journal.pone.0148049.
Tian W, Yao Y, Fan G, Zhou Y, Wu M, Xu D, et al. Changes in lipid profiles during and after (neo)adjuvant chemotherapy in women with early-stage breast cancer: a retrospective study. PLoS One. 2019;14(8). https://doi.org/10.1371/journal.pone.0221866.
Poorey VK, Thakur P. Alteration of lipid profile in patients with head and neck malignancy. Indian J Otolaryngol Head Neck Surg. 2016;68(2):135–40. https://doi.org/10.1007/s12070-015-0829-4.
Soininen P, Kangas AJ, Würtz P, Tukiainen T, Tynkkynen T, Laatikainen R, et al. High-throughput serum NMR metabonomics for cost-effective holistic studies on systemic metabolism. Analyst. 2009;134(9):1781. https://doi.org/10.1039/b910205a.
Georges P, Chaussain JL, Gendrel D, Donnadieu M, Job JC. Plasma glucose and alanine relationship after fast in normal children and children with ketotic hypoglycemia. Pediatr Res. 1981;15(1):89. https://doi.org/10.1203/00006450-198101000-00110.
Faerch K, Borch-Johnsen K, Vaag A, Jørgensen T, Witte DR. Sex differences in glucose levels: a consequence of physiology or methodological convenience? The Inter99 study. Diabetologia. 2010;53(5):858–65. https://doi.org/10.1007/s00125-010-1673-4.
Mauvais-Jarvis F. Gender differences in glucose homeostasis and diabetes. Physiol Behav. 2018;187:20–3. https://doi.org/10.1016/j.physbeh.2017.08.016.
Tiziani S, Lopes V, Günther UL. Early stage diagnosis of oral cancer using 1H NMR-based metabolomics. Neoplasia. 2009;11(3):269–IN10. https://doi.org/10.1593/neo.81396.
Banks I. No man's land: men, illness, and the NHS. BMJ. 2001;323(7320):1058–60. https://doi.org/10.1136/bmj.323.7320.1058.
Hoang T, Zhang C, Geye HM, Speer T, Yu M, Wiederholt PA, et al. Gender associations in head and neck squamous cell carcinoma (HNSCC). J Clin Oncol. 2017;30(15_suppl):e16010. https://doi.org/10.1200/jco.2012.30.15_suppl.e16010.
Torri GM, Torri J, Gulian JM, Vion-Dury J, Viout P, Cozzone PJ. Magnetic resonance spectroscopy of serum and acute-phase proteins revisited: a multiparametric statistical analysis of metabolite variations in inflammatory, infectious and miscellaneous diseases. Clin Chim Acta. 1999;279(1-2):77–96. https://doi.org/10.1016/S0009-8981(98)00166-1.
The work has been funded by National Science Centre of Poland (Grant 2015/17/B/NZ5/01387). The funding source had no influence on the design and conduct of the study, analysis and interpretation of the data and preparation of the manuscript.
Ethics approval and consent to participate
All procedures performed in the studies involving human participants were in accordance with the ethical standards of the institutional and/or national research committee and with the 1964 Helsinki declaration and its later amendments or comparable ethical standards. The study was approved by Ethical Committee of the Maria Sklodowska-Curie National Research Institute of Oncology. Informed written consent was obtained from all individual participants included in the study.
Consent for publication
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
NMR pulse sequence parameters; Table S2. - Characteristics of the external test group; Fig. S1. The 400 MHZ 1H NMR CPMG median spectra from the preCHT and postCHT groups; Fig. S2. Result of permutation testing of the obtained OPLS-DA model; Table S3. Radiologically and clinically measured tumor regression.
About this article
Cite this article
Boguszewicz, Ł., Bieleń, A., Jarczewski, J.D. et al. Molecular response to induction chemotherapy and its correlation with treatment outcome in head and neck cancer patients by means of NMR-based metabolomics. BMC Cancer 21, 410 (2021). https://doi.org/10.1186/s12885-021-08137-4
- Head and neck
- Treatment outcome