1. Introduction
Most patients with locally advanced head and neck squamous cell carcinoma (HNSCC) require multimodality treatment [1, 2], frequently consisting of surgery followed by radiotherapy (RT) or primary concurrent radiochemotherapy (RCT) [
3]. While the RT-part directly targets primary tumors and involved cervical lymph nodes, the chemotherapy-part aims at increasing radio-sensitivity and targets circulating tumor cells [2-4].
Modern photon-based radiation aims at primarily damaging primary tumors and suspicious cervical lymph nodes by more focused mechanisms almost exclusively at the site of radiation [5, 6]. Despite the introduction of three-dimensional and, conformal, intensity modulated radiation modalities, radiation-induced changes remains a challenge in RCT of HNSCC-patients [
7]. Thus, direct and desired effects of radiation observed in cancer cells [
6] may also be observed as undesired effects in otherwise healthy tissue, which lies adjacent to but is not infiltrated by HNSCC (e.g. salivary glands, mucosal membranes, or skeletal muscle) [
5]. For skeletal muscle, these alterations are mediated via tissue or stem cell injury, cellular signal pathway alterations, and (epi)genetic changes [8-13]. In contrast, the chemotherapy of a primary RCT was not previously observed to induce similar and significant tissue alterations [3, 4].
The occurrence of these tissue alterations was previously discussed to depend on the interval between the end of RCT and the time of tissue assessment [14, 15]. During the first weeks after RCT, muscular inflammation was discussed to lead to interstitial edema [
16]. In contrast, months after RCT a shift towards muscular fibrosis [14, 15] due misdirected wound healing [17, 18] can be observed. Thus, ultimately functional impairment may be observed [19-21].
Despite these functional impairments [19-21], this temporal dependence of tissue alterations after primary RCT also becomes relevant in the context of salvage surgery [22, 23]. Reports on the incidence of persistent primary tumors or cervical lymph nodes after primary RCT range from 22 up to 40% [24-26]. Irrespective of the risk-benefit ratio of salvage surgery after primary RCT (i.e. missing microscopic residual disease if not performed vs. probability of overtreatment if performed) with possible surgical complications, the optimal time window for salvage surgery is currently considered between 6 and 12 weeks after completion of primary RCT [
27]. For this specific time window, no extensive tissues fibrosis and scarring was postulated, while acute adverse events of primary RCT have already subsided [
27].
Although various experts in the field currently consider this specific time window optimal, the definition of this interval is primarily based on studies exploring overall survival and surgical complications after salvage surgery [
27]. To the best of our knowledge, no direct assessment of alterations of skeletal muscle composition in HNSCC-patients after primary RCT has been performed yet.
Currently, various methods exist to explore RT- and chemotherapy-induced alterations of skeletal muscle composition during and after RCT in HNSCC-patients [
28]. Besides subjective methods such as palpation, which are more of historical relevance, objective methods to assess skeletal muscle composition in the neck were previously performed based on routinely acquired images during routine oncologic follow-up or small trials by using ultrasound or shear wave elastography [17, 18, 21, 29]. Despite these attempts, available data applying these methods to assess RT-induced alteration of skeletal muscle composition during and after RCT of HNSCC-patients are sparse.
No attempt has yet been made to explore short-term RCT-induced alterations of skeletal muscle composition after RCT in HNSCC-patients applying a radiomic approach [16-18, 21, 28, 29]. Radiomics is an emerging data-driven approach aiming at the extraction and processing of quantitative data to analyze image-based information [
20]. The target is to treat medical imaging data of patients as data-minable sources for additional clinical information. The result can provide an ameliorated basis for the decision-making process [
19].
The primary aim of this study was to quantify tissue alterations in skeletal muscles of the head and neck after primary RCT applying radiomic feature analysis on routinely acquired CT images, ultimately aiming at defining an optimal time window for salvage surgery in patients with persistent HNSCC after first line treatment.
2. Materials and Methods
2.1. Study population and additional clinical data
This retrospective cohort study adhered to the “Strengthening the Reporting of Observation studies in Epidemiology” (STROBE) guidelines [
30]. From 2008 until 2021, patients of the institutional head and neck cancer registry at the Department of Otorhinolaryngology, Head and Neck Surgery, Medical University of Innsbruck, Austria, that had 1) incident, histologically confirmed, locally advanced HNSCC (UICC III or IV), 2) who were treated with primary RCT and 3) for whom contrast-enhanced computed tomography scans before treatment (“staging-CT”) and after treatment (“restaging-CT”) were available, were eligible.
From 1,110 potentially eligible patients, 840 did not meet the inclusion criteria. From the remaining 270 patients, 25 have been excluded due to insufficient quality (e.g., dental artifacts) of the contrast-enhanced computed tomography scan (CT) prior to primary RCT (“staging-CT”) and/or after primary RCT (“restaging-CT) (n=10), unavailability of either of the two imaging modalities (n=13) or SCM radiation dose <49 Gy. Of the remaining 245 patients, a representative random sample was drawn with SPSS 28 (IBM, Armonk, NY, USA). The study flow diagram and patient inclusion modified according to Standards for Reporting Diagnostic accuracy studies criteria (STARD) [
31] is depicted in
Figure 1.
2.2. Clinical data
Clinical data has been extracted from the institutional head-and-neck cancer registry or the electronic hospital information system (PowerChart, Cerner, Kansas City, MI, USA) including age, gender, year of first diagnosis, tumor site, UICC stage, date of the CT scans, alcohol consumption, smoking, Radiation dose, p16 status, body mass index (BMI), serum protein level, American Association of Anesthesiologists (ASA) score and functional outcome.
2.3. CT imaging acquisition
Staging and restaging contrast-CTs adhered to the head and neck CT imaging protocols of the Department of Radiology (Medical University of Innsbruck, Austira) and were acquired with a Light Speed VCT or a Light Speed 16 CT scanner (GE Medical Systems, Vienna, Austria). The scan volume ranged from the skull base to the upper mediastinum with a resolution of 512×512 pixels, 2 mm slice thickness, collimation of 24 × 1.2 mm, and 0.45 pitch. Sagittal and coronal images were reconstructed from the axial images. As a contrast agent, Jopamiro 370 (Bracco Austria GmbH, Vienna, Austria) was administered intravenously adjusted to the patient’s body weight. Staging- and restaging contrast CTs used the same protocols, making the studies of each patient comparable.
2.4. Segmentation of head-and-neck musculature
All staging and restaging-CTs were first exported to Digital Imaging and Communications in Medicine format using the IMPAX EE image archiving and communication system (Agfa HealthCare, Bonn, Germany). Thereafter, the images were exported and further processed using mint Lesion™ (Mint Medical GmbH, Heidelberg, Germany, version 3.8.1). For each patient, both sternocleidomastoideus muscles (SCM) and the paravertebral musculature (PVM) including the following muscles: trapezius, longus capitis, splenius capitis, semispinalis capitis, longissimus capitis, levator scapulae, longus colli, were segmented in staging- and restaging-CT.
The SCM, in the field of radiation and additionally affected by scatter-irradiation, was chosen to explore RT-induced alterations, while the PVM, primarily affected by chemotherapy and less by scatter-irradiation, was chosen to explore chemotherapy-induced alterations.
Segmentations were performed manually in all data-sets slice by slice in axial planes using the “paint on slices” tool provided by the software from the upper edge to the lower edge of the 3
rd cervical vertebra. This approach was previously proposed to be effective in the assessment of head-and-neck musculature in HNSCC-patients [
32].
Figure 2.
Example of segmented head-and-neck musculature in a staging-CT of a 56-year-old, male HNSCC-patients with a tumor of the hypopharynx (not depicted) staged cT3 cN0 cM0. The manual slice by slice segmentation was performed in the axial plane using the “paint on slice” tool provided by mint Lesion™ (Mint Medical GmbH, Heidelberg, Germany, version 3.8.1) at the level of the 3
rd cervical vertebra, as previously proposed [
32], for the right SCM (*), the left SCM (**) and the PVM (***), which includes the following muscles: trapezius, longus capitis, splenius capitis, semispinalis capitis, longissimus capitis, levator scapulae, longus colli.
Figure 2.
Example of segmented head-and-neck musculature in a staging-CT of a 56-year-old, male HNSCC-patients with a tumor of the hypopharynx (not depicted) staged cT3 cN0 cM0. The manual slice by slice segmentation was performed in the axial plane using the “paint on slice” tool provided by mint Lesion™ (Mint Medical GmbH, Heidelberg, Germany, version 3.8.1) at the level of the 3
rd cervical vertebra, as previously proposed [
32], for the right SCM (*), the left SCM (**) and the PVM (***), which includes the following muscles: trapezius, longus capitis, splenius capitis, semispinalis capitis, longissimus capitis, levator scapulae, longus colli.
2.5. Data analysis
At the time this study was conducted, mint Lesion™ provided 13 radiomic parameters for the 3 segmented muscle groups SCM right, SCM left, and PVM each before and after therapy. Normally distributed data were described by mean and standard deviation, and non-normally distributed parameters were described by median and 25th and 75th percentiles. Frequencies were tabulated and presented with percentages.
To reduce the number of radiomic parameters (initially 13), a principal component analysis (PCA) with varimax rotation was performed. Here, 3 components could be extracted. Attempts with more components did not yield better results. The 3 extracted components were: 1) volume, 2) mean pixel positivity (MMP) and 3) uniformity. Volume in milliliters (ml) represents the spatial dimension of the segmented muscles SCM and PVM. A change in volume was considered as surrogate for loss of musculature (i.e. sarcopenia) or inflammation and interstitial edema. Mean positivity of pixels (MPP) in Hounsfield-units (HU), represents an intensity parameter, which is a dimensionless absolute number. A change in MPP was considered as surrogate for fibrosis (i.e. increase) or edema (i.e. decrease). Uniformity, a dimensionless absolute number, represents the texture of the segmented muscles SCM and PVM. A change in uniformity was considered as surrogate for changes in texture (i.e fibrosis or edema).
These components had high loadings for the output values of volume, MPP, and uniformity provided by mint Lesion™ (
Table 1; all >0.9.). Therefore, these 3 original values provided by mint Lesion™ were used for the present analyses. The remaining 10 radiomic parameters provided by mint Lesion™ correlated closely with one of these 3 radiomic features each. Thus, they were largely redundant and were not considered in further analysis.
Table 1.
All 13 individual radiomic features and their corresponding radiomic key feature extracted with mint LesionTM (Mint Medical GmbH, Heidelberg, Germany, version 3.8.1).
Table 1.
All 13 individual radiomic features and their corresponding radiomic key feature extracted with mint LesionTM (Mint Medical GmbH, Heidelberg, Germany, version 3.8.1).
Radiomic key features |
Individual radiomic features |
Shape features |
Short axis diameter1 Long axis diameter1 Volume2
|
Texture features3
|
Entropy Uniformity |
Intensity features3
|
Maximal density Minimal density Mean density Skewness of density Standard deviation of density MPP Uniformity of distribution of positive pixels (UPP) Kurtosis |
These three key features represented the spatial dimensions, pixel intensity, and pixel uniformity of the segmented muscles.
The values of the 3 radiomic parameters volume, MPP and uniformity were available for each of the 3 muscle groups (SCM left, SCM right, PVM) before and after therapy. Only the SCM data of the irradiated side were evaluated; if the radiation dose was the same on the right and left sides, the right SCM was used. PVM was omitted from the radiation field and received the lowest possible radiation dose. The radiation dose to the irradiated SCM and the interval between diagnostic CT and restaging-CT were recorded. In addition, data on BMI and serum protein level before and after therapy were available.
First, we tested for linear correlations (Pearson) between values before therapy for the radiomic parameters volume, MPP, uniformity, and patient age, BMI, and serum protein level. The influence of gender, tumor site, UICC-stage, general health as measured by ASA score (I/II vs. III/IV), smoking status (<= or > 10 packyears), and alcohol consumption (daily or less frequently than daily) on pretherapeutic parameters was examined by variance analysis.
The differences before and after therapy for the parameters volume, MPP, uniformity, BMI, and serum protein level were tested with the paired-samples T test (two-sided ). Cohens' d with Hedges correction was used as the effect size parameter. For further mechanistic analyses, the differences after therapy minus before therapy were calculated for these data and subjected to bivariate parametric correlation analyses and variance analysis. In addition, for the radiomic parameters, the results of the paired T tests were adjusted for the effect of BMI difference by including BMI difference as a covariate.
2.6. Ethical considerations
The study has been approved by the review board of the Medical University of Innsbruck, Austria (1269/2018). All procedures conducted in these studies involving human participants are in accordance with the ethical standards of the institutional review board and with the Helsinki declaration (1964) and its later amendments or comparable ethical standards.
3. Results
3.1. Patient population
A total of 1,110 patients recorded in the institutional HNC registry were potentially eligible. Of these, 853 were excluded as they met one or more of the exclusion criteria (
Figure 1). After drawing a representative random sample of the remaining 247 patients, 98 patients were included in this study.
Of these, 73 (74.5%) were male and 25 (25.5%) were female. Mean age at initial diagnosis was 62.0 (± 9.6) years ranging from 42 to 81 years. Of the included 98 patients, the tumor site was oropharynx in 46 (46.9%), hypopharynx in 20 (20.4 %), larynx in 15 (15.3%), oral cavity in 13 (13.3%), and other site in 4 (4.1%) patients. Of the 46 patients diagnosed with oropharyngeal HNSCC, 29 (63.0%) were categorized as p16 positive, using immunohistochemistry with a positivity cut-off of 60%. A total of 93 (94.9%) patients had a UICC Stage III-IV and 5 (5.1%) a UICC stage I-II. HNSCC.
In terms of comorbidities, 53 (54.1%) patients were classified as ASA III/IV and 45 (45.9%) as ASA I/II. All 98 (100.0%) included patients were smokers and 52 patients (53.1%) drank alcohol daily; the remaining 46 patients (46.9%) did not drink daily. BMI data were collected from 79 (80.6%) patients at initial diagnosis. The mean BMI was 24.2 (± 4.9) ranging from 13.8 to 47.7. Thus, 7 (8.9%) patients are underweighted, 41 were (51.9%) normal weighted, 25 (31.6 %) were overweighed and 6 (7.6%) were defined as adipose. Additional detail about clinical characteristics of the 98 included HNSCC-patients is provide in the following table (
Table 2).
Table 2.
Clinical characteristics of the 98 included HNSCC-patients.
Table 2.
Clinical characteristics of the 98 included HNSCC-patients.
|
|
Number |
Percentages |
Sex |
Male |
73 |
74.5% |
|
Female |
25 |
25.5% |
Age |
≤50 |
11 |
11.2% |
|
51-60 |
35 |
35.7% |
|
61-70 |
33 |
33.7% |
|
≥71 |
19 |
19.4% |
Tumor site |
Oral cavity |
13 |
13.3% |
|
Oropharynx |
46 |
46.9% |
|
Hypopharynx |
20 |
20.4% |
|
Larynx |
15 |
15.3% |
|
Others |
4 |
4.1% |
UICC1 truncated |
Stage III |
19 |
19.4% |
|
Stage IV |
79 |
80.6% |
ASA |
ASA I/II |
45 |
45.9% |
|
ASA III/IV |
53 |
54.1% |
Alcohol consumption |
< daily |
52 |
53.1% |
|
daily |
46 |
46.9% |
Smoking habits |
< 10 PY |
24 |
24.5% |
|
≥ 10 PY |
74 |
75.5% |
BMI-classified |
Underweight |
7 |
8.9% |
|
Normal weight |
41 |
51.9% |
|
Overweight |
25 |
31.6% |
|
Adipose |
6 |
7.6% |
Radiation dose1 |
≤60 Gy |
26 |
26.5% |
|
>60 Gy |
72 |
73.5% |
3.2. Primary radiochemotherapy and time intervals
All 98 included HNSCC-patients were treated with primary RCT. Of these, 72 patients (73.5%) received a radiation dosage greater than 60 gray (Gy), 26 patients (26.5%) received 60 Gy or less on the irradiated SCM. Means, standard deviation and range for the imaging interval (i.e. interval between staging-CT and restaging-CT) and restaging interval (i.e. interval between end of primary RCT and restaging-CT) are provided in
Table 3 (
Table 3).
3.3. Variable reduction (principal component analysis)
For dimensional reduction of the radiomic parameters yielded by mint Lesion™, a PCA with varimax rotation was performed. From the 13 parameters of mint Lesion™, 3 principal components could be extracted. The 3 extracted components via PCA explained more than 70% of the variance of the values and were uncorrelated (Appendix 1). The Kaiser-Meyer-Olkin Measure of Sampling Adequacy was 0.61 indicating an acceptable sampling adequacy.
The component matrix allows the interpretation of the 3 extracted components based on their factor loadings (Appendix 2). The first component represents a measure of uniformity of pixels in the segmented muscle, which is inversely correlated with entropy. The second component represents a measure of pixel intensity and the third component represents a measure of spatial dimension.
3.4. Factors influencing pretherapeutic Volume, uniformity and MPP
For comparisons, the right SCM was used as reference for the 3 segmented muscle groups. A correlation (Pearson) between muscle volume, BMI (r= 0.56; p<0.001) and serum protein level (r=0.31; p=0.002) before therapy was observed. In addition, a weak inverse correlation between uniformity and age was observed (r= -0.21; p=0.035).
Further comparisons (independent t-tests) revealed differences in gender. Pretherapeutic SCM volume on the right was 10.0 (±2.8) ml in 73 men and 5.9 (±1.7) ml in 25 women (p<0.001). Differences according to gender were also observed for MPP, 58.8 (±7.9) in men and 63.0 (±9.1) in women (p=0.014), while there were no gender differences for uniformity. The other demographic and clinical factors (UICC stage, general health status (ASA-score), smoking status (<10
Py/>10Py) or alcohol consumption (daily or less frequently) did not affect the radiomic parameters explored.
3.5. Volume, uniformity and MPP before and after therapy
The radiomic parameters for the SCM on the irradiated side of the neck that was exposed to the full radiation dose did not change except for volume (
Table 4). Volume showed a decrease from 9.00 (±3.2) to 8.4 (±2.7) ml (p=0.007), with a Cohen's d of 0.28 indicating a weak effect.
The SCM, in the field of radiation and additionally affected by scatter-irradiation, was chosen to explore RT-induced alterations, while the PVM, primarily affected by chemotherapy and less by scatter-irradiation, was chosen to explore chemotherapy-induced alterations.
The PVM, primarily affected by chemotherapy and less by scatter-irradiation, was chosen to explore primarily for chemotherapy associated muscular changes. However, changes analogous to those in the SCM were observed. The volume decreased from 96.5 (±30.2) to 91.9 (±25.8) (p=0.007; Cohen's d = 0.28). MPP and uniformity showed trends toward increase but may be due to chance (p>0.05,
Table 4). Since the PVM were outside the radiation field, direct radiation exposure cannot explain the observed volume decreases.
Consequently, it was tested whether the volume difference correlates with other parameters. Significant changes before and after therapy were seen in BMI with a decrease from 23.9 (±4.2) to 20.98 (±3.59) kg/m² (p=<0.001; Cohen's d = 0.9) and serum protein levels decreased from 7.4 (±0.53) to 6.75 (±0.85) mg/% (p<0.001; Cohen's d 0.7). Therefore, Pearson correlation analyses were included between volume differences before and after therapy and BMI difference, as well as difference in serum protein levels before and after therapy. First, there was a significant correlation of the volume differences of SCM and PVM (r=0.58; p<0.001). In addition, volume decrease of SCM and decrease of BMI correlated (r=0.27; p=0.007) as well as volume decrease of PVM musculature and BMI decrease (r=0.41; p<0.001). The decrease in serum protein levels did not correlate with the explored radiomic key parameters.
To test whether the radiomic volume decreases were due to the decrease in BMI, the results of the t-tests were adjusted for the difference in BMI values. For this purpose, general linear models completely analogous to the paired t-tests were used, with the BMI difference as a covariate. Here, the irradiated SCM showed a significant interaction between volume and BMI decrease (p=0.007), whereas the difference before/after therapy was no longer significant (p=0.9). Similarly, PVM showed a significant interaction of volume decrease with BMI decrease (p<0.001) and the volume decrease before/after therapy was no longer significant (p=0.71).
4. Discussion
Multimodality treatment is frequently required in the treatment of locally advanced HNSCC [1-4]. Despite the precision of modern RT, which primarily targets the primary tumor and suspect cervical lymph nodes [5, 6], undesired tissue alterations (i.e fibrosis) especially in skeletal muscle of the head and neck were previously observed [14-18]. These alterations occurs in musculature, which was directly in the field of radiation (i.e. SCM) but also in musculature, which was only indirectly affected by scatter irradiation (i.e. PVM ) [14-18].
A link between the observed tissue alterations and the time interval between end of RCT and the time of tissue assessment was previously postulated [14, 15]. This postulated temporal linkage is of crucial importance in the context of salvage surgery, which is frequently necessary in the context of tumor or lymph node persistency after RCT [22-26]. The optimal time window to perform salvage surgery was previously postulated by experts in the field between 6 and 12 weeks after the completion of RCT [
27]. For this specific time window no extensive tissues fibrosis and scarring has to be expected, while acute adverse events of primary RCT have already subsided [
27].
Unfortunately, the definition of this specific time window for salvage surgery was mainly based on studies exploring overall survival and surgical complications after salvage surgery [
27]. Studies that directly assess tissue alterations in skeletal muscle of the head and neck after RCT are sparse [17, 18, 21, 29]. To the best of our knowledge, no attempt has yet been made to explore short-term RCT-induced alterations of skeletal muscle after RCT in HNSCC-patients applying radiomics, a data-driven approach aiming at the extraction and processing of quantitative data to analyze image-based information [19, 20].
The primary aim of this study was to quantify tissue alterations in skeletal muscles of the head and neck after primary RCT applying radiomic feature analysis, ultimately aiming at defining an optimal time window for salvage surgery in patients with persistent HNSCC after first line treatment.
From a total of 247 eligible patients with locally advanced HNSCC recorded in the institutional HNC registry, a representative sample of 98 patients was drawn. Of these 98 patients, common clinical characteristics including sex, age, tumor site, p16 status, smoking and drinking habits were comparable with previous, larger cancer registry based studies [
33] (
Table 2). Thus, the sample drawn from the originally 247 eligible patients appears representative.
The imaging interval (i.e. the time between the pretreatment staging-CT and the posttreatment restaging-CT) was approximately 21 weeks, ranging from 15 weeks to 45 weeks (
Table 3). In this specific interval diagnostic work-up, interdisciplinary tumor board presentation and pretreatment procedures including dental treatments and application of percutaneous gastrostomies were carried out. At our institution this pretreatment procedures prior to the start of primary RCT requires approximately 3 to 4 weeks and primary RCT an additional 6 to 8 weeks. Restaging was performed 6 to 12 weeks after end of RCT. Thus, without complications during pretreatment procedures (e.g. percutaneous gastrostomy wound infections) or during primary RCT (e.g. postponing a radiochemotherapy cycle due to changes in white blood count), this interval ranges from 15 to 24 weeks. In the present study, the maximum interval observed was approximately 45 weeks. In this specific patient a combination of pre- (peritonitis after percutaneous gastrostomy) and intratreament complications (multiple postponing of radiochemotherapy cycles due to neutropenia) occurred.
The restaging interval (i.e. the time between end of the RCT and the restaging-CT) was approximately 8 weeks, ranging from 4 weeks to 33 weeks (
Table 3). Thus, the mean time interval of the population explored in this study is in line with the recommended restaging interval, if additional salvage surgery is required [
27]. The patient with the shortest restaging interval of 29 days was diagnosed with cT4a cN3b cM0 laryngeal cancer and therefore required urgent salvage laryngectomy. Thus a considerably shorter restaging interval was observed for this patient. The patient with the longest restaging interval of 229 days, was the one patient with multiple pre- and intratreatment complications, which required a prolonged intensive care unit stay after completion of the primary RCT.
In summary, the intervals observed in the present study appear in line with intervals reported in previous studies [34-36]. Consequently, the observations of this study may be applied to the postulated optimal time window to perform salvage surgery by experts in the field between 6 and 12 weeks after the completion of RCT [
27].
The first key finding of the present study was that volume was the only radiomic key parameter, which was subject to significant change after primary RCT. The mean SCM-volume decreased from 9.0 to 8.4 ml and the PVM-volume from 96.5 to 91.9 ml (both p=0.007). Some crucial aspects of this observation need to be discussed: considering the restaging interval of approximately 8 weeks (
Table 3), an increase in volume due to muscular inflammation and interstitial edema in the SCM, directly affected by irradiation, was to be expected [
16]. In addition, a significant decrease in mean PVM-volumes was observed, which is primarily affected by chemotherapy but not by irradiation.
Consequently, it was considered unlikely that these observed changes in volume were primarily caused by irradiation, chemotherapy or the combination of the both treatment modalities alone. Moreover additional, significant decreases in mean BMI from 23.9 to 21.0 kg/m² and in mean serum protein levels from 7.4 to 6.6 mg/% (both p<0.001) were observed. These two parameters were included in Pearson correlation analyses, which revealed a significant and strong correlation for the BMI decrease (p<0.001; r=0.41). In a last step, the original results of the t-test applied to the volumetric changes were adjusted for the difference in BMI values ultimately resulting in insignificant volume changes for SCM and PVM after treatment (both p>0.05). Other studies, such as Choi et al., came to a similar conclusion. However, in the aforementioned work, the method of measurement differed and contained significantly fewer parameters [37-40].
Thus, the second key finding of the present study was that no significant changes in the explored radiomic key features was observed, if adjusted for changes in BMI before and after treatment. This observation has several implications.
Firstly, this observation highlights the importance of optimal assessment and optimization of nutrition before, during and after primary RCT. Regular assessments of nutritional status and BMI before, during and after primary RCT for HNSCC-patients was previously recommended by other authors [
41].
Secondly, the data of the present study supports the previously postulated optimal time window to perform salvage surgery [
27]. Neither did the volume of the muscles explored significantly change, if corrected for changes in BMI, nor did the radiomic key features MMP and uniformity, which can be considered as surrogates for tissue fibrosis. Thus, irrespective of the remaining risk-benefit ratio of salvage surgery after primary RCT (i.e. missing microscopic residual disease if not performed vs. probability of overtreatment if performed), the risk of additional surgical complications due to tissue fibrosis or volumetric change in the head and neck musculature is not supported by the present data.
Certain limitations of the present study need to be addressed. Firstly, this comparatively small numbered retrospective study exploring only patients with advanced-stage HNSCC should be supplemented by a larger, prospective investigation. Secondly, various data processing programs are available to segment anatomical structures. For the present study the commercially available software mint Lesion™ was chosen due to the advantage of providing a structured, standardized feature output available to everyone, which minimizes the risk of bias. However, at the time of the study, only 13 radiomic features were routinely extracted by mint Lesion
TM. Furthermore, these 13 features were reduced to three key features via PCA. Although, applying machine learning or deep neural networks for statistical analysis, attempts with more components did not yield better results [
42]. Thirdly, the time interval chosen for the present study only explored for short-term radiomic changes with a mean restaging interval of approximately 8 weeks. An expansion of this observation interval to months or years would be crucial to assess whether previously proposed changes of head and neck musculature due to primary RCT can be detected by the means of a radiomic approach.
Author Contributions
alphabetically ranked by surname: Idea, D.D., G.W. and H.R.; conceptualization, D.D., W.F., H.K., J.M., A.R., M.S.; methodology, D.D., G.W. and H.R.; formal analysis, T.M.G., F.K. P.Z.; investigation, D.D., H.K., H.R., M.S.; data curation, D.D., H.K., H.R., M.S.; writing—original draft preparation, reviewing and editing, D.D., B.H., J.S., H.R., M.S.; supervision, D.D., B.H., J.S., H.R.; project administration, D.D., S.M. and H.R. All authors have read and agreed to the published version of the manuscript.